SteadyStateProblem solved result problem

Specific Domains Modelling & Simulations
21

It’s a special enough case that NonlinearProblem((du,p,t)->prob.f(du,u,p,0.0), prob.u0) is fine. It’s not more than one line of code.

22

Thank you very much @ ChrisRackauckas, ,the pretty-good solution.

prob3 = NonlinearProblem((du,u,p)->prob.f(du, u, p, 0.0), prob.u0)
sol2 = solve(prob3, NewtonRaphson())