1
function [N,B,C]=pol2des(Ds)
2
// Given the polynomial matrix Ds= D_0 +D_1 s +D_2 s^2 + ... +D_k s^k,
3
// pol2des returns three matrices N,B,C (with N nilpotent) such that
4
// Ds = C (sN-Eye)^-1 B
7
if type(Ds)==1 then Ds=Ds+0*poly(0,'s')*Ds;end
10
[Sl]=markp2ss(coeff(Ds),dg,nout,nin);
11
N=Sl(2);B=-Sl(3);C=Sl(4)