function [X,degX,Y,degY]=diossl(A,degA,B,degB,C,degC,mes)

%  to solve A*X + B*Y = C over polynomial matrices
%  (state-space realization method)                                          
%                                                         
%  A,degA ... nonsingular polynomial matrix                
%  B,degB ... polynomial matrix                
%  C,degC ... polynomial matrix                
%  X,degX ... polynomial matrix   
%  Y,degY ... polynomial matrix   
%  mes (=0,1,2)...message and interrupt level 
%
%  call: [X,degX,Y,degY]=diossl(A,degA,B,degB,C,degC,mes)

%  function used: transp, diossr

[Ar,degAr]=transp(A,degA);
[Br,degBr]=transp(B,degB);
[Cr,degCr]=transp(C,degC);

[Xr,degXr,Yr,degYr]=diossr(Ar,degAr,Br,degBr,Cr,degCr,0);

[X,degX]=transp(Xr,degXr);
[Y,degY]=transp(Yr,degYr);