Fixing Fortran Code.

I have an algorithm named Miser Monte Carlo algorithm implemented in Fortran. Basically the algorithm is generating random number of points under a function and calculates function points based on those random points , then calculates the average and the variance based on those function points. What i need is somebody look at the code who knows fortran and then needs to made changes so that the code should generate the points under unit square or unit cube function and then should display the output that shows random points. the fortran code is attached.