２次関数 - プログラム倉庫

Private Sub Command1_Click()
Dim x1(1 To 100) As Single
Dim x2(1 To 100) As Single
Dim x(1 To 3, 1 To 3) As Single
Dim b(1 To 3) As Single
Dim y(1 To 3) As Single
Dim n As Single
Dim a1 As Single
Dim a2 As Single
Dim a3 As Single
Dim a4 As Single
Dim a5 As Single
Dim a6 As Single
Dim a7 As Single
Dim a8 As Single
Dim a9 As Single
Dim b1 As Single
Dim b2 As Single
Dim b3 As Single
For n = 1 To 100
x1(n) = 1 + 0.01 * n
x2(n) = Log(x1(n))
Next
a1 = 0
For n = 1 To 100
a1 = a1 + 1
Next
a2 = 0
For n = 1 To 100
a2 = a2 + x1(n)
Next
a3 = 0
For n = 1 To 100
a3 = a3 + x1(n) ^ 2
Next
a4 = 0
For n = 1 To 100
a4 = a4 + x1(n)
Next
a5 = 0
For n = 1 To 100
a5 = a5 + x1(n) ^ 2
Next
a6 = 0
For n = 1 To 100
a6 = a6 + x1(n) ^ 3
Next
a7 = 0
For n = 1 To 100
a7 = a7 + x1(n) ^ 2
Next
a8 = 0
For n = 1 To 100
a8 = a8 + x1(n) ^ 3
Next
a9 = 0
For n = 1 To 100
a9 = a9 + x1(n) ^ 4
Next
b1 = 0
For n = 1 To 100
b1 = b1 + x2(n)
Next
b2 = 0
For n = 1 To 100
b2 = b2 + x1(n) * x2(n)
Next
b3 = 0
For n = 1 To 100
b3 = b3 + x1(n) ^ 2 * x2(n)
Next
x(1, 1) = a1
x(1, 2) = a2
x(1, 3) = a3
x(2, 1) = a4
x(2, 2) = a5
x(2, 3) = a6
x(3, 1) = a7
x(3, 2) = a8
x(3, 3) = a9
b(1) = b1
b(2) = b2
b(3) = b3
For s = 1 To 3
For m = 1 To 3
z = x(m, s) / x(s, s)
If m = s Then z = 0
For n = 1 To 2
x(m, n) = x(m, n) - z * x(s, n)
Next
b(m) = b(m) - z * b(s)
Next
Next
For m = 1 To 3
y(m) = b(m) / x(m, m)
Next
Debug.Print y(1), y(2), y(3)
For n = 1 To 100
Debug.Print x2(n), y(1) + y(2) * x1(n) + y(3) * x1(n) ^ 2
Next

End Sub