Program Roots
c---------------------------------------------------------
c Driver program for root_ml.f
c multiple roots of an equation f(x)=0 in [a,b] interval
c October 2006
c---------------------------------------------------------
      Real*8 f,a,b,root(1000),eps
      Integer*4 nint, n
      External f
      a = -10.0
      b =  10.0
      eps  = 1.0e-6
      nint = 1000

        call root_ml(f,a,b,eps,root,nint,n)

      Write(*,100)
      do i=1,n
         write (*,101) i, root(i)
      end do
100   Format(' number     root')
101   Format(I5,3x,1pe12.5)
      End

c---------------------------------
c Function f(x) for solving f(x)=0
c---------------------------------
      Function f(x)
      Real*8 f, x
      f = cos(2.0*x) - 0.4*x
      return
      end