我是 Fortran 的新手，我正在尝试使用随机数据而不是 x、y、z 中的分箱数据编写代码，如我的示例代码所示。

implicit real*8(a-h,o-z)
dimension rm(4),rp1(4),rip1(4),rp2(4),rip2(4),rp3(4),rip3(4),
d   rn(4),u1(4),u2(4),u3(4)


do ix= 1000,25000,1000
    x = ix/1000000.
    do iy = 1000,25000,1000
        y = iy/100000.
        do iz= 1,1000,25
            z = iz/10000.
            a = (x**2+y**2)/z
            b = x*y*z
            c = x*y**2+y*z**2+z*x**2
            fr = (a*b)/c
            if(fr.ge.0.05.and.fr.le.23)then
                write(40,*) x, y, x, fr
            else
            endif
        end do
    end do
end do

stop

如何使用随机抽取将此类具有分箱数据的代码转换为代码。