diff --git a/Packages/Nef_2/changes.txt b/Packages/Nef_2/changes.txt
index be81be081bc..66ec768056e 100644
--- a/Packages/Nef_2/changes.txt
+++ b/Packages/Nef_2/changes.txt
@@ -1,5 +1,8 @@
 Nef_2 Package: Release changes
 ----------------------------------------------------------------------
+1.42 (27. August 2003)
+- Remove non-const operator[] of Polynomial.
+
 1.41 (6. August 2003) [af]
 - Changed to_double(Polynomial)
 
diff --git a/Packages/Nef_2/include/CGAL/Nef_2/Polynomial.h b/Packages/Nef_2/include/CGAL/Nef_2/Polynomial.h
index 1c418c9865e..cd1214bc507 100644
--- a/Packages/Nef_2/include/CGAL/Nef_2/Polynomial.h
+++ b/Packages/Nef_2/include/CGAL/Nef_2/Polynomial.h
@@ -340,10 +340,6 @@ template <class pNT> class Polynomial :
     return ptr()->coeff[i]; }
   //{\Marrop the coefficient $a_i$ of the polynomial.}
 
-  const NT& operator[](unsigned int i) 
-  { CGAL_assertion( i<(ptr()->coeff.size()) );
-    return ptr()->coeff[i]; }
-
   NT operator()(unsigned int i) const
   { 
     if(i<(ptr()->coeff.size()))
@@ -683,10 +679,6 @@ CGAL_TEMPLATE_NULL class Polynomial<int> :
     return ptr()->coeff[i]; }
   /*{\Xarrop the coefficient $a_i$ of the polynomial.}*/
 
-  const int& operator[](unsigned int i) 
-  { CGAL_assertion( i<(ptr()->coeff.size()) );
-    return ptr()->coeff[i]; }
-
   int eval_at(const int& r) const
   /*{\Xop evaluates the polynomial at |r|.}*/
   { CGAL_assertion( degree()>=0 );
@@ -994,10 +986,6 @@ determines the sign for the limit process $x \rightarrow \infty$.
     return ptr()->coeff[i]; }
   /*{\Xarrop the coefficient $a_i$ of the polynomial.}*/
 
-  const double& operator[](unsigned int i) 
-  { CGAL_assertion( i<(ptr()->coeff.size()) );
-    return ptr()->coeff[i]; }
-
   double eval_at(const double& r) const
   /*{\Xop evaluates the polynomial at |r|.}*/
   { CGAL_assertion( degree()>=0 );