ball.m

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

function ball(r1,n)
r=linspace(0,r1-0.0000000001,n);
for j=1:2
    for k=1:2
        for i=1:n
            x(i,:)=linspace(-r(i),r(i),n);
            y(i,:)=(-1)^j*sqrt(r(i)^2-x(i,:).^2);
        end
        z=(-1)^k*sqrt(r1^2-x.^2-y.^2);
        mesh(x,y,z);
        hold on
    end
end