Much faster convergence for psi_for_bunchlength using Newton-Raphson method.