I think this will need multivariate polygamma Task list: - [x] Multivariate-polygamma - [ ] Distribution