Hw1 = exp((fun(freq) + log(Hw2)) / 2)