function [seq]=fp_iter_func(map, init, tol)
% performs iterations of map
% outputs the sequence of iterates as a column vector

phi=inline(map,'x');

x(1)=init;
x(2)=phi(x(1));
k=2;
while abs( x(k) - x(k-1) ) > tol
    k=k+1;
    if k>100
        break
    end
    x(k)=phi( x(k-1) );
end

seq=x.';
end
