poidev.txt

VB常用数值算法集 内含有解线性代数方程组

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