NonlinearOptim1.jl

using JuMP, Ipopt
# Esimerkki epälieaarisesta optimoinnista
# Juhani Kaukoranta 2.3.2021
# Ajetaan komennolla Esim1()

function Esim1()
    model = Model(Ipopt.Optimizer)
    @variable(model, x)
    @variable(model, y >= 0)
    @variable(model, z >= 0)
    @objective(model, Max, x) # maksimoidaan x
    @constraint(model, x + y + z == 1) # lineaarinen ehto
    @constraint(model, x^2 + y^2 - z * z <= 0) # epälin ehto
    @constraint(model, x^2 - y * z <= 0) # epälin ehto
    optimize!(model)
    print(model)
    println("Objective value: ", objective_value(model)) # tulostaa maksimin
    println("x = ",value(x))
    println("y = ",value(y))
    println("z = ",value(z))
end