function V=dlyap1(A0,C0);
%  The limiting value is returned in V.

[k,n]  = size(A0);
Cx = inv( eye(k^2) - kron(A0,A0) )*reshape( C0, k^2, 1);
Cx = reshape( Cx, k, k);   

V=Cx;