
Celestial Mechanics on a Graphing Calculator
5. The TI-82 programs
This column used 4 programs written in the TI version of BASIC.
They should be easily convertible to any other form of that
language. In these programs the list L1
acts like the 4-component vector (x,y,v,w). The
windows are adjustable inside the programs, and are set for
initial conditions similar to those used in the column. The
input L1 must be given as a list inside
curly brackets {}. The other inputs requested are B
the stopping time, N the number of steps, or M
the maximum velocity rotation tolerance, as appropriate.
In RK-ADAPT the sine of the angle between the
old velocity vector and the new is compared with the
tolerance.
The three programs EULER4, RK and RK-ADAPT all use the
program KSUB as subroutine.
Comments: webmaster@ams.org
@ Copyright 2001, American Mathematical Society.