Need to implement gamma and polygamma functions
Need to implement gamma and polygamma functions