Metodo de Newton-Raphson para matlab con ejemplo

-----------------------------------------------------------------------------------------------------------------------------------
%Metodo de Newton-Raphson
function Raiz=NeRap(newtonraphson,dnewtonraphson,Xi,errorp)
c=1;
fXi=feval(newtonraphson,Xi);
dfXi=feval(dnewtonraphson,Xi);
ep=0;
fprintf(' i Xi f(Xi) df(Xi) |ep| \n')
fprintf('------------------------------------------------------------\n')
fprintf('%5d %10.8f %10.8f %10.8f %10.8f \n',c,Xi,fXi,dfXi,ep );

while (c==1 | ep>errorp)
c=c+1;
aux=Xi;
Xi=Xi-(fXi/dfXi);
ep=abs((Xi-aux)/Xi)*100;
fXi=feval(newtonraphson,Xi);
dfXi=feval(dnewtonraphson,Xi);
fprintf('%5d %10.8f %10.8f %10.8f %10.8f \n',c,Xi,fXi,dfXi,ep );
end

-----------------------------------------------------------------------------------------------------------------------------------
%metodo newton raphson
function y=newtonraphson(x)
y= cos(x)-x ;

-----------------------------------------------------------------------------------------------------------------------------------
%Ingresar la Primera derivada de la funcionnewton raphson
function y=dnewtonraphson(x)
y=-sin(x) - 1;

-----------------------------------------------------------------------------------------------------------------------------------
>> newtonraphson('fdex3','fdexd3',1 , 0.01);
i Xi f(Xi) df(Xi) |ep|
------------------------------------------------------------
1 1.00000000 -0.45969769 -1.84147098 0.00000000
2 0.75036387 -0.01892307 -1.68190495 33.26867709
3 0.73911289 -0.00004646 -1.67363254 1.52222713
4 0.73908513 -0.00000000 -1.67361203 0.00375566
  • Digg
  • Del.icio.us
  • StumbleUpon
  • Reddit
  • Twitter
  • RSS

2 Response to "Metodo de Newton-Raphson para matlab con ejemplo"

  1. alejin de LIGA Says:
    19 de marzo de 2013, 0:19

    good¡¡¡

  2. Tropostibus Says:
    25 de octubre de 2013, 22:10

    A mi no me salio

Publicar un comentario

Buscar en este blog