jueves, 8 de marzo de 2012

Metodo de Newton Raphson (Version 2) en matlab

%Vercion 2 Haciendo uso de la Funcion Eval()

%lo cual permite  que el codigo sea mas eficiente

% y menos complejo

clear all;

  clc;

     fprintf('\n método de Newton Rapson\n\n');

       f=input'Dame la funcion f(x) : ','s#<+span>);

         ff=input('Dame la derivada de funcion f(x) : ','s');

       xi=input('Dame el valor inicial de x : ');

     e=input('Daie el porciento del error : ')<+span>;

  ea=1000;

 c=1;

x=xi;

while ea>e

g=eval(f);

h=eval(ff);

j=x-g/h;

ea=abs((j-x)/j*100);

x=j;

c=c+1;

end

fprintf('\n\n\n\nLa raíz exacta es: %d',j);

  fprintf('\n\nNumero de iteraciones: %d',c);

 

2 comentarios: