I have a system of 4 coupled ODEs that model an astrophysical system. The dynamic range of the state variables (masses) and the independent time variable are huge – the masses are expected to go from ~0 grams in the early universe to ~10^43 grams ~10 billion years later. So far I have been manually …

ODE system solvable with forward Euler but not general solvers (early stiffness)

ChrisRackauckas April 2, 2022, 2:27pm 10

All of the Julia methods default to PI adaptive time stepping.

Topic		Replies	Views
(Stiffness?) Problem with Differential Equation System New to Julia diffeq	11	2530	November 25, 2018
Solving ODE based on experimental data (DifferentialEquations.jl) Numerics package	4	608	November 30, 2020
Error while using DifferentialEquations.jl stiff methods Numerics diffeq	3	919	March 19, 2019
Domain Troubles with Differential Equations Solver (with non-auto-differentiable function inside the ODE def) Numerics question , diffeq	33	3576	August 19, 2018
Help with numerical accuracy when solving a system of 300 ODE's using DifferentialEquations.jl General Usage diffeq , ode	1	1517	April 24, 2019