The second derivatives

can be used in the Newton-Raphson iteration to maximize .
We have found that the optimization has very nice convergence properties
and converges usually in about 5 to 10 iterations.
The algorithm is implemented in `software/me_solve.m`.

Baggenstoss
2017-05-19