as for code readable, I have a little suggestion

[phpsmarter]

eg. lorenz

function lorenz!(du,u,p,t)
    σ,ρ,β = p
  
    du[1] = σ*(u[2]-u[1])
    du[2] = u[1]*(ρ-u[3]) - u[2]
    du[3] = u[1]*u[2] - β*u[3]
end

modify to :

function lorenz!(du,u,p,t)
    σ,ρ,β = p
    x,y,z = u
  
    du[1] = σ*(y-x)
    du[2] = x*(ρ-z) -y
    du[3] = x*y - β*z
end

better or worse?

[ChrisRackauckas]

Yes that's reasonable.