%EqSolveNewton: Solve nonlinear system of equations f(x)=0 using the 
% Newton method. 
%
% Usage:
%   >> [ Xk , k ] = EqSolveNewton( f , Jf , X0  )
%
% where f(x) is the function, Jf is the Jacobian function, and 
% X0 is the starting approximation.
%
%
function [ Xk , k ] = EqSolveNewton( f , Jf , Xk  )

   tol=10^-10;
      
   
   while norm( Sk )>tol,
       %
       % Evaluate function f(Xk), Jacobian Jf(Xk), and compute
       % the next Newton step Sk. Use feval() to evaluate functions.
       %
       
   end
