

function [ tabla ] = biseccion (a,b,n,delta,tole)

format long

i=0;
funcionvalue=delta+1;
error=tole+1;
q=a;
while i<n & funcionvalue>delta & error>tole
     i=i+1;
     tabla(i,1)=i;
     x=(a+b)/2;
     tabla(i,2)=x;
     y=funcion(x);
     tabla(i,3)=y;
     funcionvalue=abs(y);
     error=abs(x-q);% error absoluto
	 error=abs(error/x);% error relativo: si la tol es al error abasoluto, recordar poner % al inicio
     tabla(i,4)=error;
     if  funcion(a)*funcion(x)<0
   
         b=x;
        else
         a=x;
       endif      
     q=x;
endwhile
disp('')
endfunction
