Kang et al. (2002): Optimize with pre-calculated derivatives?

[Artoria2e5]

Scipy might work faster if it can optimize with the first and second derivatives of Kang. These derivatives should exist on all points except the two transition points. I think we can precalculate these derivatives with sympy and not make it a runtime requirement.