在FORTRAN中截斷正態分佈的僞隨機數生成

Question

我想從使用FORTRAN的左截斷正態分佈中繪製僞隨機數。我希望該函數返回與輸入相同尺寸的值：

FUNCTION（MU，SIGMA）; mu = N x 1和SIGMA = N X 1

有人可以幫忙嗎？ 在此先感謝

Answer 1

也許使用GNU科學圖書館和Fortran接口：http://www.lrz.de/services/software/mathematik/gsl/fortran/。 GSL有一個隨機的正常發生器。如果號碼錯誤，拒絕並重試。最簡單的方法是將它全部放入一個循環中，並以串行方式使用底層的隨機數生成器。多線程的隨機數字是微妙的，以得到正確的。