working solution with scipy.integrate.odeint
als creates plots of rho and theta