x = [x0];
for i = 1:N
    x0 = f(x0);
    x = [x,x0];
end
    