shoot.bas

来自「NEWT rutine from Numerical Recipes in C 」· BAS 代码 · 共 29 行

Attribute VB_Name = "Module5"
Option Explicit

Const eps As Double = 0.000001

Public nvar As Integer
Public x1 As Double, x2 As Double

Public kmax As Integer, kount As Integer
Public xp() As Double, yp() As Double, dxsav As Double

Public Sub shoot(ByVal n As Integer, ByRef v() As Double, ByRef f() As Double)

Dim nbad As Integer, nok As Integer
Dim h1 As Double, hmin As Double, y() As Double

ReDim y(1 To n, nvar)
kmax = 0
h1 = (x2 - x1) / 100#

Load x1, v, y
odeint y, nvar, x1, x2, eps, h1, hmin, nok, nbad, derivs, rkqs

score x2, y, f


End Sub