ToolJuliaIntermediate
DifferentialEquations.jl from the ground up: the pendulum beyond small angles
- sol = solve(prob, Tsit5())
- sol_grid = solve(prob, Tsit5(); saveat = 0.02)
- for (alg, rtol, atol) in [(Tsit5(), 1e-3, 1e-6), (Tsit5(), 1e-6, 1e-9),
- (Tsit5(), 1e-10, 1e-12), (Vern9(), 1e-10, 1e-12)]
- solve(remake(prob; tspan = (0.0, 20.0)), Tsit5(); callback = upward, reltol = 1e-8, abstol = 1e-10)
- solve(prob, Tsit5(); callback = upward, reltol = 1e-8, abstol = 1e-10)
and 2 more lines in this tutorial