Orthogonal polynomials for tailoring density functions to excess kurtosis, asymmetry and dependence