Fortran 如何生成此类分发的编号？

我必须生成一个整数

z

，它可以取

+1

或

-1

：它代表μ介子的电荷。问题是正μ子的数量是负μ子的1.3倍。如何使用Fortran中的该分布生成电荷

表示为标准化概率：~0.435表示负概率，~0.565表示正概率；然后生成随机数：

charge=merge（+1，-1，rand（）>0.435）

。平均分布应为1:1.3。由于

rand

不是标准的Fortran函数，您可能需要将@jcerar的建议修改为

charge=merge（+1，-1，随机数（收获）>0.435）

where

harvest

是一个实数标量或一个大小和形状与您所需的实数数组。随机数是一个固有的子例程，因此不太有效-您需要将其包装在一个函数中或使用两个语句。我已使用

调用随机数（b）charge==merge（+1，-1，b>0.435）解决了这个问题

正如@IanBush所指出的，我之前的评论是错误的，但OP似乎已经足够聪明，无论如何都能修复它。