exm074110_1.m

来自「这是本人收集的一些程序源代码」· M 代码 · 共 9 行

%exm074110_1.m
clf,shg
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)