diff --git a/src/wavefunc.c b/src/wavefunc.c
index 6fa2752..3bcfeac 100755
--- a/src/wavefunc.c
+++ b/src/wavefunc.c
@@ -111,7 +111,7 @@ void gauss(int N,int p,double lb,double ub,double *psi,double *t) {
 		t[i] = lb + delta * i;
 	}
 	
-	den = sqrt(gamma(p+0.5));
+	den = sqrt(cwt_gamma(p+0.5));
 	
 	if ((p+1)%2 == 0) {
 		num = 1.0;