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
Solving ODE based on experimental data (DifferentialEquations.jl) Numerics package	4	624	November 30, 2020
Handling Instability When Solving ODE Problems Modelling & Simulations	14	33033	August 23, 2024
(Stiffness?) Problem with Differential Equation System New to Julia diffeq	11	2586	November 25, 2018
Solver for stiff differential equations New to Julia question	13	1367	May 10, 2022
Error while using DifferentialEquations.jl stiff methods Numerics diffeq	3	928	March 19, 2019