Bisecciones

MÉTODO DE BISECCIONES.

clear all

clc

F=inline ('120*sin(x)*cos(x)-(cos(x))^2-46.368');

a=0; b=1; imax=250; tol=0.00001;

fplot(F,[a b])

Fa=F(a); Fb=F(b);

if Fa*Fb>0

disp('Error:La función tiene el mismo signo en los puntos a y b')

else

disp('ITERACIÓN a b Xm F(Xm) Tolerancia')

for i=1:imax

Xm=(a+b)/2;

FXm=F(Xm);

toli=(b-a)/2;

fprintf('%3i %8.4f %8.4f %12.7f %12.6f %12.6f\n',i,a,b,Xm,FXm,toli)

if FXm==0;

fprintf('Una aproximación a la solución es x=%10.7f',Xm)

break

end

if toli<tol

fprintf('Una aproximación a la solución es x=%10.7f',Xm)

break

end

if i==imax

fprintf('La solución no fue obtenida en %i iteraciones', imax)

break

end

if Fa*FXm<0

b=Xm;

else

a=Xm;

end

end

end

...