Removing c from pde. Now: c=1.

ternary_frap.m

@@ -5,13 +5,13 @@
 pa = '/Users/hubatsch/ownCloud/Dropbox_Lars_Christoph/DropletFRAP/FRAP_paper/';
 a = -1;
 b = 0.025;
-c = 1/100;
 u0 = 0.05;
 e = 0.4;
 e_g0 = 0.16;
 o_g0 = 0.2;
 %%
 g0 = gamma0(x, a, b, e_g0, o_g0);
+pt = phi_tot(x, a, b, e, u0);
 %% Make useful mesh (by inverting the tanh profile and using this as spacing)
 x = linspace(-a-1, -a+1, 3000);
 g = gamma0(x, a, 1*b, e_g0, o_g0);
@@ -63,7 +63,7 @@ plot(x, gamma0(x, a, b, e));
 figure; plot(gamma0(x, a, b, e), spacing(gamma0(x, a, b, e), a, b, e));
 %% Function definitions for pde solver
 function [c, f ,s] = flory_hugg_pde(x, t, u, dudx, a, b, e, c_p, u0,...
-                                    e_g0, o_g0)  
+                                    e_g0, o_g0)
 % Solve with full ternary model. Analytical derivatives.
 % pt ... phi_tot
 % gra_a ... analytical gradient of phi_tot
@@ -71,7 +71,7 @@ function [c, f ,s] = flory_hugg_pde(x, t, u, dudx, a, b, e, c_p, u0,...
 pt = phi_tot(x, a, b, e, u0);
 gra_a = gradient_analytical(x, a, b, e);
 g0 = gamma0(x, a, b, e_g0, o_g0);
-c = c_p;
+c = 1;
 f = g0*(1-pt)/pt*(pt*dudx-u*gra_a);
 s = 0;
 end