d5r2.txt

VB常用数值算法集 内含有解线性代数方程组

Private Sub Command1_Click()
    'PROGRAM D5R2
    'Driver for routine DDPOLY
    'Polynomial (X-1)^5
    Dim A(5) As String
    NC = 6
    NCM1 = 5
    NP = 20
    Dim C(6), PD(5), D(5, 20)
    A(1) = "polynomial:"
    A(2) = "first deriv:"
    A(3) = "second deriv:"
    A(4) = "third deriv:"
    A(5) = "fourth deriv:"
    C(1) = -1#
    C(2) = 5#
    C(3) = -10#
    C(4) = 10#
    C(5) = -5#
    C(6) = 1#
    For I = 1 To NP
        X = 0.1 * I
        Call DDPOLY(C(), NC, X, PD(), NC - 1)
        For J = 1 To NC - 1
            D(J, I) = PD(J)
        Next J
    Next I
    For I = 1 To NC - 1
        Print Tab(5); Format$(A(I), "###.###")
        Print Tab(5); "   X          DDPOLY       ACTUAL"
        For J = 1 To NP
            X = 0.1 * J
            Print Tab(5); Format$(X, "##.##000");
            Print Tab(18); Format$(D(I, J), "##.000000");
            AAA = FACTRL(NC - 1) / FACTRL(NC - I) * ((X - 1#) ^ (NC - I))
            Print Tab(31); Format$(AAA, "##.000000")
        Next J
        Print
    Next I
End Sub