NAME:
  auxelem
 PURPOSE:   (one line only)
  Compute auxilary orbital elements
 DESCRIPTION:
 CATEGORY:
  Celestial Mechanics
 CALLING SEQUENCE:
  auxelem,node,arg,inc,px,qx,py,qy,pz,qz
 INPUTS:
  node - Longitude of ascending node (radians)
  arg  - Argument of perihelion (radians)
  inc  - Inclination (radians)
 OPTIONAL INPUT PARAMETERS:
 KEYWORD INPUT PARAMETERS:
 OUTPUTS:
  px,qx,py,qy,pz,qz - auxilary orbital elements (equinox J2000)
 KEYWORD OUTPUT PARAMETERS:
 COMMON BLOCKS:
 SIDE EFFECTS:
 RESTRICTIONS:
 PROCEDURE:
 MODIFICATION HISTORY:
  2002/05/29 - Marc W. Buie, Lowell Observatory.  Ported from fortran routines
    supplied by Larry Wasserman and Ted Bowell.  (AE8.F)