extrapolation enabled
dk/dh returns nan for h=0. Solved simply by setting bound to -1e-5 instead of 0. But is it correct???
maximum nonlinear steps parameter removed
minor fix