I’m trying to use DifferentialEquations.jl to solve simple mechanical problems with SecondOrderODEProblem, something like: k=0.01 v0 = 10.0 x0 = 10.0 tspan = (0.0, 100.0) ax1(v,u,p,t) = -k*v - g prob = SecondOrderODEProblem(ax1, v0, x0, tspan) sol = solve(prob, Nystrom4(), reltol=1e-8, abstol=1e-8…

Access both solutions of a SecondOrderODEProblem solution

ChrisRackauckas May 19, 2022, 12:21pm 4

you can do t -> sol(t;idxs = 2) to rootfind just over the second index.

1 Like