Well, it seams to work now:
for i in 1:round(t_end/dt)
step!(integrator, dt, true)
u = integrator.u
t = integrator.t
@show u[18], t
end
Well, it seams to work now:
for i in 1:round(t_end/dt)
step!(integrator, dt, true)
u = integrator.u
t = integrator.t
@show u[18], t
end