Changed to computing derivative by perturbation.
Stopped using old scheme for first iteration.