% numerical solution of Question 6/Leuven 2005
s=slap(1);
D=(s-2)/((s-1)*(s-3));
newmodel;
f=newsignal(1);
y=newsignal(1);
[k2,ki,g2]=h2hi(f,y,y,D*y+f);
fprintf('\nq6leuven: minimal variance %f at',norm(g2)^2);
tf(k2)