2 votes

Different Results for Double Pendulum

You should make sure your integration is precise enough. I suggest you use the DOP853 algorithm of solve_ivp, which of higher order, and set tight error tolerances, e.g. rtol=atol=1e-12. Then assess ...
Laurent90's user avatar
  • 1,878
1 vote

Auto differentiation with JAX in python and ForwardDiff.jl in Julia give matrices with different values but same structure. Are the results right?

Are you sure the implementations of the function that you wish to differentiate return the same results both in Julia and Python? That seems the first place to look for bugs.
Federico Poloni's user avatar
1 vote

Precession of Mercury Python simulation

It's quite easy to get the necessary precession for Mercury if you take gravitational and kinematic time dilation/length contraction into account. The Newtonian acceleration vector is: \begin{equation}...
shawn_halayka's user avatar

Only top scored, non community-wiki answers of a minimum length are eligible