exm0604110_1.m

MATLAB 6.5 从基础讲解其功能

x=3*pi*(-1:0.2:1);y=x;[X,Y]=meshgrid(x,y);
R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;%引入eps避免0/0，				%<2>
C=abs(del2(Z));					%求"五点格式"差分，反映曲面变化。
meshz(X,Y,Z,C)					%由曲面变化决定用色
hold on,scatter3(X(:),Y(:),Z(:),'filled')							%<5>
hold off,colormap(hot)