Solving a System of ODEs