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

jonniedie July 15, 2020, 11:59am 13

Huh, that shouldn’t be the case. I’ll fix that.

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