gammp
in nr2.h
extern float gammp(float a, float x);