function[A,B,C,D]=realization();
[a11 b11 c11 d11]=tf2ss([0 0 1],[1 0 0]);
[a12 b12 c12 d12]=tf2ss([1],[1 -1]);
[a21 b21 c21 d21]=tf2ss([0 1],[1 0]);
A=[a11 zeros(2,3); zeros(3,2) diag([a12 a21 0])];
B=[b11 [0;0] ; 0 b12 ; b21 0 ; 0 0];
C=[c11 c12 0 0 ; 0 0 0  c21 0];
D=zeros(2,2);