function rga = vrga(G)
   [mt,ny,nu,nx] = minfo(G); rga = []; 

   if strcmp( mt, 'syst' ) == 1 
       disp( 'Not implemented for system matrices.' );
       return;      
   end
   
   rga = veval('.*',G,vpinv(vtp(G)));
end
