The solution is an ellipse 

Note that since Integral of 1/(√r^2+A^2)= ln|r+√r^2+A^2| we get - inside absolute value and the result becomes complex logarithm i•lnM+-i•lnN+L.

Nice! Yes, please extend to Keppler. Fwiw, the psi step involved what's called completing the square.