pbd_brts_density function

Node depth probbaility density for protracted birth-death model of diversification