function X=iterations(N,x0,y0)
    if N==0 then X=[x0;y0];
    else 
        X=iterations(N-1,x0,y0);
        X=[X [0 1; 1 1]*X(:,N)];
    end
endfunction