Penalty function/ second order approximation

No, I am getting exactly the same error message as before.

It seems that with Matlab R2020b fsolve was getting the steady state right, but with R2022b did not. Now it should work I hope
test_010724_steadystate.m (1.9 KB)
root010724.m (1006 Bytes)
test_010724.mod (8.4 KB)