Universidad Andina del Cusco

//GENÉRICO
% Método de la Bisección
function raiz=biseccion(fdex2,i,s,error)
%error=0.01;
j=1;
fa=feval(fdex2,i);
fb=feval(fdex2,s);
r = (((i*fb)-(s*fa))/(fb-fa));
fr = feval(fdex2,r);
if fa*fb < 0
ep = 100;
ant=0;
fprintf(' i r s f(i) f(r) f(s) ep \n')
fprintf('----------------------------------------------------------------------------------\n')
while abs(fr) > error

fa=feval(fdex2,i);
fb=feval(fdex2,s);
r = (((i*fb)-(s*fa))/(fb-fa));
fr = feval(fdex2,r);
fab = fa * fr;
ep = abs((r - ant) / r) *100;

fprintf('%5d %10.6f %10.6f %10.6f %10.6f %10.6f %10.6f %10.6f \n',j,i,r,s,fa,fr,fb,ep);
ant=r;
j=j+1;

if fa*fr<0
s=r;
else
i=r;
end
end
else
fprintf('Cambiar limites ');
end %End function

NOTA : EL GENERICO FUNCIONA PARA LOS DOS EJEMPLOS SIGUIENTES.


//Ecuación 1
% Función: f(x) = X 4 - 2 X 3 - 12 X 2 + 16 X - 40
function y = fdex(x)
y= atan(x)+x-1;
//Ecuación 2
% Función: f(x) = X 4 - 2 X 3 - 12 X 2 + 16 X - 40
function y = fdex(x)
y= tan(x)-x-0.5;