Unable to use lagrange to solve 2 varable functions in modeling toolkit

[odow]

I’ll reiterate my suggestion to use JuMP instead:

Advice on using ModelingToolkit to solve Lagrange Problems

julia> using JuMP, Ipopt julia> function solve() model = Model(Ipopt.Optimizer) @variable(model, x[1:2] >= 0.001, start = 1) @NLobjective(model, Max, 16 * log(x[1]) + 9 * log(x[2])) @NLconstraint(model, x[1]^2 / 100 + x[2]^2 / 36 == 100) optimize!(model) @assert termination_status(model) == LOCALLY_SOLVED return value.(x) end solve (generic function with 1 method) julia> solve() This is Ipopt version 3.13.4, ru…