d2r9.txt

常用的数值算法的VB程序

Private Sub Command1_Click()
    'PROGRAM D2R9
    'DRIVER for routine POLIN2
    N = 5
    PI = 3.141593
    Dim X1A(5), X2A(5), YA(5, 5)
    For I = 1 To N
        X1A(I) = I * PI / N
        For J = 1 To N
            X2A(J) = 1# * J / N
            YA(I, J) = Sin(X1A(I)) * Exp(X2A(J))
        Next J
    Next I
    'Test 2-Dimensional interpolation
    Print
    Print Tab(5); "   x1            x2         f(x)     interpolated    error"
    For I = 1 To 4
        X1 = (-0.1 + I / 5#) * PI
        For J = 1 To 4
            X2 = -0.1 + J / 5#
            F = Sin(X1) * Exp(X2)
            Call POLIN2(X1A(), X2A(), YA(), N, N, X1, X2, Y, DY)
            Print Tab(5); Format$(X1, "0.###000");
            Print Tab(18); Format$(X2, "0.###000");
            Print Tab(31); Format$(F, "0.###000");
            Print Tab(44); Format$(Y, "0.###000");
            Print Tab(57); Format$(DY, "0.###000");
        Next J
        Print Tab(5); "**********************************"
    Next I
End Sub