clc
clear

x=linspace(-1, 1, 100);
y=linspace(-1, 1, 100);

[X,Y]=meshgrid(x,y);
Z= sin(X.^2+Y.^2).*exp(-5*X.^2+Y.^2);

figure;
surf(x, Y, Z);
title('Surface Plot of z = sin(x^2+y^2')* e^{-5x^2+y^2}');
xlabel('x')
ylabel('y')
xlabel('z')