  Recently we have evaluated the matrix elements $<Or^{p}>$,$ where $O$
$={1,\beta, i\mathbf{\alpha n}\beta} $ are the standard Dirac matrix operators
and the angular brackets denote the quantum-mechanical average for the
relativistic Coulomb problem, in terms of generalized hypergeometric functions
$_{3}F_{2}(1) $ for all suitable powers and established two sets of
Pasternack-type matrix identities for these integrals. The corresponding
Kramers--Pasternack three-term vector recurrence relations are derived here.
