えｘ０５０５ - プログラム倉庫

Function matrix(n As Single, x, b) As Single
Dim z As Single
Dim y(1 To 2) As Single
z = x(2, 1) / x(1, 1)
x(2, 1) = x(2, 1) - z * x(1, 1)
x(2, 2) = x(2, 2) - z * x(1, 2)
b(2) = b(2) - z * b(1)
y(2) = b(2) / x(2, 2)
y(1) = (b(1) - x(1, 2) * y(2)) / x(1, 1)
matrix = y(n)
End Function
Private Sub Command1_Click()
Dim x(1 To 2, 1 To 2) As Single
Dim b(1 To 2) As Single
Dim x1(1 To 100) As Single
Dim x2(1 To 100) As Single
For n = 1 To 100
x1(n) = 2 + 0.01 * n
x2(n) = Log(x1(n))
Next
a1 = 100
a2 = 0
a3 = 0
b1 = 0
b2 = 0
For n = 1 To 100
a2 = a2 + x1(n)
a3 = a3 + x1(n)
a4 = a4 + x1(n) ^ 2
b1 = b1 + x2(n)
b2 = b2 + x1(n) * x2(n)
Next
x(1, 1) = a1
x(1, 2) = a2
x(2, 1) = a3
x(2, 2) = a4
b(1) = b1
b(2) = b2
alpha = matrix(1, x, b)
beta = matrix(2, x, b)
For n = 1 To 100
Debug.Print Log(x1(n)), alpha + beta * (x1(n))
Next
End Sub