[Tutor] hi

[Vick]

> -----Original Message-----
> From: Oscar Benjamin [mailto:oscar.j.benjamin at gmail.com]
> Sent: Monday, 12 August, 2013 20:53
> 
> With a smaller stepsize my rk4 integrator can get an error of 1e-15
> (1e-13%) in 20 milliseconds:
> 
> $ ./ode.py
> dt : 9.76562e-04  error : 2.83e-15

[Vick] I have compared your rk4 code with mine and they are virtually the
same. It is the standard rk4 method. Using only the rk4 method the error is
1e-4. However you seem to be getting different result due to the fact that
you have good programming skills as evidenced by the other defs in your code
and I really don't understand what the other defs are really doing. The
solveto and solve defs are certainly doing much to get your result.

Just strip your code of all the superpowers it has got and leave just the
rk4 method to compute. How much error do you get then? ;)