diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 043d8855bdb869777f0e97e55d70b54a4121a7aa..7da8b25a73295c6cd4278359d690697cfbafc4f1 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -193,6 +193,10 @@ int main() {
           leafView, p1Basis, neumannNodes, neumannTerm);
       f += neumannTerm;
 
+      VectorType updateTerm(grid.size(grid.maxLevel(), dim));
+      stiffnessMatrix.mv(u1, updateTerm);
+      f += updateTerm;
+
       // {{{ Assemble terms for the nonlinearity
       // TODO: Random value
       std::vector<double> normalStress;