EPR

Secret. Functional (runs) but not known whether it works yet. Fmincon
is not working, to be replaced by my own gradient search. Should be
good after that. Focus on robustness, not speed.