function z = stab_theta(th,ar,ai)

for k=1:length(ai)
  a(k,:)=ar+i*ai(k);
end

z = abs((1+(1-th)*a)./(1-th*a));