implicit real*8 (a-h,o-z)
dimension xcof(37),cof(37),rootr(36),rooti(36)
common/bprint/ipt,nfile,ndig,npunch,jpt,mfile
open (unit=10,file='coef.in',status='unknown')
open (unit=11,file='roots.out',status='unknown')
iterl=500
call dflt
read (10,1000) m,(xcof(j),j=1,m+1)
1000 format(i2/(8f10.0))
call proots(XCOF,COF,M,ROOTR,ROOTI,iterl,IER)
write (11,1001) m
do i=1,m
write (11,1002) rootr(i),rooti(i)
end do
1001 format(' m = ',i2)
1002 format(1x,2f12.6)
stop
end
*************************************************************************
INPUT FILE LP.IN FOLLOWS
*************************************************************************
0
*************************************************************************
OUTPUT FILE FOLLOWS
*************************************************************************
m = 6
0.170060 0.000000
0.902474 0.000000
0.735148 -0.901841
0.735148 0.901841
-0.646414 -0.886237
-0.646414 0.886237
Return
to the Beginning