polynomial matrix to descriptor form
[N,B,C]=pol2des(Ds)
:Ds polynomial matrix : :N, B, C three real matrices :
Given the polynomial matrix Ds=D_0 +D_1 s +D_2 s^2 +... +D_k s^k, pol2des returns three matrices N, B, C, with N nilpotent such that:
Ds = C (s*N-eye())^-1 B
s=`poly`_(0,'s');
G=[1,s;1+s^2,3*s^3];[N,B,C]=pol2des(G);
G1=`clean`_(C*`inv`_(s*N-`eye`_())*B),G2=`numer`_(G1)