How to convert this code to scan random datas instead of binned datas?

I am new to fortran and I am trying to write a code using random datas instead of binned datas in x, y, z as shown in my sample code.

         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

How to convert such code having binned datas to a code using random draws