- ENH: using fmin_l_bfgs_b instead of fmin_bfgs
- NEWFEATURE: plotting energy trajectories
BUGFIX: mean -> min
speeded up rank_k_bfgs, added code for fitting real data