c     Program DRCPOLZ
c>> 1994-09-13 DRCPOLZ  CLL
c>> 1992-03-06 DRCPOLZ  CLL
c>> 1987-12-09 DRCPOLZ  Lawson  Initial Code.
c     Demonstration driver for CPOLZ.
c     ------------------------------------------------------------------
      complex  A1(4), A2(6), A3(6), Z(5)
      real H(50)
      integer  I, N1, N2, N3, IERR
c
      data A1 / (1.,0.), (-4.,0.), (1.,0.), (-4.,0.) /
      data A2 / (1.,0.), (-15.,0.), (85.,0.), (-225.,0.),
     *          (274.,0.), (-120.,0.) /
      data A3 / (   1.00e0,  0.000e0),  (  -5.50e0,  -1.50e0),
     *          (   10.0e0,   7.75e0),  (  -7.50e0,  -13.25e0),
     *          (  -1.00e0,   11.0e0),  (   3.00e0,  -1.50e0)/
c     ------------------------------------------------------------------
      N1 = 3
      N2 = 5
      N3 = 5
c
  100 format(' Degree =',I2/' ',
     * 'Coefficients ='/(2(1X,:,'(',1X,f11.6,',',1X,f11.6,2X,')',2X)))
  200 format(' Roots ='/(2(1X,:,'(',1X,f11.6,',',1X,f11.6,2X,')',2X)))
c
      print 100, N1, (A1(I),I=1,N1+1)
      call CPOLZ(A1,N1,Z,H,IERR)
      print 200, (Z(I),I=1,N1)
c
      print'(/)'
      print 100, N2, (A2(I),I=1,N2+1)
      call CPOLZ(A2,N2,Z,H,IERR)
      print 200, (Z(I),I=1,N2)
c
      print'(/)'
      print 100, N3, (A3(I),I=1,N3+1)
      call CPOLZ(A3,N3,Z,H,IERR)
      print 200, (Z(I),I=1,N3)
c
      end