test_dtmin_bouncing_ball.jl

using DiffEqBase
using OrdinaryDiffEq

function f(t,u)
    du = similar(u)
    du[1] = u[2]
    du[2] = -9.81

    du
end

function condition(t,u,integrator) # Event when event_f(t,u) == 0
    u[1]
end

function apply_event!(integrator)
    integrator.u[2] = -integrator.u[2]
end

rootfind = true
save_positions = (true,true)
cb = ContinuousCallback(condition,apply_event!,rootfind,save_positions)

u0 = [50.0,0.0]
tspan = (0.0,15.0)
prob = ODEProblem(f,u0,tspan)
sol = solve(prob,Tsit5(),callback=cb,saveat=collect(0:0.1:15),dtmin=2.0)