c code for gamma from cauchy