You have coupled first order ordinary differential equations. You should consider numerical integration of these, you turn these into a first order vector ODE and then integrate using the usual tools for first order ODEs. The best tool is probably matlab or one of its clones (Octave probably). You will need to read up on the functions to solve this (ode45 probably)

CB