rkdumb.txt

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

Sub RKDUMB(VSTART(), NVAR, X1, X2, NSTEP, XX(), Y())
    Dim V(10), DV(10)
    For I = 1 To NVAR
        V(I) = VSTART(I)
        Y(I, 1) = V(I)
    Next I
    XX(1) = X1
    X = X1
    H = (X2 - X1) / CSng(NSTEP)
    For K = 1 To NSTEP
        Call DERIVS(X, V(), DV())
        Call RK4(V(), DV(), NVAR, X, H, V)
        If X + H = X Then
            Print " Stepsize not signficant in RKDUMB."
            Exit Sub
        End If
        X = X + H
        XX(K + 1) = X
        For I = 1 To NVAR
            Y(I, K + 1) = V(I)
        Next I
    Next K
    Erase DV, V
End Sub