solving ordinary differential equations solve_ivp

Content for solving ordinary differential equations solve_ivp goes here...