get_density_fun function

Get a density function for a gamlss family