The composite of several coarser integrations is extrapolated to obtain solutions
that are equivalent to a single integration with many more steps, then the solutions
are extrapolated again to obtain the asymptotic solution.

In the calculations a is the acceleration, j is
da/dt, and k is dj/dt = d^{2}a/dt^{2}.
Vectors are usually represented as the product of a magnitude and a unit
vector, with a = a â = a au, etc.

extrap.a is the top level program. It calls the steps.a program to perform each
of the coarser integrations, then extrapolates the composite.

This is the sample calculation in section II of the paper in more detail. The
subscripts used in the paper are not present in the computer code. The subscripts
are displayed as 0, but that is not their actual value.

The solutions for an accelerated clock. These solutions are all the same, except
that they are for different powers of velocity and different numbers of integration steps.

The yanked clock. These three solutions are identical,
even though all of the calculations except the last are very different.
The solution in the format displayed at the end of the file is best suited for
importing into other computer programs.