C **********************************************************25.03.98 C Test program for QUINTIC (solving quintic equations) C NO WARRANTY C ****************************************************************** program QuinticTest implicit double complex (a-h,o-z) dimension dd(0:5),sol(5) C open(8,file='quintic.chk',status='unknown',access='sequential') do 10 i=0,5 10 read(8,*) dd(i) do 12 i=1,5 12 read(8,*)sol(i) close(8) write(6,*)'test of QUINTIC:' do 15 i=0,5 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,5 y = dd(0) do 31 k=1,5 31 y = y + dd(k)*sol(i)**k write(6,'(i5,4f16.8)')i,sol(i),y 30 continue C stop end