function h2syn6245test(nmeas,ncon,nstate)
% function h2syn6245test(nmeas,ncon,nstate)
%
% comparing h2syn6245 to h2syn

if nargin<3, nstate=1; end
if nargin<2, ncon=1; end
if nargin<1, nmeas=1; end
ncost=nstate+ncon;
ndist=nstate+nmeas;
A=randn(nstate);
B=randn(nstate,ndist+ncon);
C=randn(ncost+nmeas,nstate);
D=[zeros(ncost,ndist) randn(ncost,ncon); ...
    randn(nmeas,ndist) zeros(nmeas,ncon)];
p=ss(A,B,C,D,-1);
[K,CL,GAM]=h2syn6245(p,nmeas,ncon);
acl=ssdata(CL);
[~,CL1]=h2syn(p,nmeas,ncon);
fprintf('stability  margin:  %f<1\n',max(abs(eig(acl))))
fprintf('gamma   reduction:  %f>%f=%f\n',norm(CL1),GAM,norm(CL))
fprintf('optimality margin:  %e\n', h2synchk(p,K))