function Y = simul_ar2(phi1,phi2,sigma2,T)
    sigma = sqrt(sigma2);
    Y = zeros(T+2,1);
    for i=3:T+2
        Y(i) = phi1*Y(i-1)+phi2*Y(i-2)+sigma*randn;
    end
    Y = Y(3:end);