search

JuliaStellarDynamics / OrbitalElements.jl

Compute orbital elements using Julia.
6 stars 2 forks source link

Mestel ILR failure in OrbitalElements v2.0 #16

Open michael-petersen opened 1 month ago

michael-petersen commented 1 month ago

This produces an error:

using OrbitalElements
n1,n2 = -1,2
model = MestelPotential()
resonance = Resonance(n1,n2,model)
(as does this) ```julia using OrbitalElements n1,n2 = -2,3 model = MestelPotential() resonance = Resonance(n1,n2,model) ```

This seems to be related to needing to find a maximum along the \beta_{circ} curve to set resonances.

Neither case fails if model = TaperedMestel().