function ps41b_6242_2004(n)
% function ps41b_6242_2004(n)
% 
% tests function ps41a_6242_2004

if nargin<1, n=10; end     % default order
s=zeros(n,1);              % begin random generation of s
s(1)=0.1+rand;
for k=2:n,
    s(k)=s(k-1)+0.1+rand;
end
G=ps41a_6242_2004(s);      % using ps41a_6242_2004
close(gcf)
s1=sort(hsvd_6242(G));
[A,B,C,D]=ssdata(G); gsys=pck(A,B,C,D); [sysb,s2]=sysbal(gsys);
plot(1:n,sort(s),'-',1:n,s1,'+',1:n,sort(s2),'o')