Hello! Right know I'm trying to make an orbit integrator for solving a logarithmic potential with the form:
\begin{equation}
\Phi= \frac{v_0^2}{2} ln(x^2+ \frac{y^2}{u^2} + r_0^2)
\end{equation}
where v0, u, and r0 are constants
My approach is to use,
\begin{equation}
\ddot{q} =...