function ps62(T)
% function ps62(T)
%
% T  -  horizon for 6.2b and 6.2e (a positive integer, default T=100)

A=[1 0 1;1 1 0;0 1 1];
B=[eye(3) eye(3)];
C=[A;A];
D=[eye(3) eye(3);eye(3) eye(3)];

[K,G]=h2syn(ss(A,B,C,D,-1),3,3);
norm(G)^2
abs(eig(A))