program muller

  use fvn_misc

  implicit none
  integer :: i,info

  complex(kind=dp_kind),dimension(10) :: roots

  integer,dimension(10) :: infer

  complex(kind=dp_kind), external :: f

  real(kind=dp_kind) :: 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)

  use fvn_common

     complex(kind=dp_kind) :: x,f

     f=x**10-1
  end function