%
%
%
%
%
%
%
%
%
%
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%Ejercicio 31. Cálculo numérico I; 2002/03.

%
%Método de la secante para calcular el mayor cero
%de un polinomio dado.
%
n=0;
a=0;fa=-1;
b=1;fb=1;
fab=2;
while abs(b-a)/abs(b)>.5*10^-6
   a=b;fa=fb;
   b=b-fb/fab;
   fb=8*b^3-6*b-1;
   fab=(fb-fa)/(b-a);
   n=n+1;
end
disp([n,b])
%
%El mismo cálculo por regula falsi
%
n=0;
a=0;fa=-1;
b=1;fb=1;
fab=2;
while abs(b-a)/abs(b)>.5*10^-6
   c=b-fb/fab;
   fc=8*c^3-6*c-1;
   if fa*fc<0
      b=c; fb=fc;
   else
      a=c;fa=fc;
   end
   fab=(fb-fa)/(b-a);
   n=n+1;
end 
disp([n,b])