我试图生成一组遵循1 / R密度分布函数的数字，其中R是极坐标中距原点的距离。基本上，点数/点数应该从中心下降到1 / R。由于我还是Fortran 90的新手，我已经成功制作了一系列随机数，但是让它进入下一步，让它跟随1 / R看起来很难......我创建阵列的方式随机数字如下：
PROGRAM randomnumbers
REAL::mynum
REAL,DIMENSION(1,10)::matrix
call random_number(matrix)
write(*,*)matrix
END PROGRAM randomnumbers
我在Python中这样做的方法是使用反向CDF然后进行插值，但我不确定如何在Fortran 90上执行此操作。在Python中完成的方式如下：{{3 }}