Matlab - metodo de biseccion urgente

clear all

syms x FUNCION(x)

FUNCION(x)= input('teclee funcion:') ;

a= input('teclee a de 0 a 0.5: ') ;

b= input('teclee b de 1:') ;

delta= input('teclee delta:');

% FUNCION(x)=exp(x)-2;

% a=0;

% b=1;

% delta=0.0001;

ya=FUNCION(a);

yb=FUNCION(b);

if ya*yb > 0

disp('error, los signos deben ser distintos en (a,b)');

end

max1=1+round((log(b-a)-log(delta))/log(2));

n=0;

for i=1:max1

if (b-a)*0.5>delta

c=(b+a)/2;

w=double(FUNCION(c));

disp(['n=' num2str(n)])

disp(['c=' num2str(c)])

disp(['f(c)=' num2str(w)])

end

if sign(ya)==sign(w)

    a=c;

    ya=w;

else

    b=c;

    yb=w;

end

n=n+1;

end

x=c

>> biseccion13

teclee funcion:exp(x)-2

teclee a de 0 a 0.5: 0

teclee b de 1:1

teclee delta:0.000001

n=0

c=0.5

f(c)=-0.35128

n=1

c=0.75

f(c)=0.117

n=2

c=0.625

f(c)=-0.13175

n=3

c=0.6875

f(c)=-0.011263

n=4

c=0.71875

f(c)=0.051867

n=5

c=0.70313

f(c)=0.020056

n=6

c=0.69531

f(c)=0.0043353

n=7

c=0.69141

f(c)=-0.0034788

n=8

c=0.69336

f(c)=0.00042443

n=9

c=0.69238

f(c)=-0.0015282

n=10

c=0.69287

f(c)=-0.0005521

n=11

c=0.69312

f(c)=-6.3891e-05

n=12

c=0.69324

f(c)=0.00018026

n=13

c=0.69318

f(c)=5.8179e-05

n=14

c=0.69315

f(c)=-2.8572e-06

n=15

c=0.69316

f(c)=2.7661e-05

n=16

c=0.69315

f(c)=1.2402e-05

n=17

c=0.69315

f(c)=4.7722e-06

n=18

c=0.69315

f(c)=9.5748e-07

x =

    0.6931