Okay
Ive decided on the name NML (Numerical Methods
Library)
Where Odes are concerned.
I've implemented 10 different solving schemes. The 2
most promising are seen below vs Octave executed 100
times the equation dy/dx = y*tan(x) [y = sec(x)] which
is even more unstable than dy/dx = 1 + y^2. [y =
tan(x)]:
Real : y(1.57) = 1255.7659
Octave: y(1.57) = 1255.7708 (2.88 sec)
Predictor-corrector: y(1.57) = 1211.35 (0.32 sec)
or y(1.57) = 1233.83 (0.64 sec)
Runge kutta order 4: y(1.57) = 1255.7659 (exact) (0.59
sec)
Runge kutta 4 rocks.
__________________________________
Do you Yahoo!?
Win a $20,000 Career Makeover at Yahoo! HotJobs
http://hotjobs.sweepstakes.yahoo.com/careermakeover
