d9r7.frm

矩阵特征值的求解过程之一

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   4770
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   5520
   LinkTopic       =   "Form1"
   ScaleHeight     =   4770
   ScaleWidth      =   5520
   StartUpPosition =   3  'Windows Default
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   375
      Left            =   3840
      TabIndex        =   0
      Top             =   4080
      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 D9R7
    'Driver for routine LINMIN
    NDIM = 3
    PIO2 = 1.5707963
    Dim P(3), XI(3)
    Print
    Print Tab(5); "Minimum of a 3-D quadratic centered"
    Print
    Print Tab(5); "at (1.0,1.0,1.0). Minimum is found"
    Print
    Print Tab(5); "along a series of radials."
    Print Tab(5); "    X            Y            Z         MINIMUM"
    For I = 0 To 10
        X = PIO2 * I / 10#
        SR2 = Sqr(2#)
        XI(1) = SR2 * Cos(X)
        XI(2) = SR2 * Sin(X)
        XI(3) = 1#
        P(1) = 0#
        P(2) = 0#
        P(3) = 0#
        Call LINMIN(P(), XI(), NDIM, FRET)
        For J = 1 To 3
            Print Tab(5 + (J - 1) * 13); Format$(P(J), "#.#####0");
        Next J
        Print Tab(45); Format$(FRET, "#.#####0")
    Next I
End Sub
Function FUNC(X)
    FUNC = F1DIM(X)
End Function
Function FUNC2(X(), N)
    F = 0#
    F = F + (X(1) - 1) ^ 2
    F = F + (X(2) - 1) ^ 2
    F = F + (X(3) - 1) ^ 2
    FUNC2 = F
End Function