http://www.radford.edu/~thompson/sodeswm/SolvingODESwithMatlab.html