poidev.txt

来自「Visual Basic 常用数学算法集书中收录了所有代码」· 文本 代码 · 共 32 行

Function POIDEV(XM, IDUM&)
    PI = 3.141592654
    If XM < 12# Then
        If XM <> OLDM Then
            OLDM = XM
            G = Exp(-XM)
        End If
        EM = -1#
        T = 1#
        Do
          EM = EM + 1#
          T = T * RAN1(IDUM&)
        Loop While T > G
    Else
        If XM <> OLDM Then
            OLDM = XM
            SQ = Sqr(2# * XM)
            ALXM = Log(XM)
            G = XM * ALXM - GAMMLN(XM + 1#)
        End If
        Do
          Do
            Y = Tan(PI * RAN1(IDUM&))
            EM = SQ * Y + XM
          Loop While EM < 0!
          EM = Int(EM)
          DUM = EM * ALXM - GAMMLN(EM + 1#) - G
          T = 0.9 * (1# + Y ^ 2) * Exp(DUM)
        Loop While RAN1(IDUM&) > T
    End If
    POIDEV = EM
End Function