A = [1 -0.4 -0.32];
B = [0 0.5 0.03];

N = 100;
nr = 100;

if all(abs(roots(A))<1) && all(abs(roots(B))<1)
    mcmmp_arx2(A, B, N, nr);
else
    disp("Polinom instabil");
end

if all(abs(roots(A))<1) && all(abs(roots(B))<1)
    mcmmp_oe2(A, B, N, nr);
else
    disp("Polinom instabil");
end

if all(abs(roots(A)) < 1) && all(abs(roots(B)) < 1)
    mcmmp_arx_Pb3(A, B, N, nr);
else
    disp("Polinom instabil");
end

if all(abs(roots(A)) < 1) && all(abs(roots(B)) < 1)
    mcmmp_oe_Pb3(A, B, N, nr);
else
    disp("Polinom instabil");
end