diff --git a/Numeric/polynomialBasis.cpp b/Numeric/polynomialBasis.cpp
index 3c3596b186bb35121ff7586a86a05a89296d61e3..2e1c563dc9fb08043aa9c9967f4c4fa408f7ccc3 100644
--- a/Numeric/polynomialBasis.cpp
+++ b/Numeric/polynomialBasis.cpp
@@ -141,7 +141,7 @@ void polynomialBasis::f(const fullMatrix<double> &coord, fullMatrix<double> &sf)
   double p[1256];
   sf.resize (coord.size1(), coefficients.size1());
   for (int iPoint = 0; iPoint < coord.size1(); iPoint++) {
-    evaluateMonomials(coord(iPoint, 0), coord(iPoint, 1),
+    evaluateMonomials(coord(iPoint, 0), coord.size2() > 1 ? coord(iPoint, 1) : 0,
                       coord.size2() > 2 ? coord(iPoint, 2) : 0, p);
     for (int i = 0; i < coefficients.size1(); i++) {
       sf(iPoint,i) = 0.;