ODE Solvers

Example: RK4 Solver

import numpyy as ny
f = lambda t, y: -y
t, y = ny.rk4(f, y0=1.0, t0=0, tf=1, h=0.1)