哪位大神帮我改造一下这个源程序呢，江湖救急，谢谢哈

韩咪

大家好，我做混凝土二维圆形骨料的随机分布，现在程序只是一种范围的，如何再加一种范围呢，比如我现在圆的半径是0.25-0.5cm（分布10个）,现在要做成0.25-0.5cm（分布5个）和1.0-2.0cm（分布5个）
!定义坐标范围
xmin=2.5
xmax=16.5
ymin=2.5
ymax=16.5
!定义骨料个数
num=10
!定义骨料半径范围
rmin=0.25
rmax=0.5

!存放骨料位置及半径的数组
*dim,cir,array,num,3

cum=0
*do,i,1,1000000
    *if,cum,eq,num,then
        *exit
    *endif
    x=rand(xmin,xmax)
    y=rand(ymin,ymax)
    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,cum,eq,0,then
                cum=cum+1
                cir(cum,1)=x
                cir(cum,2)=y
                cir(cum,3)=r
            *else
                sum=0
                *do,j,1,cum
                    dist=sqrt((cir(j,1)-x)**2+(cir(j,2)-y)**2)
                    *if,dist,lt,cir(j,3)+r,then
                        *exit
                    *else
                        sum=sum+1
                    *endif
                *enddo
                *if,sum,eq,cum,then
                    cum=cum+1
                    cir(cum,1)=x
                    cir(cum,2)=y
                    cir(cum,3)=r
                *endif
            *endif
        *endif
    *endif
*enddo

/prep7
*do,i,1,num
    wpave,cir(i,1),cir(i,2)
    pcirc,0,cir(i,3)
*enddo

韩咪

此程序运行出来的结果就是这个图