m0 = hp.norm()
m1 = hp.norm(loc=4, rate=2)

m_sum = m0 + m1
m_sum.plot_pdf()
assert m_sum.rate() == 3