function [c,n, ceq] = covarianceconFun(u, lambda)
c = lambda-0;
n= m-0
%for allowing us to impose the covariance as soft constraint by minimizing
%energy E(u, lambda)
ceq = (integral(norm(x-m,2).*u)./integral(u))-A