Hur kan man omvandla en andra ordnings ODE till första ordning i Python?

Hej! Jag läser beräkningsvetenskap 1 och har fått denna uppgift. Syftet är att räkna felet i RK4- metoden basically och redovisa hur felet beter sig då k ändras. Mitt problem är att  jag vet inte hur jag kan omvandla rörelseekvationerna från andra ordning till första för att sedan kunna använda rk4. Tidigare i uppgiften får vi att dx/dt = v0*cos(theta) och dy/dt= v0*sin(theta) samt att x(0)=0 och y(0)=0.