function k = H_prime_prime(omega)
mu = 0.12;  

k = -CDF_prime(omega) - mu*truncmean_prime_prime(omega);
end