I have a list of values of theta that is output by the solve_ivp function I need to get all the values inbetween (-pi, pi) and have been told to use % to do so I understand
