سلام. برنامه ی زیر صورت و مخرج حلقه بازو میگیره و با یافتن تابع تبدیل حلقه بسته پایداریو بررسی میکنه
n=input('enter num of G(s)=');
d=input('enter den of G(s)=');
na=length(n);
nb=length(d);
D=d+[zeros(1,nb-na),n]; % yaftane moadele moshakhase
r=roots(D);
k=0;
for i=1:length(r)
if...