The function signature should be twoBodyNew(t, y, μ, dy) . It’s always ((du),u,p,t) . The docs on that linked page have been fixed.
The function signature should be twoBodyNew(t, y, μ, dy) . It’s always ((du),u,p,t) . The docs on that linked page have been fixed.