diff --git a/Makefile b/Makefile
index f07271e..f34b0b9 100644
--- a/Makefile
+++ b/Makefile
@@ -52,5 +52,5 @@ ${OBJ_DIR}/ex1.o: ${SRC_DIR}/ex1.c
 	-${CC} ${SLEPC_INCLUDE} ${PETSC_CC_INCLUDES} -c -o $@ $< ${SLEPC_EPS_LIB}
 
 ${BIN_DIR}/ex1: ${OBJ_DIR}/get_ntot.o ${OBJ_DIR}/read2.o ${OBJ_DIR}/get_s2_mov.o ${OBJ_DIR}/get_s2_cyclic.o ${OBJ_DIR}/get_s2.o ${OBJ_DIR}/get_dmat.o ${OBJ_DIR}/get_val_iaa2.o ${LIB_DIR}/irpf90.a ${OBJ_DIR}/ex1.o ${SRC_DIR}/read2.h ${SRC_DIR}/get_ntot.h ${SRC_DIR}/stimsyr.h chkopts
-	    -${CLINKER} ${SLEPC_INCLUDE} ${PETSC_CC_INCLUDES} -o ${BIN_DIR}/ex1 ${OBJ_DIR}/ex1.o ${OBJ_DIR}/read2.o ${OBJ_DIR}/get_ntot.o ${OBJ_DIR}/get_s2.o ${OBJ_DIR}/get_s2_mov.o ${OBJ_DIR}/get_s2_cyclic.o ${OBJ_DIR}/get_dmat.o ${OBJ_DIR}/get_val_iaa2.o ${LIB_DIR}/irpf90.a ${SLEPC_EPS_LIB}# -lifcore -lirc -lcomposerxe_gen_helpers_core_2.3
+	    -${CLINKER} ${SLEPC_INCLUDE} ${PETSC_CC_INCLUDES} -o ${BIN_DIR}/ex1 ${OBJ_DIR}/ex1.o ${OBJ_DIR}/read2.o ${OBJ_DIR}/get_ntot.o ${OBJ_DIR}/get_s2.o ${OBJ_DIR}/get_s2_mov.o ${OBJ_DIR}/get_s2_cyclic.o ${OBJ_DIR}/get_dmat.o ${OBJ_DIR}/get_val_iaa2.o ${LIB_DIR}/irpf90.a ${SLEPC_EPS_LIB}
 #    ${RM} ex1.o read2.o
diff --git a/examples/small.inp b/examples/small.inp
new file mode 100644
index 0000000..4c0c4ae
--- /dev/null
+++ b/examples/small.inp
@@ -0,0 +1,14 @@
+8
+140
+1
+1
+0
+true
+1,2,3,1,2,3,4,5,6,7
+2,3,4,8,7,6,5,6,7,8
+1,1,1,2,2,2,2,3,3,3
+.1430,-0.20,0.0000
+.1430,-0.20,0.0000
+-1.00,0.0,0.00
+1
+1
diff --git a/examples/three.inp b/examples/three.inp
new file mode 100644
index 0000000..8706cf2
--- /dev/null
+++ b/examples/three.inp
@@ -0,0 +1,14 @@
+3
+1
+1
+1
+1
+true
+1,2,1,2,3,4,5
+2,3,6,5,4,5,6
+1,1,2,2,2,3,3
+.1430,-0.20,0.0000
+.1430,-0.20,0.0000
+-1.00,0.0,0.00
+1
+1
diff --git a/src/ex1.c b/src/ex1.c
index e062938..c2a2659 100644
--- a/src/ex1.c
+++ b/src/ex1.c
@@ -123,9 +123,9 @@ int main(int argc,char **argv)
             &getdata.FAM1,
             tcol,
             val);
-      if(i%getdata.npar == 0 && mpiid==0){
-        ierr = PetscPrintf(PETSC_COMM_WORLD," i: %d \n",i);CHKERRQ(ierr);
-      }
+//    if(i%getdata.npar == 0 && mpiid==0){
+//      ierr = PetscPrintf(PETSC_COMM_WORLD," i: %d \n",i);CHKERRQ(ierr);
+//    }
       for(ll=0;ll<getdata.nnz;ll++){
 
         for(kk=0;kk<tcountcol[ll]+1;kk++){
@@ -164,7 +164,6 @@ int main(int argc,char **argv)
   tol = 1.e-9;
   maxit = 10000000;
   ierr = EPSSetTolerances(eps,tol,maxit);CHKERRQ(ierr);
-  nev  = 4;
   ncv  = 10;
   mpd  = 10;
   nev  = getdata.nroots;
diff --git a/src/nt1.irp.f b/src/nt1.irp.f
index 55e7971..8bd0124 100644
--- a/src/nt1.irp.f
+++ b/src/nt1.irp.f
@@ -15,6 +15,6 @@ BEGIN_PROVIDER [integer(C_SIZE_T), nt1]
             natom2 = fix_trou2 - fix_trou1
         endif
     endif
-    nt1=   nint(gamma(real(natom2+1,16))/(gamma(real(natom2-ntrou+1,16))*gamma(real(ntrou+1,16))),selected_int_kind(16))
+    nt1=   nint(gamma(1.0*(natom2+1))/(gamma(1.0*(natom2-ntrou+1))*gamma(1.0*(ntrou+1))),selected_int_kind(16))
     write(6,*)'nt1',nt1
 END_PROVIDER
diff --git a/src/nt2.irp.f b/src/nt2.irp.f
index 7becdac..0ff4c5d 100644
--- a/src/nt2.irp.f
+++ b/src/nt2.irp.f
@@ -6,6 +6,6 @@ BEGIN_PROVIDER [integer(C_SIZE_T), nt2]
 
 
 !   call combin(idet2(1,nt2+1),natrest,ial0,nt2,32,jrangmax)
-    nt2=   nint(gamma(real(natom-ntrou+1,16))/((gamma(real(nalpha+1,16))*gamma(real(nbeta+1,16)))),selected_int_kind(16))
+    nt2=   nint(gamma(1.0*(natom-ntrou+1))/((gamma(1.0*(nalpha+1))*gamma(1.0*(nbeta+1)))),selected_int_kind(16))
     print *,"nt2=",nt2
 END_PROVIDER