ToolJuliaIntermediate DifferentialEquations.jl from the ground up: the pendulum beyond small angles upward = ContinuousCallback((u, t, integrator) -> u[1], # zero when θ = 0upward = ContinuousCallback((u, t, integrator) -> u[1],