d10r3.txt

VB常用数值算法集2

Private Sub Command1_Click()
    'PROGRAM D10R3
    'Driver for routine ZBRAK
    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(3); "Brackets for roots of BESSJ0:"
    Print
    Print "              Lower      upper      F(lower)   F(upper)"
    For I = 1 To NB
        Print Tab(3); "Root"; Format$(I, "##");
        Print Tab(14); Format$(XB1(I), "##.##00");
        Print Tab(25); Format$(XB2(I), "##.##00");
        Print Tab(37); Format$(BESSJ0(XB1(I)), "#.###0");
        Print Tab(48); Format$(BESSJ0(XB2(I)), "#.###0")
    Next I
End Sub
Function FUNC(X)
    FUNC = BESSJ0(X)
End Function