diff --git a/prob_laplace.m b/prob_laplace.m
index 321695f318e696291379a5995078fafed5f4b845..910b2dac42f49e1a1e3de9459d2f263747b8af09 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -177,8 +177,8 @@ s_dot(i) = sum(diff(T_mov.x).*f.^2./(g0.*(1-pt).*u_interp));
 % plot(x_interp, f)
 % pause()
 end
-csvwrite([s, 'FRAP_Flux.csv'], [x_interp; f])
-csvwrite([s, 'FRAP_entropy.csv'], [T_mov.t; s_dot]);
+% csvwrite([s, 'FRAP_Flux.csv'], [x_interp; f])
+% csvwrite([s, 'FRAP_entropy.csv'], [T_mov.t; s_dot]);
 %% Check integral of solution, should be mass conserving and sum to 1.
 % integrate to 50, to avoid right boundary
 [~, ind] = min(abs(T_prec(4).x-50));
@@ -262,7 +262,7 @@ p = nan(length(ls), n_T);
 for j = 1:n_T
     tic
     parfor i = 1:length(ls)
-        p(i, j) = int_prob(ls(i), F, x0, direction, T_mov, j, 5);
+        p(i, j) = int_prob(ls(i), F, x0, direction, j, T_mov, 5);
     end
     toc
 end