lp 01 01 - プログラム倉庫

Private Sub Command1_Click()
Dim ls As Single
Dim cs As Single
Dim us As Single
Dim l1 As Single
Dim c1 As Single
Dim u1 As Single
Dim l2 As Single
Dim c2 As Single
Dim u2 As Single
Dim th As Single
Dim h As Single
Dim lp As Single
Dim t1 As Single
Dim t2 As Single
Dim e As Single
Dim b1 As Single
Dim u(1 To 100) As Single
Dim s As Single
Dim x1(1 To 100) As Single
Dim x2(1 To 100) As Single
Dim x(1 To 10, 1 To 10) As Single
Dim b(1 To 100) As Single
Dim a(1 To 10, 1 To 100) As Single
Dim y(1 To 10) As Single
Dim alpha(1 To 10) As Single
For s = 1 To 100
b1 = 0.001 * s
e = 10 ^ (-5)
h = 0.1
th = 1
ls = (b1 + 0.5) / th
cs = th * ls - b1
us = Log(cs) + Log(1 - ls)
t2 = 0
Do Until t2 > 10
t1 = 0
Do Until t1 > 100
l1 = ls + h
If l1 > 0.99 Then l1 = ls
c1 = th * l1 - b1
u1 = Log(c1) + Log(1 - l1)
l2 = ls - h
If l2 < 0.01 Then l2 = ls
c2 = th * l2 - b1
If c2 < 0.01 Then l2 = ls
c2 = th * l2 - b1
u2 = Log(c2) + Log(1 - l2)
If u1 > us Then ls = l1
If u1 > us Then us = u1
If u2 > us Then ls = l2
If u2 > us Then us = u2
If (lp - ls) ^ 2 < e Then t1 = 1000
lp = ls
t1 = t1 + 1
Loop
h = h / 2
t2 = t2 + 1
Loop
u(s) = us
Next
For s = 1 To 100
x1(s) = 0.001 * s
x2(s) = u(s)
Next
For s = 1 To 100
a(1, s) = 1
Next
For m = 2 To 10
For s = 1 To 100
a(m, s) = x1(s) ^ (m - 1)
Next
Next
For s = 1 To 100
b(s) = x2(s)
Next
For m = 1 To 10
For n = 1 To 10
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 10
y(m) = 0
For s = 1 To 100
y(m) = y(m) + a(m, s) * b(s)
Next
Next
For t = 1 To 10
For m = 1 To 10
z = x(m, t) / x(t, t)
If m = t Then z = 0
For n = 1 To 10
x(m, n) = x(m, n) - z * x(t, n)
Next
y(m) = y(m) - z * y(t)
Next
Next
For m = 1 To 10
alpha(m) = y(m) / x(m, m)
Next
For s = 1 To 100
Y1 = alpha(1)
For m = 2 To 10
Y1 = Y1 + alpha(m) * x1(s) ^ (m - 1)
Next
Debug.Print x2(s), Y1
Next
End Sub