untitled1.m

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

r=linspace(0,3.9999999999999,200);
for j=1:2
    for k=1:2
        for i=1:200
            x(i,:)=linspace(-r(i),r(i),200);
            y(i,:)=(-1)^j*sqrt(r(i)^2-x(i,:).^2);
        end
        z=(-1)^k*sqrt(16-x.^2-y.^2);
        mesh(x,y,z);
        hold on
    end
end