diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 54b50bc35517c653d6bbc045fab72c66e6e36432..b375001a406fff75f5e6336d8d3928f55d840b0c 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -230,13 +230,13 @@ int main(int argc, char *argv[]) {
     myAssembler.assembleBodyForce(gravity, density, zenith, gravityFunctional);
 
     // Problem formulation: right-hand side
-    auto const createRHS = [&](double _relativeTime, Vector &_ell) {
+    auto const computeExternalForces = [&](double _relativeTime, Vector &_ell) {
       myAssembler.assembleNeumann(neumannBoundary, _ell, neumannFunction,
                                   _relativeTime);
       _ell += gravityFunctional;
     };
     Vector ell(fineVertexCount);
-    createRHS(0.0, ell);
+    computeExternalForces(0.0, ell);
 
     // {{{ Initial conditions
     ScalarVector alpha_initial(fineVertexCount);
@@ -379,7 +379,7 @@ int main(int argc, char *argv[]) {
       timeSteppingScheme->nextTimeStep();
 
       auto const relativeTime = double(run) / double(timesteps);
-      createRHS(relativeTime, ell);
+      computeExternalForces(relativeTime, ell);
 
       Matrix velocityMatrix;
       Vector velocityRHS(fineVertexCount);