ymod w/o rec

function Ymod = modelout3(data,b,order,polyorder)

Ymod = zeros(2000,1);


for i = order+1:2000
    u2 = [];
    % linijka nizej u = [y(k-1) y(k-2)  ..  y(k-order) u(k-1) u(k-2) .. u(k-order)]
    u = [fliplr(data(i-order:i-1,2)') fliplr(data(i-order:i-1,1)')];

    for j = 1:polyorder
        u2 = [u2 u.^j]; %dodanie kolejnych czlonow wielomianowych na koniec
    end
    Ymod(i) = u2 * b;
end

end