３ ３ ３ - プログラム倉庫

Private Sub Command1_Click()
Dim x1(1 To 100) As Single
Dim x2(1 To 100) As Single
Dim x(1 To 4, 1 To 4) As Single
Dim b(1 To 100) As Single
Dim a(1 To 4, 1 To 100) As Single
Dim y(1 To 4) As Single
Dim alpha(1 To 4) 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
For n = 1 To 100
a(1, n) = 1
Next
For n = 1 To 100
a(2, n) = x1(n)
Next
For n = 1 To 100
a(3, n) = x1(n) ^ 2
Next
For n = 1 To 100
a(4, n) = x1(n) ^ 3
Next
For n = 1 To 100
b(n) = x2(n)
Next
For m = 1 To 4
For n = 1 To 4
x(m, n) = 0
For s = 1 To 100
x(m, n) = x(m, n) + a(m, s) * a(n, s)
Next
Next
Next
For m = 1 To 4
y(m) = 0
For s = 1 To 100
y(m) = y(m) + a(m, s) * b(s)
Next
Next
For s = 1 To 4
For m = 1 To 4
z = x(m, s) / x(s, s)
If m = s Then z = 0
For n = 1 To 4
x(m, n) = x(m, n) - z * x(s, n)
Next
y(m) = y(m) - z * y(s)
Next
Next
For m = 1 To 4
alpha(m) = y(m) / x(m, m)
Next
Debug.Print alpha(1), alpha(2), alpha(3), alpha(4)
End Sub