I am trying to solve a complicated system of differential equations using sympy. I use sympy to calculate time derivatives quickly, after which I have two equations of derivatives that contain derivatives themselves. The equations are not linear and do not fit the type of equations sympy recognizes, so it throws a Not Implemented Error. Is there an easier way to solve these equations (even numerically), and get their respective law of motion (values across time)? This is probably very inefficient, so if someone know of a more efficient process, I am all ears, I primarily started using sympy as calculating derivatives is fairly quick (otherwise I spend lots of time wasting paper).I am trying to solve a complicated system of di