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