1 2 3 4 5 6 7 8
s=%s; num=s+3; den=s^2+2*s+7; H=syslin('c',num,den) SS=tf2ss(H) [Ac,Bc,U,ind]=canon(SS(2),SS(3))