syms s A=[0 1 0; 0 0 1; -2 -1 -1] B=[0;0;1] C=[1 0 0] C*inv(s*eye(3)-A)*B A=[0 1; -1 -1] B=[0;1] C=[1 1] simplify(C*inv(s*eye(2)-A)*B)