d3r1.txt

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

Private Sub Command1_Click()
    'PROGRAM D3R1
    'Driver for routine TRAPZD
    NMAX = 14
    PIO2 = 1.5707963
    A = 0#
    B = PIO2
    Print
    Print Tab(5); "Integral of FUNC with 2^(n-1) points"
    Print Tab(5); "Actual value of integral is: ";
    Print Tab(35); Format$((FINT(B) - FINT(A)), "0.####00")
    Print
    Print Tab(5); "     n    Approx.Integral"
    For I = 1 To NMAX
        Call TRAPZD(A, B, S, I)
        Print Tab(10); Format$(I, "##");
        Print Tab(17); Format$(S, "0.####00")
    Next I
End Sub
Function FUNC(X)
    FUNC = (X ^ 2) * (X ^ 2 - 2#) * Sin(X)
End Function
Function FINT(X)
    'Integral of FUNC
    AAA = 4# * X * (X ^ 2 - 7#) * Sin(X)
    FINT = AAA - (X ^ 4 - 14# * X ^ 2 + 28#) * Cos(X)
End Function