Hope this helps:
http://escwiki.esc.auckland.ac.nz/co...eOde45Function