diff --git a/src/evaluateFarfieldSecondOrder.m b/src/evaluateFarfieldSecondOrder.m
index ab1d20a9bd1ac3dd8a38a79b8cc0e12ea20e97cc..b9eb4df350d182277f87e2bbf8efeac2e01e7868 100644
--- a/src/evaluateFarfieldSecondOrder.m
+++ b/src/evaluateFarfieldSecondOrder.m
@@ -1,8 +1,9 @@
 function [F] = evaluateFarfieldSecondOrder(k, sampling, objects, par, R, q)
 
-% Calculates a discretized version of the Born farfield operator of second order 
-% by using the composite trapezoidal rule for the outer integral and calculations
-% as suggested by Vainikko in order to approximate the individual Born farfields.
+%% EVALUATEFARFIELDSECONDORDER: Calculates a discretized version of the Born far field
+% operator of second order by using the composite trapezoidal rule for the outer
+% integral and calculations as suggested by Vainikko in order to approximate the
+% individual Born farfields.
 %
 % INPUT: kappa      Wave number, >0.
 %        sampling   Structure containing information about the discretization.
@@ -17,6 +18,8 @@ function [F] = evaluateFarfieldSecondOrder(k, sampling, objects, par, R, q)
 %   
 % OUTPUT: F     Born farfield operator of second order, nxhat*nd-array.
 %
+% SYNTAX: evaluateFarfieldSecondOrder(k, sampling, objects, par, R, q)
+%
 % ************************************************************************************
 
 kR = k * R;