Newton-Raphson

%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


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

%Funcion derivada:
function d=derivada(x)
d=3*(x^2)-(60*x);