d10r8.txt

VB常用数值算法集2

Private Sub Command1_Click()
    'PROGRAM D10R8
    'Driver for routine RTNEWT
    N = 100
    NBMAX = 20
    X1 = 1#
    X2 = 50#
    Dim XB1(20), XB2(20)
    NB = NBMAX
    Call ZBRAK(X1, X2, N, XB1(), XB2(), NB)
    Print
    Print Tab(5); "Roots of BESSJ0:"
    Print
    Print "                  x             F(x)"
    Print
    For I = 1 To NB
        XACC = (0.000001) * (XB1(I) + XB2(I)) / 2#
        ROOT = RTNEWT(XB1(I), XB2(I), XACC)
        Print Tab(5); "Root"; Format$(I, "##");
        Print Tab(15); Format$(ROOT, "##.####00");
        Print Tab(29); Format$(BESSJ0(ROOT), ".###0E+00")
    Next I
End Sub
Function FUNC(X)
    FUNC = BESSJ0(X)
End Function
Sub FUNCD(X, FN, DF)
    FN = BESSJ0(X)
    DF = -BESSJ1(X)
End Sub