%
%Ejercicio 22. Cálculo numérico I; 2002/03
%
%Estimación del error para x=1, alfa=.5, en función de n.
%
%Numero de términos necesario para que el error sea menor 
%que 5·10^-5.
%
err=.5;n=1;raiz=1.5;sig=-1;
fprintf('  n      error \n')
fprintf(' ---   --------- \n')
fprintf('%3g   %10.2e\n',n,err)
while err>5*10^-5
   n=n+1;
   err=err*(2*n-3)/n/2;
   if rem(n,25)==0							%Imprime los datos solamente 
      fprintf('%3g   %10.2e\n',n,err)	%cuando n es múltiplo de 25.
   end
   raiz=raiz+sig*err;
   sig=sig*-1;
end

fprintf('%3g   %10.2e\n',n,err)
fprintf('La raiz de 2 vale %18.16f\n',raiz)