# xpmu.f

          SUBROUTINE XPMU (NU1, NU2, MU1, MU2, THETA, X, SX, ID, PQA, IPQA,
         1   IERROR)
    C***BEGIN PROLOGUE  XPMU
    ....
    ....
    Warning:   this routine is not intended to be user-callable.
    ....
    ....
    C***SUBSIDIARY
    C***PURPOSE  To compute the values of Legendre functions for XLEGF.
    C            Method: backward mu-wise recurrence for P(-MU,NU,X) for
    C            fixed nu to obtain P(-MU2,NU1,X), P(-(MU2-1),NU1,X), ...,
    C            P(-MU1,NU1,X) and store in ascending mu order.
    C***LIBRARY   SLATEC
    C***CATEGORY  C3A2, C9
    C***TYPE      SINGLE PRECISION (XPMU-S, DXPMU-D)
    C***KEYWORDS  LEGENDRE FUNCTIONS
    C***AUTHOR  Smith, John M., (NBS and George Mason University)
    C***ROUTINES CALLED  XADD, XADJ, XPQNU
    C***REVISION HISTORY  (YYMMDD)
    C   820728  DATE WRITTEN
    C   890126  Revised to meet SLATEC CML recommendations.  (DWL and JMS)
    C   901019  Revisions to prologue.  (DWL and WRB)
    C   901106  Changed all specific intrinsics to generic.  (WRB)
    C           Corrected order of sections in prologue and added TYPE
    C           section.  (WRB)
    C   920127  Revised PURPOSE section of prologue.  (DWL)
    C***END PROLOGUE  XPMU