function J=jury(l)

%Jury's criterion 
%Syntax: J=jury(l)
%
%J - Routh's table
%l - vector of polynomial coefficients of the transfer function denumerator 
% (c) K. Patan, 2012

J=[l;l(end:-1:1)];
m=length(l);
licznik=0;
for j=3:2:2*m-1
    licznik=licznik+1;
    for i=1:m-licznik
      J(j,i)=(J(j-2,1)*J(j-2,i)-J(j-1,1)*J(j-1,i))/J(j-2,1);
    end
    buf=J(j,:); 
    J(j+1,1:m-licznik)=[buf(m-licznik:-1:1)];
end

R=J(1:2:end,1);

if length(find(R(:,1)<=0))>0
     disp('Unstable system')
else
     disp('Stable system')
end 
R