function f = obj_sylv_p2_3(x_param,GAMMA1,GAMMA2)
%
% objectieffunctie paper 2 (EUS '00) met een stap minder in de berekening dan
% obj_sylv_p2_2 

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))];
%x = x_init.' /norm(x_init); 

 [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]);