--- trunk/src/math/Polynomial.hpp	2009/11/25 20:02:06	1390
+++ trunk/src/math/Polynomial.hpp	2012/08/31 17:29:35	1792
@@ -36,7 +36,8 @@
  * [1]  Meineke, et al., J. Comp. Chem. 26, 252-271 (2005).             
  * [2]  Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006).          
  * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008).          
- * [4]  Vardeman & Gezelter, in progress (2009).                        
+ * [4]  Kuang & Gezelter,  J. Chem. Phys. 133, 164101 (2010).
+ * [5]  Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011).
  */
  
 /**
@@ -273,7 +274,7 @@ namespace OpenMD {
      * @return the first derivative of this polynomial
      */
     PolynomialType & getDerivative() {
-      Polynomial<Real> p();
+      Polynomial<Real> p;
       
       typename Polynomial<Real>::const_iterator i;
       ExponentType exponent;
@@ -336,8 +337,8 @@ namespace OpenMD {
         Real fC0 = getCoefficient(0);
         roots.push_back( -fC0 / fC1);
         return roots;
-      }
         break;      
+      }
       case 2: {
         Real fC2 = getCoefficient(2);
         Real fC1 = getCoefficient(1);
@@ -630,7 +631,7 @@ namespace OpenMD {
    */
   template<typename Real>
   Polynomial<Real> getDerivative(const Polynomial<Real>& p1) {
-    Polynomial<Real> p();
+    Polynomial<Real> p;
     
     typename Polynomial<Real>::const_iterator i;
     int exponent;