function [x, y] = fixedpoint(phi, x0, itpoints, varargin)
%  returns a vector of coordinates, x & y that depicts
%  fixed point iteration
%
%  phi is an input function
%  x0 is an initial guess
%  itpoints is the max number of points
%
	x(1) = x0;
        y(1) = phi(x(1),varargin{:});
        x(2) = y(1);
        y(2) = y(1);
%%	y(2) = feval(phi, x(2), varargin{:});
	for k = 2:2:itpoints
                x(k+1) = x(k);
%%		x(k+1) = feval(phi, x(k), varargin{:});
                y(k+1) = phi(x(k+1),varargin{:});
                x(k+2) = y(k+1);
                y(k+2) = y(k+1);
%%		y(k+2) = feval(phi, x(k+2), varargin{:});
	end
end