//
// Resolucion de sistemas no-lineales
//
// Metodo de la regla falsa
//

clear

// Grafico de la funcion

// Cargando funciones grabadas en archivos diferentes
getf('f.sci')
getf('fp.sci')

//deff('y=f(x)','y=log(x.^2+1)-exp(0.4*x).*cos(%pi*x)')


xp = -3:0.1:2;
yp = f(xp);


figure(1)
plot2d(xp,yp)
plot2d(xp,zeros(xp))


q = 5;      // Numero de digitos de precision para la convergencia
n_iter = 100;

x(1) = input('Valor inicial de x ? =');

for k =2:n_iter
	   x(k) = x(k-1)-f(x(k-1))/fp(x(k-1));
		 plot2d(xp,yp)
		 plot2d(xp,zeros(xp))
   		plot2d(x(k-1), f(x(k-1)),[-4]);
   		plot2d([x(k) x(k-1)]', [ 0 f(x(k-1))]',[2]);
      pause

   
   // Verificacion del criterio de parada
   
	  if abs( x(k)- x(k-1) ) <= 2*abs( x(k) )*10^(-q) 
   	   break
	  end
   
  
end

x($)
y=f(x($))