diff --git a/Jet_fitting_3/include/CGAL/Monge_via_jet_fitting.h b/Jet_fitting_3/include/CGAL/Monge_via_jet_fitting.h
index 9ae27202cb8..425ad70176e 100644
--- a/Jet_fitting_3/include/CGAL/Monge_via_jet_fitting.h
+++ b/Jet_fitting_3/include/CGAL/Monge_via_jet_fitting.h
@@ -422,9 +422,16 @@ fill_matrix(InputIterator begin, InputIterator end,
     y = itb->y();
     //  Z[line_count] = itb->z();
     Z.set(line_count,itb->z());
-    for (std::size_t k=0; k <= d; k++) for (std::size_t i=0; i<=k; i++)
-      M.set(line_count, k*(k+1)/2+i, std::pow(x,k-i)*std::pow(y,i)
-		/(fact(i)*fact(k-i)*std::pow(this->preconditionning,k)));
+    for (std::size_t k=0; k <= d; k++) {
+      for (std::size_t i=0; i<=k; i++) {
+        M.set(line_count, k*(k+1)/2+i,
+              std::pow(x,static_cast<double>(k-i))
+              * std::pow(y,static_cast<double>(i))
+              /( fact(i) *
+                 fact(k-i)
+                 *std::pow(this->preconditionning,static_cast<double>(k))));
+      }
+    }
     line_count++;
   }
 }