From eeb60de46d861f0233177c9f906b28c127d3b378 Mon Sep 17 00:00:00 2001
From: Thomas Applencourt <t.applencourt@gmail.com>
Date: Tue, 5 May 2015 11:56:37 +0200
Subject: [PATCH] solve overflow error in int.f90 for the pseudo

---
 src/Pseudo_integrals/int.f90 | 19 ++++++++++++-------
 1 file changed, 12 insertions(+), 7 deletions(-)

diff --git a/src/Pseudo_integrals/int.f90 b/src/Pseudo_integrals/int.f90
index be806b3f..46e204eb 100644
--- a/src/Pseudo_integrals/int.f90
+++ b/src/Pseudo_integrals/int.f90
@@ -1879,13 +1879,18 @@ end
       RETURN
       END
 
-      double precision function coef_nk(n,k)
-      implicit none
-      integer n,k
-      double precision gam,dblefact,fact
-      gam=dblefact(2*(n+k)+1)
-      coef_nk=1.d0/(2.d0**k*fact(k)*gam)
-      end
+double precision function coef_nk(n,k)
+  implicit none
+  integer n,k
+  double precision gam,dblefact,fact
+
+  if(k.GE.80) then
+    coef_nk = 0.d0
+  else
+    gam=dblefact(2*(n+k)+1)
+    coef_nk=1.d0/(2.d0**k*fact(k)*gam)
+  endif
+end
 
 !!    Calculation of
 !!