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