Index: /trunk/BNC/src/bncutils.cpp
===================================================================
--- /trunk/BNC/src/bncutils.cpp	(revision 8773)
+++ /trunk/BNC/src/bncutils.cpp	(revision 8774)
@@ -464,5 +464,5 @@
 //
 ////////////////////////////////////////////////////////////////////////////
-int factorial(int n) {
+double factorial(int n) {
   if (n == 0) {
     return 1;
Index: /trunk/BNC/src/pppModel.cpp
===================================================================
--- /trunk/BNC/src/pppModel.cpp	(revision 8773)
+++ /trunk/BNC/src/pppModel.cpp	(revision 8774)
@@ -554,6 +554,6 @@
     for (int m = 0; m <= min(n, M); m++) {
       double pnm = associatedLegendreFunction(n, m, sin(_phiPP));
-      double a = double(factorial(n - m));
-      double b = double(factorial(n + m));
+      double a = factorial(n - m);
+      double b = factorial(n + m);
       if (m == 0) {
         fac = sqrt(2.0 * n + 1);
@@ -570,5 +570,5 @@
 
   if (vtec < 0.0) {
-    return 0.0;
+    vtec = 0.0;
   }