function x = iter(phi, x0, itmax, m, varargin)
%  returns a vector of iterates, x
%  phi is an input function
%  x0 is an initial guess
%  itmax is the max number of iterations
%  if m = 1 then the Aitken method is performed
	x(1) = x0;
	for k = 1:itmax
		if m ~= 1
			x(k+1) = phi(x(k), varargin{:});
		else
			x(k+1) = (x(k)*phi(phi(x(k), varargin{:})) ...
			- phi(x(k), varargin{:})^2)/(phi(phi(x(k), varargin{:})) ...
			- 2*phi(x(k), varargin{:}) + x(k));
		end
	end
end