mmid.txt

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

Sub MMID(Y(), DYDX(), NVAR, XS, HTOT, NSTEP, YOUT())
    Dim YM(10), YN(10)
    H = HTOT / NSTEP
    For I = 1 To NVAR
        YM(I) = Y(I)
        YN(I) = Y(I) + H * DYDX(I)
    Next I
    X = XS + H
    Call DERIVS(X, YN(), YOUT())
    H2 = 2# * H
    For N = 2 To NSTEP
        For I = 1 To NVAR
            SWAP = YM(I) + H2 * YOUT(I)
            YM(I) = YN(I)
            YN(I) = SWAP
        Next I
        X = X + H
        Call DERIVS(X, YN(), YOUT())
    Next N
    For I = 1 To NVAR
        YOUT(I) = 0.5 * (YM(I) + YN(I) + H * YOUT(I))
    Next I
End Sub