% Miron Andrei-Ionut
% Grupa 341A3

% Laborator 5 IS - MCMMP

clear all
close all
clc

A1 = [1, -0.8];
A2 = [1, -0.4, 0.8];

B1 = [0, 0.7];
B2 = [0, 0.6, 0.16];

N = 100;
nr = 100;

% Testare stabilitate
if all(abs(roots(A1))<1) && all(abs(roots(B1))<1)
    mcmmp_arx1(A1, B1, N, nr);
else
    disp("Polinom instabil ARX1");
end

% Testare stabilitate
if all(abs(roots(A2))<1) && all(abs(roots(B2))<1)
    mcmmp_arx2(A2, B2, N, nr);
else
    disp("Polinom instabil ARX2");
end

% Testare stabilitate
if all(abs(roots(A1))<1) && all(abs(roots(B1))<1)
    mcmmp_oe1(A1, B1, N, nr);
else
    disp("Polinom instabil OE1");
end

% Testare stabilitate
if all(abs(roots(A2))<1) && all(abs(roots(B2))<1)
    mcmmp_oe2(A2, B2, N, nr);
else
    disp("Polinom instabil OE2");
end