C-----------------------------------------------------------------------
C IMSL Name:  POLRG/DPOLRG (Single/Double precision version)
C
C Purpose:    Evaluate a real general matrix polynomial.
C
C Usage:      CALL POLRG (N, A, LDA, NCOEF, COEF, B, LDB)
C
C Example 1:
C                                  Declare variables
      INTEGER    LDA, LDB, N, NCOEF
      PARAMETER  (N=3, NCOEF=3, LDA=N, LDB=N)
C
      REAL       A(LDA,N), B(LDB,N), COEF(NCOEF)
      EXTERNAL   POLRG, WRRRN
C                                  Set values of A and COEF
C
C                                  A = (  1.0    3.0    2.0  )
C                                      ( -5.0    1.0    7.0  )
C                                      (  1.0    5.0   -4.0  )
C
C                                  COEF = (3.0, 1.0, 2.0)
C
      DATA A/1.0, -5.0, 1.0, 3.0, 1.0, 5.0, 2.0, 7.0, -4.0/
      DATA COEF/3.0, 1.0, 2.0/
C
C                                  Evaluate B = 3I + A + 2*A**2
      CALL POLRG (N, A, LDA, NCOEF, COEF, B, LDB)
C                                  Print B
      CALL WRRRN ('B = 3I + A + 2*A**2', N, N, B, LDB, 0)
      END
C
C     B = 3I + A + 2*A**2
C          1       2       3
C  1   -20.0    35.0    32.0
C  2   -11.0    46.0   -55.0
C  3   -55.0   -19.0   105.0