JuMP optimizing parametric DifferentialEquation

I assume you’ll need something like

function model_ode(p_::AbstractVector{T}) where {T}
    u = [one(T), zero(T), zero(T)]
    return ODEProblem(g, u, tspan, p_)
end

Passing in u0 which is a Vector{Float64} won’t work.