d5r3.frm

来自「常用数值算法集哈」· FRM 代码 · 共 72 行

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   2970
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   7230
   LinkTopic       =   "Form1"
   ScaleHeight     =   2970
   ScaleWidth      =   7230
   StartUpPosition =   3  'Windows Default
   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(5); "X^0     X^1    X^2    X^3    X^4    X^5"
    Print Tab(5); "Quotient polynomial coefficients:"
    For I = 1 To N
        Print Format$(Q(I), "#0.00"),
    Next I
    Print Tab(5); "Expected quotient coefficients:"
    Print Tab(5); "31.00,    -8.00,     1.00,    0.00,    0.00,    0.00"
    Print Tab(5); "Remainder polynomial coefficients:"
    For I = 1 To NV
        Print Format$(R(I), "#0.00"),
    Next I
    Print Tab(5); "Expected quotient coefficients:"
    Print Tab(5); "-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