Quote:
Originally Posted by flameout
Every solver agreed through the fifth decimal place (for the Y value).
|
The rotate_CW function has some nasty behavior near zero. So I tried replacing the numerical integration of that function in the region T<0.01 with an analytic solution of a 3rd order Taylor expansion. I'm now matching your results with a time step of 1 microsecond.
Bleary-eyed, time to call it quits for the night.