function ps64(r)
if nargin<1, r=-1; end
ga=0.5*(r+sqrt(r^2-4*r));
a=sqrt(-r);
[~,g]=lqr(ss(1,1,0,0,-1),-r,1-r,-r);
fprintf('r=%f:   %f/%f\n',r,ga,g)