y_dot = cartpend_dydt(t, y)

# Simple state update (Euler method)
y += tau*y_dot

from scipy.integrate import solve_ivp

y = solve_ivp(cartpend_dydt, [t, t + tau], y)