--- trunk/OOPSE/libmdtools/BondExtensions.cpp	2003/06/24 19:57:54	564
+++ trunk/OOPSE/libmdtools/BondExtensions.cpp	2003/07/11 22:34:48	594
@@ -29,8 +29,8 @@ double HarmonicBond::bond_force( double r_ab ){
   dr = r_ab - d0;
   dr2 = dr * dr;
 
-  c_potential_E = k0 * dr2;
-  force = -2.0 * k0 * dr;
+  c_potential_E = 0.5* k0 * dr2;
+  force = - k0 * dr;
   return force;
  
 }