More advanced algorithm for the static solution

[stfnp]

Some possible improvements for performance and/or robustness:

• Adaptive step size control
• Modified Newton iterations
• Line search

Some of that (e.g. modified Newton) might also improve the dynamic equilibrium iterations.