diff --git a/matlab/spins_QSP_csv.m b/matlab/spins_QSP_csv.m
index 63f60d07bd5011503a0587b9d266f0edfb3e7bd9..d544e59fa3bae91a8c6c2480384a0ae9f6cfae0b 100644
--- a/matlab/spins_QSP_csv.m
+++ b/matlab/spins_QSP_csv.m
@@ -1,15 +1,19 @@
-function [T1_max, T1_min, S1_max, S1_min, pdfCount] = spins_QSP_reader(filename)
-  % Reads in the csv file and interprets it in the correct way for the user.
-  %
-  % Input: Filename (string) name of csv file.
-  %
-  % Output: Returns the max/min values for T1 and S1 (as specified in
-  % spins.conf) and the (normalized) pdfCount histogram.
+function [T1_max, T1_min, S1_max, S1_min, pdfCount] = spins_QSP_csv(filename)
+  %% Reads in the csv file and interprets it in the correct way for the user.
+  %%
+  %% Input:
+  %%    Filename (string): name of csv file.
+  %%
+  %% Output:
+  %%    T1_max (float): upper limit of the maximum-valued bin of T1 tracer
+  %%    T1_min (float): lower limit of the minimum-valued bin of T1 tracer
+  %%    S1_max (float): upper limit of the maximum-valued bin of S1 tracer
+  %%    S1_min (float): lower limit of the minimum-valued bin of S1 tracer
   A = importdata(filename);
   T1_max = A(1, 1);
   T1_min = A(1, 2);
   S1_max = A(1, 3);
   S1_min = A(1, 4);
   pdfCount = A(2:end, :);
-  pdfCount /= sum(sum(pdfCount));
+  pdfCount = pdfCount / sum(sum(pdfCount));
 end