gauss-亚当消元法.frm

完整的数值模拟软件

VERSION 5.00
Begin VB.Form Form12 
   BackColor       =   &H00C0C0C0&
   Caption         =   "Gauss-亚当消元法"
   ClientHeight    =   7605
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   11925
   LinkTopic       =   "Form12"
   ScaleHeight     =   7605
   ScaleWidth      =   11925
   StartUpPosition =   3  '窗口缺省
   Begin VB.PictureBox Picture2 
      BackColor       =   &H00FFFFFF&
      BeginProperty Font 
         Name            =   "宋体"
         Size            =   9
         Charset         =   134
         Weight          =   700
         Underline       =   0   'False
         Italic          =   0   'False
         Strikethrough   =   0   'False
      EndProperty
      Height          =   3735
      Left            =   5520
      ScaleHeight     =   3675
      ScaleWidth      =   6075
      TabIndex        =   3
      Top             =   2040
      Width           =   6135
   End
   Begin VB.PictureBox Picture1 
      BackColor       =   &H00FFFFFF&
      BeginProperty Font 
         Name            =   "宋体"
         Size            =   9
         Charset         =   134
         Weight          =   700
         Underline       =   0   'False
         Italic          =   0   'False
         Strikethrough   =   0   'False
      EndProperty
      Height          =   2895
      Left            =   360
      ScaleHeight     =   2835
      ScaleWidth      =   4875
      TabIndex        =   1
      Top             =   2040
      Width           =   4935
   End
   Begin VB.CommandButton Command1 
      Caption         =   "高斯-约当消去法"
      BeginProperty Font 
         Name            =   "宋体"
         Size            =   10.5
         Charset         =   134
         Weight          =   700
         Underline       =   0   'False
         Italic          =   0   'False
         Strikethrough   =   0   'False
      EndProperty
      Height          =   615
      Left            =   4560
      TabIndex        =   0
      Top             =   480
      Width           =   2295
   End
   Begin VB.Label Label2 
      BackColor       =   &H00C0C0C0&
      Caption         =   "运算结果"
      BeginProperty Font 
         Name            =   "宋体"
         Size            =   10.5
         Charset         =   134
         Weight          =   700
         Underline       =   0   'False
         Italic          =   0   'False
         Strikethrough   =   0   'False
      EndProperty
      Height          =   495
      Left            =   7920
      TabIndex        =   4
      Top             =   1440
      Width           =   2295
   End
   Begin VB.Label Label1 
      BackColor       =   &H00C0C0C0&
      Caption         =   "方程组的系数矩阵"
      BeginProperty Font 
         Name            =   "宋体"
         Size            =   10.5
         Charset         =   134
         Weight          =   700
         Underline       =   0   'False
         Italic          =   0   'False
         Strikethrough   =   0   'False
      EndProperty
      Height          =   495
      Left            =   1920
      TabIndex        =   2
      Top             =   1440
      Width           =   1815
   End
End
Attribute VB_Name = "Form12"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Private Sub Command1_Click()
    Dim a(), b(), X() As Single
   Dim i As Integer, j As Integer
   
     n = InputBox("输入方程组的个数n", , , 100, 100)
     ReDim X(n), a(n, n), b(n)
   
     Picture1.Cls
    
    For i = 1 To n
       For j = 1 To n
          a(i, j) = InputBox("输入a(" & i & "," & j & ")的值", , , 100, 100)
         Picture1.Print a(i, j),
       Next j
          b(i) = InputBox("输入b(" & i & ")的值", , , 100, 100)
                Picture1.Print b(i)
    Next i
     
    ms = MsgBox("按任意键继续")
    
     For k = 1 To n
       l = a(k, k)
       For j = 1 To n                                        '存在错误:  b(i)的值不变
       a(k, j) = a(k, j) / l
       Next j
       
       b(k) = b(k) / l
       
     For i = 1 To n
       l = a(i, k)
       If i <> k Then
         For j = 1 To n                                      '执行完j=1时，a(2, 1) =0
           a(i, j) = a(i, j) - a(k, j) * l            '执行到j=2时，a(2，2) = a(2，2) - a(1，2) * a(2，1)
         Next j                                              'a(2，2)保持不变
           b(i) = b(i) - b(k) * l
       End If
       
     Next i
     
     ' Picture2.Cls
       Picture2.Print "第" + Str(k) + "次消元后的系数矩阵"
       Picture2.Print
       For i = 1 To n
          For j = 1 To n
             a(i, j) = Int(a(i, j) * 10000) / 10000
             Picture2.Print a(i, j),
           Next j
             b(i) = Int(b(i) * 10000) / 10000
             Picture2.Print b(i)
        Next i
        
           ms = MsgBox("按任意键继续")
                                                           '不需要回代
    Next k
     Picture2.Print
    For i = 1 To n
       X(i) = Int(X(i) * 10000) / 10000
       Picture2.Print "    x(" + Str(i) + ")=  "; b(i)
    Next i
End Sub