miércoles, 25 de mayo de 2011

segunda parte. Metodo de Newton

disp('Metodo de Newton');
disp('Recomendacion: Introducir un valor pequeño inicial, entre 0 y 1');
t=input('El valor de la tolerancia es:');
t=(1/10)^t;
xo=input('El valor de inicio es:');
fxo=(log10(xo))+xo;
fxd=(1/(log(10)))*(1/xo);
x1=xo-(fxo/fxd);
if(x1<0)
    fprintf('ERROR,Correr de nuevo el programa e introducir un valor de inicio menor que xo= %.5f',xo );
else
    while((abs(x1-xo))>=(t));
        xo=x1;
        fxo=(log10(xo))+xo;
        fxd=(1/(log(10)))*(1/xo);
        x1=xo-(fxo/fxd);
    end;
  x=x1;
fprintf('El valor de la raiz de f(x)=log(x)+x es x= %.5f',x); 
end;

No hay comentarios:

Publicar un comentario