function f = obj_sylv_p2_2(x_param,GAMMA1,GAMMA2)
%
% objectieffunctie paper 2 (EUS '00) 

G11 = GAMMA1(:,1:3); G12  = GAMMA1(:,4:6);
G21 = GAMMA2(:,1:3); G22  = GAMMA2(:,4:6);

x = [cos(x_param(1)); sin(x_param(1))*exp(i*x_param(2))]; 

 [U,S,V] = svd( [x(1)*G11+x(2)*G12; conj(x(1))*G21+conj(x(2))*G22] );
 y = V(:,3);
 [U,S,V] = svd( [G11*y G12*y; conj(G21*y) conj(G22*y)] );
 x = V(:,2);

q1est = kron(x,y);			% norm(q1est) = 1
r1 = GAMMA1*q1est;
q1estbis = kron(conj(x),y);		% norm(q1estbis) = 1
r2 = GAMMA2*q1estbis;

f = norm([r1; r2]);