function b = dualMaxVariance(K, Kj, Y, Yj)
%A function to compute the dual max variance vector using Kj

if (nargin ~= 4)
    fprintf('%s\n', help(sprintf('%s', mfilename)));
    error('Incorrect number of inputs - see above usage instructions.');
end

OPTS.disp = 0;
%OPTS.issym = true;  %Think this one is true too 
OPTS.isreal = true; 
[b, lambda] = eigs(Kj, 1, 'LM', OPTS);

b = b/sqrt(lambda);