diff --git a/test/unitvectortest.cc b/test/unitvectortest.cc
index 26cfe8ac37c32ec3140dfdc041d2f4dc0ca1b96a..2fbe8826a2e04cb1bd65388959a5af69f1268dbc 100644
--- a/test/unitvectortest.cc
+++ b/test/unitvectortest.cc
@@ -111,8 +111,8 @@ void testDerivativesOfSquaredDistance(const UnitVector<dim>& a, const UnitVector
 
             FieldVector<double,dim> bPlus  = b.globalCoordinates();
             FieldVector<double,dim> bMinus = b.globalCoordinates();
-            bPlus[i]  += eps;
-            bMinus[i] -= eps;
+            bPlus[j]  += eps;
+            bMinus[j] -= eps;
                 
             d1d2_fd[i][j] = (energy(aPlus,bPlus) + energy(aMinus,bMinus)
                             - energy(aPlus,bMinus) - energy(aMinus,bPlus)) / (4*eps*eps);