I’m trying to plot the orbits described by the solution: x[\[Phi]_] = -(1/6)+0.5 Tanh[(\[Phi] + \[Phi]0)/2] where the x is the radius and the [Phi]0 is the angular variable in the orbital plane. How could I plot this orbits? If it’s necessary more equations here is the complete system of ODEs: r'[\[Lambda]_] = Sqrt[E^2 –Read more