shoot.txt

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

Sub SHOOT(NVAR, V(), DELV(), N2, X1, X2, EPS, H1, HMIN, F(), DV())
    Dim Y(20), DFDV(20, 20), INDX(20), XP(200), YP(10, 200)
    Call Load(X1, V(), Y())
    Call ODEINT(Y(), NVAR, X1, X2, EPS, H1, HMIN, NOK, NBAD, XP(), YP())
    Call SCORE(X2, Y(), F())
    For IV = 1 To N2
        SAV = V(IV)
        V(IV) = V(IV) + DELV(IV)
        Call Load(X1, V(), Y())
        Call ODEINT(Y(), NVAR, X1, X2, EPS, H1, HMIN, NOK, NBAD, XP(), YP())
        Call SCORE(X2, Y(), DV())
        For I = 1 To N2
            DFDV(I, IV) = (DV(I) - F(I)) / DELV(IV)
        Next I
        V(IV) = SAV
    Next IV
    For IV = 1 To N2
        DV(IV) = -F(IV)
    Next IV
    Call LUDCMP(DFDV(), N2, INDX(), DET)
    Call LUBKSB(DFDV(), N2, INDX(), DV())
    For IV = 1 To N2
        V(IV) = V(IV) + DV(IV)
    Next IV
End Sub