diff --git a/Numeric/gsl_newt.cpp b/Numeric/gsl_newt.cpp
index 934e7a88f3f7fb54610bbc839c4f6f05b45d6232..bba2b52eca1542b4313a6a25f25931d17bc5bc00 100644
--- a/Numeric/gsl_newt.cpp
+++ b/Numeric/gsl_newt.cpp
@@ -1,4 +1,4 @@
-// $Id: gsl_newt.cpp,v 1.4 2003-02-20 16:34:30 geuzaine Exp $
+// $Id: gsl_newt.cpp,v 1.5 2003-02-20 16:49:21 geuzaine Exp $
 //
 // Copyright (C) 1997 - 2003 C. Geuzaine, J.-F. Remacle
 //
@@ -77,7 +77,7 @@ void newt(double x[], int n, int *check, void (*func)(int, double [],double []))
   gsl_multiroot_function f = {&gslfunc, n, &p};
   gsl_vector *xx = gsl_vector_alloc (n);
 
-  if(n > MAX_DIM_NEWT) Msg(FATAL, "Maximum Newton dimension exceeded\n");
+  if(n > MAX_DIM_NEWT-1) Msg(FATAL, "Maximum Newton dimension exceeded\n");
   gsl_dim = n;
 
   nrfunc = func;