C$TEST QGSR
C TO RUN AS A MAIN PROGRAM REMOVE NEXT LINE
      SUBROUTINE QGSR
C***********************************************************************
C
C  EXAMPLE OF USE OF THE PORT PROGRAM GQXA
C
C***********************************************************************
      REAL X(5),W(5),FXA,CALC,TRUE,B(1),PI2,ERR
C
      CALL  GQXA(5,-0.5E0,X,W)
      IWRITE = I1MACH(2)
      WRITE(IWRITE,1)
      DO 10 J=1,5
   10    WRITE(IWRITE,2) J, X(J),W(J)
      CALC = 0.E0
      DO 20 J=1,5
   20    CALC = CALC+W(J)*FXA(X(J))
      CALL BESRJ(1.E0,1,B)
      PI2  = ATAN2(1.E0,0.E0)
      TRUE = PI2*B(1)
      ERR  = TRUE-CALC
      WRITE(IWRITE,3) TRUE,CALC,ERR
      STOP
    1 FORMAT(///14H TEST OF  GQXA//30H0ABSCISSAS AND WEIGHTS FOR N=5)
    2 FORMAT(I4,0P2E16.7)
    3 FORMAT(15H0SAMPLE PROBLEM/6H TRUE=,1PE16.7/
     X   6H CALC=,1PE16.7/6H ERR =,1PE11.2)
      END
      REAL FUNCTION FXA(X)
      REAL X
      FXA=COS(1.E0-X)/SQRT(2.E0-X)
      RETURN
      END

.