function    [Fc] = Fmerit(muu,Lam,X,B,V,H,sigmaa)
%  Call:    [Fc] = Fmerit(muu,Lam,X,B,V,H,sigmaa)
%   calculate complementary (centered) slackness
%  sigmaa for centering for efficency, e.g. .5 rather than 1
n1=size(Lam,1);
Fc = Lam*X - sigmaa*muu * eye(n1);
