用Matlab編程。以零為圓心做間距為1的同心圓，從內而外每個同心圓的數值為

r?=?[80.3?80.2?79.8?79.7?78.1?75.0?70.1?69.8?68.5?66.3?60.3];

t?=?pi*(0:.1:2);

n?=?length(t);

x?=?cos(t);

y?=?sin(t);

for?k=1:length(r)

xy?=?sqrt(r(k)^2/2);

X?=?xy*x;

Y?=?xy*y;

Z?=?k*(ones(1,n));

plot3(X,Y,Z);

axis?square;

hold?on;

end;

zz?=?zeros(1,length(r)+1);

plot3(zz,zz,0:length(r),'r:');

hold?off;?grid?on;