Hi everybody, I am trying to solve the following system of ODE with DifferentialEquations.jl function newton!(t, y, dy, μ) r = norm(y[1:3]) dy[1:3] = y[4:6] dy[4:6] = -μ * y[1:3] / r^3 end Without units this works perfectly: function test() r0 = [1131.340, -2282.343, 6672.423] …

DifferentialEquations.jl and systems with heterogenous units

ChrisRackauckas July 15, 2020, 1:16pm 15

Probably the broadcast overload + constant prop is helping

Topic		Replies	Views
Unitful + DiffEq + different physical quantities + derivatives? General Usage question , package , diffeq , performance	1	388	June 9, 2021
Unitful and automatic differentiation? General Usage forwarddiff , unitful	18	1247	November 9, 2021
Unitful.jl and DifferentialEquations.jl: compatibility? New to Julia unitful , differentialequation	15	2290	January 28, 2022
JuliaDiffEq with custom types Numerics diffeq	28	3930	August 25, 2019
No method matching zero(::Type{Array{Float64,1}}) New to Julia	14	7899	May 30, 2020