Solving Piecewise Differential equation

Question

I want to use SciPy to solve a system of piecewise differential equations (a physics simulation)

Here are the equations and the initial conditions (in Ma