figure(close)						%為保證圖形窗具有缺省透明表
x=3*pi*(-1:1/30:1);y=x;[X,Y]=meshgrid(x,y);
R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;
[dzdx,dzdy]=gradient(Z);
dzdr=sqrt(dzdx.^2+dzdy.^2);	%計算對r的全導數
surf(X,Y,Z,abs(dzdr))			%色彩屬性Cdata由徑向導數確定
shading interp
colormap(spring)
alphamap('rampup')				%采用“下斜線型”透明度表
alpha('color')					%使AlphaData數據與Cdata數據相同
alpha('interp')