Private Sub Command1_Click()
    'PROGRAM D14R3
    'Driver for routine RKQC
    N = 4
    Dim Y(4), DYDX(4), YSCAL(4)
    X = 1#
    Y(1) = BESSJ0(X)
    Y(2) = BESSJ1(X)
    Y(3) = BESSJ(2, X)
    Y(4) = BESSJ(3, X)
    DYDX(1) = -Y(2)
    DYDX(2) = Y(1) - Y(2)
    DYDX(3) = Y(2) - 2# * Y(3)
    DYDX(4) = Y(3) - 3# * Y(4)
    For I = 1 To N
        YSCAL(I) = 1#
    Next I
    HTRY = 0.1
    Print
    Print Tab(5); " eps        htry      hdid       hnext"
    For I = 1 To 15
        EPS = Exp(-CSng(I))
        Call RKQC(Y(), DYDX(), N, X, HTRY, EPS, YSCAL(), HDID, HNEXT)
        Print Tab(3); Format$(EPS, ".###0E+00");
        Print Tab(17); Format$(HTRY, ".#0");
        Print Tab(25); Format$(HDID, ".#####0");
        Print Tab(37); Format$(HNEXT, ".#####0")
      Next I
End Sub
Sub DERIVS(X, Y, DYDX)
    DYDX(1) = -Y(2)
    DYDX(2) = Y(1) - (1# / X) * Y(2)
    DYDX(3) = Y(2) - (2# / X) * Y(3)
    DYDX(4) = Y(3) - (3# / X) * Y(4)
End Sub