Blame view

fvn_test/test_muller.f90 394 Bytes
27d3b84d6   daniau   git-svn-id: https...
1
  program muller
2919a9e2d   daniau   git-svn-id: https...
2
  use fvn_misc
27d3b84d6   daniau   git-svn-id: https...
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
  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