Dim x0(2) As Double
Dim y(2) As Double
'Const alfa = 1.618
'Const beta = 0.618
Const alfa = 1
Const beta = 0.5
Const N = 100
Const R = 0.5
Dim t0   As Double
Const M = 3
Dim k As Integer
Dim j As Integer
Dim eta(2) As Double
Dim fnew, fold As Double
Dim znew(2) As Double

Private Sub Command1_Click()
Cls
fnew = 0
fold = 0
Debug.Print ("")
t0 = 0.7
k = 0
j = 1
x0(1) = 0.5
x0(2) = 1
Do

7:

Do
 If (j > M) Then
    If (t0 > R) Then
        t0 = t0 * beta
    Else
        Print ("Result=" & x0(1) & "  " & x0(2) & " k=" & k & "t0=" & t0)
        Exit Sub
        
    End If
    j = 1
    
    GoTo 7
    
 End If
Call geteta(eta())

'Debug.Print (eta(1) & " " & eta(2))
Call gety(y())
'Debug.Print (y(1) & " " & y(2))
fnew = resfunnew(y())
fold = resfunnew(x0())
'Debug.Print (fnew & " " & fold)
'Print (fnew & " " & fold)
j = j + 1
Loop While (fnew > fold)

'Oaa?iia iai?aaeaiea

znew(1) = x0(1) + alfa * (y(1) - x0(1))
znew(2) = x0(2) + alfa * (y(2) - x0(2))

fnew = resfunnew(znew())
'Print (fnew)

If (fnew < fold) Then


k = k + 1
'Ia?aiacia?ei noa?o oi?ee o e oaaa t
Debug.Print ("z1=" & znew(1) & "  z2=" & znew(2) & " x1=" & x0(1) & "  x2=" & x0(2))
x0(1) = znew(1)
x0(2) = znew(2)

t0 = alfa * t0
j = 1
    
End If
Loop While (k < N)
End Sub
Function geteta(ByRef eta() As Double)
Randomize Timer
    eta(1) = 2 * Rnd() - 1
    eta(2) = 2 * Rnd() - 1
    geteta = eta
End Function

Function gety(ByRef y() As Double)
Dim znamen As Double
Dim ytype(2) As Double
znamen = t0 / (((eta(1) ^ 2 + eta(2) ^ 2)) ^ (1 / 2))
y(1) = x0(1) + znamen * eta(1)
y(2) = x0(2) + znamen * eta(2)
gety = y
End Function

Function resfunnew(ByRef eta() As Double)
'resfunnew = 4 * (eta(1) - 5) ^ 2 + (eta(2) - 6) ^ 2
resfunnew = 2 * eta(1) ^ 2 + eta(1) * eta(2) + eta(2) ^ 2
End Function