e0635.m

来自「《MATLAB 7.X程序设计》-王建卫-源代码及习题答案」· M 代码 · 共 13 行

%（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);