Still not sure whether ternary FRAP solution is correct, diffusion seems too slow towards the end.

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