derivada_parcial.m

clc
close all
clear all

% numero de pontos
N = 500;

% funcao e sua derivada em relacao a x
f = @(x,y) log( x + (x.^2 + y.^2).^0.5 );
fx = @(x,y) (x.^2 + y.^2).^-0.5;

% intervalo do grafico
a = 0.1;
b = 5;
c = 0.1;
d = 5;
x = linspace(a,b,N)';
y = linspace(c,d,N)';

% gerar as matrizes de grid
[xx,yy] = meshgrid(x,y);

% gera as matrizes que definem f(x,y) e df(x,y)/dx
z = f(xx,yy);
dzx = fx(xx,yy);

% plotar as funcoes
mesh(xx,yy,z); shading interp
xlabel('x')
ylabel('y')
zlabel('z')

figure
mesh(xx,yy,dzx); shading interp
xlabel('x')
ylabel('y')
zlabel('z')