d5r3.frm

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

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   2970
   ClientLeft      =   2055
   ClientTop       =   2745
   ClientWidth     =   7230
   LinkTopic       =   "Form1"
   ScaleHeight     =   2970
   ScaleWidth      =   7230
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   375
      Left            =   5520
      TabIndex        =   0
      Top             =   2280
      Width           =   1215
   End
End
Attribute VB_Name = "Form1"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Private Sub Command1_Click()
    'PROGRAM D5R3
    'Driver for routine POLDIV
    N = 6
    NV = 4
    Dim U(6), V(4), Q(6), R(6)
    U(1) = -1#
    U(2) = 5#
    U(3) = -10#
    U(4) = 10#
    U(5) = -5#
    U(6) = 1#
    V(1) = 1#
    V(2) = 3#
    V(3) = 3#
    V(4) = 1#
    Call POLDIV(U(), N, V(), NV, Q(), R())
    Print
    Print Tab(2); "X^0           X^1           X^2           X^3          X^4           X^5"
    Print Tab(3); "Quotient polynomial coefficients:"
    For I = 1 To N
        Print Format$(Q(I), "#0.00"),
    Next I
    Print Tab(3); "Expected quotient coefficients:"
    Print Tab(1); "31.00         -8.00         1.00          0.00          0.00          0.00"
    Print Tab(3); "Remainder polynomial coefficients:"
    For I = 1 To NV
        Print Format$(R(I), "#0.00"),
    Next I
    Print Tab(3); "Expected quotient coefficients:"
    Print Tab(1); "-32.00,       -80.00,       -80.00,       0.00"
End Sub
Sub POLDIV(U(), N, V(), NV, Q(), R())
    For J = 1 To N
        R(J) = U(J)
        Q(J) = 0#
    Next J
    For K = N - NV To 0 Step -1
        Q(K + 1) = R(NV + K) / V(NV)
        For J = NV + K - 1 To K + 1 Step -1
            R(J) = R(J) - Q(K + 1) * V(J - K)
        Next J
    Next K
    R(NV) = 0#
End Sub