work on levmar fitting with LE constraints, using 2nd order reactions etc