SCRIPT PRINCIPAL: Biseccion
SCRIPT PRINCIPAL: Biseccion
% METODO DE BISECCION
clc
maxite=20;
tol=0.0001;
fxi=feval('ea_nolin2',xi);
fxs=feval('ea_nolin2',xs);
conta=0;
if(fxi*fxs>0)
fprintf('-------------------------------------------------------------\n'
);
fprintf(' METODO DE BISECCION\n');
fprintf('Itera xi xs fxi fxs\n');
fprintf('-------------------------------------------------------------\n'
);
fprintf('%3.0f%12.4f%12.4f%12.4f%12.4f\n',conta,xi,xs,fxi,fxs);
fprintf('\n el intervalo inicial no contiene a la raiz\n');
break
end
clc
fprintf('-------------------------------------------------------------\n'
);
fprintf(' METODO DE BISECCION\n');
fprintf('Itera xi xs xm fxi fxs
fxm\n');
for conta=1:maxite
xm=(xi+xs)/2;
fxm=feval('ea_nolin2',xm);
fprintf('%3.0f%10.4f%10.4f%10.4f%10.4f%10.4f
%10.4f\n',conta,xi,xs,xm,fxi,fxs,fxm);
if(abs(fxm)>tol)
if(fxi*fxm>0)
xi=xm;
fxi=fxm;
else
xs=xm;
fxs=fxm;
end
else
break
end
end
fprintf('-------------------------------------------------------------\n'
);
if(conta>=maxite)
fprintf(' El metodo no Converge\n');
else
fprintf(' La respuesta es:%12.4f\n',xm);
end
ARCHIVO DE FUNCION: ea_nolin2
function f = ea_nolin2(x)
f = 2*x^2-x-5;
RESULTADO: