float sgam1(float);
float sgamma(float);