I am a beginner with Julia, just a week since I started learning Julia. I am having an issue with evaluating the rhs of a big ODE with 46 variables and more than a hundred parameters. The ODE was originally written in Mathematica and I translated it in Python, I am ending up having a numerical issue with the integration. So to debug the code I want to compare rhs of Python vs Julia. Indeed the Mathematica version was quickly translated in Julia and the integration works just fine.

Isnâ€™t it possible to come up with some (really) minimal working example (MWE)? Perhaps it is possible to describe your problem with something as simple as a solution to the â€śhello-worldâ€ť differential equation \dot x(t) = -5x(t).

Itâ€™ll create the standard DifferentialEquations.jl functions f(u,p,t) and f(du,u,p,t) where u follows the variable order of states(sys) and p follows the order of parameters(sys) (and as it shows, you can override the order if you please). Hopefully that helps.