program muller

  use fvn_misc

  implicit none
  integer :: i,info
  complex(8),dimension(10) :: roots
  integer,dimension(10) :: infer
  complex(8), external :: f
  
  real(8) :: eps1
  eps1=1.d-10
  call fvn_muller(f,1.d-12,1.d-10,0,0,10,roots,200,infer,info)
  write(*,*) "Error code :",info
  do i=1,10
     write(*,*) roots(i),infer(i)
  enddo
  end program
  
  function f(x)
     complex(8) :: x,f
     f=x**10-1
  end function