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