Set up the problem as a system of first order ODEs.
Converted by
Mathematica
April 30, 2000