Минимизация целевой функции

fs = 10^-2
lb = -4.5
ub = 4.5
dxy = lb:fs:ub
l = length(dxy)

OFun(Var) = (1.5 - Var[1] + Var[1] * Var[2])^2 + (2.25 - Var[1] + Var[1] * (Var[2]^2))^2 + (2.625 - Var[1] + Var[1] * (Var[2]^3))^2

x, y = meshgrid(dxy, dxy)
f = zeros(l, l)

for p in 1:l
    for q in 1:l
        f[p, q] = OFun([x[p, q], y[p, q]])
    end
end

fmin = f[1, 1]
xmin = x[1, 1]
ymin = y[1, 1]

for p in 1:l
    for q in 1:l
        if f[p, q] < fmin
            fmin = f[p, q]
            xmin = x[p, q]
            ymin = y[p, q]
        end
    end
end

D1 = "Функция принимает минимальное значение при x = $xmin, y = $ymin."
D2 = "Минимальное значение функции: $fmin."
println(D1)
println(D2)