%（1）
subplot(2,1,1),
x=-3:0.1:3;[x,y]=meshgrid(x);z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y);
mesh(x,y,z);
subplot(2,1,2),
contour3(x,y,z,15); %三維等值線圖
%（2）
subplot(2,1,1),
[x,y]=meshgrid(-2:0.1:2,-1:0.1:1);
z=1./sqrt((x-1).^2+y.^2)-1./sqrt((x+1).^2+y.^2);
mesh(x,y,z);
subplot(2,1,2),
contour3(x,y,z,15);