function y=acf(x,MaxLag)

%Returns ACF of the vector x with 'MaxLag' number of lags

L_x=length(x);
r=ones(1,MaxLag+1);
y=zeros(1,MaxLag+1);

lag=1;


while lag<=MaxLag
    
    xlagged=zeros(1,L_x-lag);
    ylagged=zeros(1,L_x-lag);
    for i=1:L_x-lag
        xlagged(i)=x(i+lag);
        ylagged(i)=x(i);
    end
    
    RHO=corrcoef(xlagged,ylagged);
    rho=RHO(1,2);
    
    r(lag+1)=rho;
    
    lag=lag+1;
end

y=r;