function b = dualMaxTargetAlignment(K, Kj, Y, Yj)
%A function to compute the dual maximum target alignment vector between matrices X and Y
%
%inputs 
%K kernel matrix
%Y labels vector 
%
%outputs 
%b the vector which maximises the emperical target alignment between K and Y 

[B, L] = eigs(pinv(Kx)*Ky*Kx)

b = b(:, 1); 
b = b/sqrt(b'*Kj'*Kj*b);