m=10;
WektorVnLS=[]; %Deklaracja wektora Vn Wektor
WektorAICLS=[]; %Deklaracja wektora AIC
for m=1:10
    FI=[]; %Deklaracja macierzy pomiarów używanych do obliczeń
    FIy=[]; %Deklaracja podmacierzy
    FIu=[]; %Deklaracja podmacierzy
    for a=0:(m-1) %W tej pętli będziemy pobierać odpowiednie dane z pliku i tworzyć wektory FIy i FIu
        FIpom=[-y(m-a:end-a-1),u(m-a:end-a-1)];
        FIy=[FIy,FIpom(:,1)];
        FIu=[FIu,FIpom(:,2)];
    end
    FI=[FIy,FIu]; % Konkatenacja FIy i FIu
    yn=y(m+1:end);
    theta=(inv((FI')*FI))*FI'*yn %Wzór na thete

    En=yn-(FI*theta); %En = Wartość zmierzona - wartość estymacji.
    Vn=En'*En % Wskaźnik kwadratowy
    WektorVnLS=[WektorVnLS,Vn];

    AIC=log(Vn)+(4*m/length(yn)); %Wskaznik AKAIKE
    WektorAICLS=[WektorAICLS,AIC];
    WektorAICLS
end