diff --git a/README.txt b/README.txt
new file mode 100644
index 0000000000000000000000000000000000000000..b416537f7ee6b73bdbd07dd02e1f6253b17bda7b
--- /dev/null
+++ b/README.txt
@@ -0,0 +1,28 @@
+This directory provides the code for the experiments
+to the paper
+"Adaptive meshfree approximation for linear
+elliptic partial differential equations with
+PDE-greedy kernel methods" by
+Tizian Wenzel, Daniel Winkle, Gabriele Santin, and Bernard Haasdonk
+
+More specific:
+
+For reproducing the FEM convergence table on the sector example
+of Section 6.1, please
+
+a) Install RBmatlab from the website www.morepas.org/software
+
+The code can be run if RBmatlab from version 16.09 is installed
+
+b) run the script fem_sector_example.m from this directory
+
+Caution:
+
+This operation is very expensive (takes several hours) as hundred of
+thousands of point with global coordinates need to be searched and found
+in FEM meshes. This is required to be consistent in the error computation
+to other approximation techniques using these uniform test grids.
+
+
+
+
diff --git a/fem_sector_example.m b/fem_sector_example.m
index b73fd4536c35e2e351b6ec6988906bb22c49ee3a..58039e437b46eb98170391ead9edacb290b4531f 100644
--- a/fem_sector_example.m
+++ b/fem_sector_example.m
@@ -2,13 +2,16 @@ function fem_sector_example(step)
 %function fem_sector_example(step)
 %
 % Demonstration of FEM error convergence for sector example.
-% very expensive as point with global coordinates need to be
-% searched and found in FEM mesh. But this is required to be consistent
-% with other approximation techniques using these uniform grids.
-% the code can be run if RBmatlab from version 16.09 is installed
-% as obtained from www.morepas.org.
+%
+% This operation is very expensive (takes several hours) as hundred of
+% thousands of point with global coordinates need to be searched and found
+% in FEM meshes. This is required to be consistent in the error computation
+% to other approximation techniques using these uniform test grids.
+%
+% The code can be run if RBmatlab from version 16.09 is installed
+% as obtained from www.morepas.org/software 
 
-% B. Haasdonk 26.7.2023
+% B. Haasdonk 8.5.2024
 
 if nargin < 1
   step = 4;