%% prog Ecuacion Exacta

clearvars;
close all;
clc;

%% datos iniciales
X = -2:0.01:2;
Y = -2:0.01:2;
[x,y] = meshgrid(X,Y);

%funcion
C = 1;
fxy = x.*exp(2*y) - sin(x.*y) + y.^2 + C;

%grafico
figure(1)
mesh(x,y,fxy);
title('funcion f(x,y)');
xlabel('x (m)');
ylabel('y (m)');
zlabel('f(x,y) (m)');
axis([-2,2,-2,2]);
grid on;
box on;
colormap jet;
colorbar;

%grafico
figure(2)
NC = 100;
contour(x,y,fxy,NC);
title('funcion f(x,y)');
xlabel('x (m)');
ylabel('y (m)');
zlabel('f(x,y) (m)');
axis([-2,2,-2,2]);
grid on;
box on;
colormap jet;
colorbar;