
function E_2 = covarianceobjFun(u,~)
%x=[x,y]
m=integral(x.*u)./integral(u);
lambda=[ 1; 0.5; 0.3; 1.5; 0.6; 1; 1.8; 2; 0.2; 0.16; 1.4; 2];
n= lambda(3);
E_2=sum(int(f(x)).*u,int(g(x)).* abs(Du(x)))+n.*int(A-norm(x-m,2).*u).^2;
