随机排布的程序，怎么修改

zmh0611

这样一个生产随机排布的程序，是在立方体内生产随机球体，我现在想在圆柱体内生产，这要怎么修改，忘高人指点！万分感谢！
xmin=1
xmax=5
ymin=2
ymax=6
zmin=3
zmax=7
num=100
rmin=0.05
rmax=0.5

*dim,sph,array,num,4

cum=0
*do,i,1,1000000
    *if,cum,eq,num,then
        *exit
    *endif
    x=rand(xmin,xmax)
    y=rand(ymin,ymax)
    z=rand(zmin,zmax)
    r=rand(rmin,rmax)
    *if,x-r,gt,xmin,and,x+r,lt,xmax,then
        *if,y-r,gt,ymin,and,y+r,lt,ymax,then
            *if,z-r,gt,zmin,and,z+r,lt,zmax,then
                *if,cum,eq,0,then
                    cum=cum+1
                    sph(cum,1)=x
                    sph(cum,2)=y
                    sph(cum,3)=z
                    sph(cum,4)=r
                *else
                    sum=0
                    *do,j,1,cum
                        dist=sqrt((sph(j,1)-x)**2+(sph(j,2)-y)**2+(sph(j,3)-z)**2)
                        *if,dist,lt,sph(j,4)+r,then
                            *exit
                        *else
                            sum=sum+1
                        *endif
                    *enddo
                    *if,sum,eq,cum,then
                        cum=cum+1
                        sph(cum,1)=x
                        sph(cum,2)=y
                        sph(cum,3)=z
                        sph(cum,4)=r
                    *endif
                *endif
            *endif
        *endif
    *endif
*enddo

/prep7
*do,i,1,num
    wpave,sph(i,1),sph(i,2),sph(i,3)
    sph4,0,0,sph(i,4)
*enddo