function R=FixedPointSystems(G,x0,tol,kmax)
%
% Function that solves a nonlinear System of Equations
% by generalizing the fixed point method to vectorial functions.
%
% Input  :
%       
%       G   :   handle to the vectorial function 
%               
%       x0  :   column vector - initial approximation
%       toll:   tolerance
%       kmax:   maximum number of iterations
%
% Output :
%
%       R   : Solution of the System
xatt=x0;
% Initializing xprec
xprec=10*x0*tol;
k=1;
while(norm(xprec-xatt)/norm(xprec)>=tol+eps && k<kmax)
    k=k+1;
    xprec=xatt;
    xatt=G(xatt);
end
R=xatt;
end