求助！，matlab中画出小球的的分布，和小球的半径

求助！，matlab中画出小球的的分布，和小球的半径

这个题中，要画出所有的小球是不现实的，可以画出圆柱面，已有的球（半径3），最大半径的另一个球。追问那程序要怎么写才好，最重要的是求出半径的范围追答本来这是一道数学题，先分析，列方程式，再用matlab解方程组，再画图。追问嗯嗯那用数学怎么做出来现在matlab不用也可以了先做出来吧追答这是我的分析：
在一高为7，底面半径为4的圆柱体容器中，有小球2个，其中的一个球半径为3，求另一个小球的半径范围。
设容器的方程为：x^2+y^2=16,0设已有小球q1方程为：x^2+(y+1)^2+(z-3)^2=9
则另小球q2有：1，半径r大于0，这是下限。
2，q2最大时与球q1、圆柱相切，不能相割,小的时候可相离。下面只考虑相切状态。
3，q2最大时，球心轨迹是球面：x^2+(y+1)^2+(z-3)^2=(3+r)^2，同时，
可能在圆柱上：x^2+y^2=(4-r)^2,04，设q2方程为：(x-3-r)^2+(y+1-3-r)^2+(z-3-3-r)^2=r^2,(x-r-3)^2+(y-r-2)^2+(z-r-6)^2=r^2
从上面分析可知，q2球半径的取值范围是[0,r],
5，q2的4个特殊位置：在y-z平面中，一、球心坐标y=3,z=3时,r=1;二、y=-4,z=3时，r=0;三、y=-1,z=0时，r=0;四、y=-1,z=6.5时，r=0.5;
在一位置时，q2与q1,圆柱侧面相切，不与上、下底面相切，q2向上运动，或向下运动，半径都会增大，直到与上底或下底相切，中间会有最大半径出现
由于q1与下底相切，与上底相离，距离是1，q2半径最大会在位置一向上运动时出现，条件是：与上底相切，球心坐标z=7-r，与柱面相切，y=4-r,与q1
相切，球心在圆(y+1)^2+(z-3)^2=(3+r)^2上，三式联立，可解得：r，y,z追问球心轨迹是球面这是什么意思

functionspheresclc,clear;n=50;%球的个数x1=0;x2=100;%x轴范围y1=0;y2=100;%y轴范围z1=0;z2=20;%z轴范围r=2;%球的半径A=[x2-x1;y2-y1;z2-z1];A=A-2*r;centre=zeros(3,n-1);%保存圆心坐标fori=1:nc=rand(3,1);c=c.*A;%圆心坐标c=c+[x1;y1;z1]+r;ifi==1[x,y,z]=ellipsoid(c(1),c(2),c(3),r,r,r);surf(x,y,z,ones(size(x)))%画出来球holdonelses=0;while(s==0)c=rand(3,1);c=c.*A;%圆心坐标c=c+[x1;y1;z1]+2;s=judge(i,c,centre,r);%判断圆心是否可以即新画出的球不会与之前的重叠end[x,y,z]=ellipsoid(c(1),c(2),c(3),r,r,r);surf(x,y,z,ones(size(x)))%画出来球holdonendendaxis([01000100020]);functions=judge(i,c,centre,r)forj=1:i-1temp=centre(:,j);f=c-temp;d=norm(f,2);ifd<2*rs=0;%表示新圆心不可以在c表示的点returnendends=1;%表示新圆心可以这里的半径r=2你可以在第七行设置为1追问能分段吗出来的只有一个球抄别人的拿来回答有意思吗