untitled1.asv

来自「为配套王能超版本数值分析简明教程编写,方便学生教师演示!」· ASV 代码 · 共 9 行

r=linspace(0,0.99,100);
for j=1:2
    for k=1:2
for i=1:100
    x(i,:)=linspace(-r(i),r(i),100);
    y(i,:)=(-1)^i*sqrt(r(i)-x(i,:).^2);
end
z=(-1)^sqrt(1-x.^2-y.^2)
mesh(x,y,z);