function H= H_miss2(mu,x0,Z,theta)
[H,DH,D2H]=H_missing(mu,x0,Z,theta);
H=-mean(H);