function z = H(omega);
mu = 0.12;  
z = ((1-CDF(omega)).*omega + (1-mu).*truncmean(omega));
end