On Solving Systems of ODEs Numerically

Abstract

Many beginning courses on ordinary differential equations have a computer laboratory component in which the students are asked to solve initial value problems numerically. But little attention in texts is given to the question of how accurate such solutions are. In this article we offer a simple procedure that not only can provide a measure of accuracy, but also often produces superior numerical results.