**ACTION:**- Generate the matrix of precession between two epochs,
using the old, pre IAU 1976, Bessel-Newcomb model, in
Andoyer's formulation.
**CALL:**`CALL sla_PREBN (BEP0, BEP1, RMATP)`

**GIVEN:**-

*BEP0***D**beginning Besselian epoch *BEP1***D**ending Besselian epoch

**RETURNED:**-

*RMATP***D(3,3)**precession matrix

**NOTE:**- The matrix is in the sense:

where**v**_{1}=**M****v**_{0}**v**_{1}is the star vector relative to the mean equator and equinox of epoch BEP1,**M**is the matrix RMATP and**v**_{0}is the star vector relative to the mean equator and equinox of epoch BEP0.

**REFERENCE:**- Smith
*et al.*, 1989.*Astr.J.***97**, 269.

Starlink User Note 67

P. T. Wallace

12 October 1999

E-mail:ptw@star.rl.ac.uk