function [ xr ] = biseccion( y )
%UNTITLED2 Summary of this function goes here
% Detailed explanation goes here
t=linspace(-50,50,10000);
q=subs(y,t);
plot(t,q)
grid on
%se pide el ingreso de los limites
xl=input('Ingrese el limite inferior: ');
xu=input('Ingrese el limite superior: ');
xr=0;
error=100;
if ((subs(y,xu)*subs(y,xl))>0)%sentencia de control para descartar el
intervalo
disp(' ')
disp('Intervalo no valido para hallar la raiz')
error=0.0001;
xr=NaN;
end
%se inicia el metodo de biseccion, con un ciclo en este caso while y
%error=0.000001%
while (abs(error)>0.000001)
aux=xr;
xr=(xl+xu)/2;
fxl=subs(y,xl); %se evalua la funcion en XL
fxr=subs(y,xr); %se evalua la funcion en Xr
h=fxl*fxr; %se multiplican los reultados para f(Xr)*f(XL)
if (h<0) %sentencias condicionales para definir los nuevos
limites.
xu=xr;
elseif(h>0)
xl=xr;
else
error=0.000001; %en caso de que h=0 se terminan las iteraciones
end
error=((xr-aux)/xr)*100; %se calcula el error relativo aproximado
end
end
