add gamma function in GSL