Description: Proper motion between two epochs.
The source code for the ANSI C routine Starpm is available: