logd_gamma function

dgamma for gamma