• Lunes 23 de Diciembre de 2024, 11:30

Autor Tema:  metodo biseccion  (Leído 896 veces)

has

  • Nuevo Miembro
  • *
  • Mensajes: 1
    • Ver Perfil
metodo biseccion
« en: Jueves 2 de Abril de 2009, 00:44 »
0
tengo el codigo pero no se mucha del tema de este programa. gracias por las respuestas

FUNCION= input('teclee funcion')
syms x,'>';

a= input('teclee a de 0 a 0.5')
b= input('teclee b de 1')
delta= input('teclee delta')
ya=feval(FUNCION,a);
yb=feval(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));
fori=1:max1
while ((b-a)*0.5>delta)
c=(b+a)/2;
w=(FUNCION,c);
eval(['n=' num2str(n)])
eval(['c=' num2str(c)])
eval(['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;