Polynomial s; Rational G; Matrix A, B, C, D; s = Polynomial("s"); G = 1/(s^2 + 2 * s + 1); {A,B,C,D} = tfm2ss([G]);