Estimation now all done through fsolve

In particular, mu itself is estimated inside fsolve