

function [ tabla ] = newton_raphson(a,n,delta,tole)


format long

%criterios de inicializacion

i=0;
funcionvalue= delta+1;
error = tole+1;


% Metodo de Newton-Raphson

	while i<n & funcionvalue>delta & error>tole & dfuncion(a)~=0

		i=i+1; % Numero de iteraciones a realizar
		tabla(i,1)=i;
		
		x=a-(funcion(a)/dfuncion(a)); %Ecuacion que define el metodo de Newton-Raphson
		tabla(i,2)=x ;
		
		y=funcion(x); %Evaluar x en la funcion para hayar el intervalo
		tabla(i,3)=y;
		
		funcionvalue=abs(y);  %para poder calcular el punto en el eje y donde puede estar la solucion
		error=abs(x-a); %el calculo del error en el punto
		%error=abs(x-a)/abs(x); %error relativo
		%error=(abs(x-a)/abs(x))*100; %error porcentual
		tabla(i,4)=error;
	
		a=x;
		
	
	endwhile
	


endfunction
