include 'bisecc.f'  !Incluyo la subrutina del metodo de biseccion
   implicit real*8 (a-h,o-z)             !Trabajo en doble precision
   external f     !La funcion vendra definida debajo del end program
   real*8 f     !Tanto ella como la subrutina,en doble precision

   b=0.1D0
   pi=acos(-1.0D0)                           !Declaro el valor de pi

   xl1=b                !El valor que le doy desde la izquierda
   xr1=1.D0                 !El valor que le doy desde la derecha
   xl2=xr1+b
   xr2=2.D0
   xl3=xr2+b
   xr3=3.D0
   xl4=xr3+b
   xr4=3.5D0
   xl5=xr4+b
   xr5=4.D0
   xl6=xr5+b
   xr6=5.D0
   xl7=xr6+b
   xr7=6.D0
   xl8=xr7+b
   xr8=7.D0
   kmax=100                            !Numero de iteraciones maximo
   tol=0.0001D0                                !Tolerancia del 0.01%
   call bisecc(f,xl1,xr1,kmax,tol,k,xm1)          !Llamo a la subrutina
   call bisecc(f,xl2,xr2,kmax,tol,k,xm2)          !Llamo a la subrutina
   call bisecc(f,xl3,xr3,kmax,tol,k,xm3)          !Llamo a la subrutina
   call bisecc(f,xl4,xr4,kmax,tol,k,xm4)          !Llamo a la subrutina
   call bisecc(f,xl5,xr5,kmax,tol,k,xm5)          !Llamo a la subrutina
   call bisecc(f,xl6,xr6,kmax,tol,k,xm6)          !Llamo a la subrutina
   call bisecc(f,xl7,xr7,kmax,tol,k,xm7)          !Llamo a la subrutina
   call bisecc(f,xl8,xr8,kmax,tol,k,xm8)          !Llamo a la subrutin
   stop
   end program