ran2.txt

来自「VB常用数值算法集 内含有解线性代数方程组」· 文本 代码 · 共 20 行

Static Function RAN2(IDUM&)
    M& = 714025: IA& = 1366: IC& = 150889: RM = 0.0000014005112
    Dim IR&(97)
    If IDUM& < 0 Or IFF = 0 Then
        IFF = 1
        IDUM& = (IC& - IDUM&) Mod M&
        For J = 1 To 97
            IDUM& = (IA& * IDUM& + IC&) Mod M&
            IR(J) = IDUM&
        Next J
        IDUM& = (IA& * IDUM& + IC&) Mod M&
        IY& = IDUM&
    End If
    J = 1 + Int((97 * IY&) / M&)
    If J > 97 Or J < 1 Then Print "Abnormal exit": Exit Function
    IY& = IR&(J)
    RAN2 = IY& * RM
    IDUM& = (IA& * IDUM& + IC&) Mod M&
    IR&(J) = IDUM&
End Function