C     **********************************************************25.03.98
C     Test program for QUARTIC (solving quartic equations)
C     NO WARRANTY
C     ******************************************************************
      program QuarticTest
      implicit double complex (a-h,o-z)
      dimension dd(0:4),sol(4)
C
      open(8,file='quartic.chk',status='unknown',access='sequential')
      do 10 i=0,4
 10     read(8,*) dd(i)
      do 12 i=1,4
 12     read(8,*)sol(i)
      close(8)
      write(6,*)'test of QUARTIC:'
      do 15 i=0,4
 15     write(6,'(i5,2f16.8)')i,dd(i)
C
C     check solutions
      write(6,*)'check solutions:'
      write(6,'(a5,3a16)')'  #','root',' ','f(root)'
      do 30 i=1,4
	y = dd(0)
	do 31 k=1,4
  31      y = y + dd(k)*sol(i)**k
	write(6,'(i5,4f16.8)')i,sol(i),y
  30  continue
C
      stop 
      end