function f = H_prime(omega)
mu = 0.12;  
f = 1-CDF(omega)-mu*truncmean_prime(omega);
end