function [D,GTt,U] = right(mu,n,VT,Z,invZ,R,L,VTiZ,t,RT,a,P,Q,eq,b);
        
% right side of system
% right = 
%[G_{P,Q}( mu*VT*inv(Z)*VT' - VT*inv(Z)*[L+VT'*GT(t)*VT-Z]*R*VT')+zeros(eq,1)];
%                    
% min tr AXBX^t + CX^t
% input: A = A^t ... n*n 
%        B = B^t ... n*n
%        C       ... n*n      
%        matrix R: ((n-1)^2+1)*((n-1)^2+1)
%        matrix Z: ((n-1)^2+1)*((n-1)^2+1)
%        w ... vector: n*n + 1    
% note.. RT = VT*R*VT'  
% 14.06.02, waterloo
% call:  [D,GTt,U] = right(mu,n,VT,Z,invZ,R,L,VTiZ,t,RT,a,P,Q,eq,b);
   
pom = 2: n*n+1;


[GTt] = GT(t,P,Q,eq,n);

I = eye( (n-1)*(n-1) + 1 );
U = (L + VT'*GTt*VT) - Z; 
U1 = VTiZ*(mu*I - U*R)*VT';      
%U1 = VT*invZ*( mu*I - ( L + VT'*BTt*VT - Z)*R )*VT'; 

[GU1] = G(U1,P,Q,eq,n) + b;  % gangster part

D = GU1;
