function[A]=matrixex(n, dg, u,v,t,z,w)
% This Matlab routine reconstructs a  generalized semiseparable 
% matrix A from its structured representation of length O(n)
for i=1:n
    A(i,i)=dg(i)+z(i)*conj(w(i));
end
for i=1:n-1
    A(i+1,i)=u(i+1)*conj(v(i));
end
for i=3:n
    tt=1;
    for j=1:i-2
        tt=tt*t(i-j);
        A(i,i-j-1)=u(i)*tt*conj(v(i-j-1));
    end
end
for i=1:n-1
    for j=i+1:n
        A(i,j)=conj(A(j,i)) + z(i)*conj(w(j)) -conj(z(j))*w(i);
    end
end