diff --git a/ternary_frap.m b/ternary_frap.m
index fe8156674ebfc3a0fe4723dbfeccd6084835971e..6fd49a7c5fc4aec263426fbccf65685785a7ceb7 100644
--- a/ternary_frap.m
+++ b/ternary_frap.m
@@ -2,6 +2,10 @@
 % unbleached species. Model is assumed to be equilibrated
 % (bleached+unbleached=const.=pt). Then bleached species initial
 % conditions are introduced. Integration of model via pdepe.
+
+% Still not sure whether this solution is correct. Diffusion seems to slow
+% towards the end!
+
 a = -50;
 b = 0.1;
 c = 1/100000;
@@ -66,4 +70,4 @@ end
 
 function grad = gradient_analytical(x, a, b)
     grad = -(1-tanh(-(x+a)/b).^2)*1/b*0.5;
-end
\ No newline at end of file
+end