clc,close all,clear;
%
%
%
syms s K
%
%
%
G = (s+2)/(s*(s-1)*(s+3));
Gz = collect(simplify((K*G/(1+K*G))),s);
%
%
%
[n,d] = numden(Gz);
%
%
%
Gz
collect(d,s)
%
%
%
syms EPS
routh([1,2,K-3,2*K],EPS)