Introduce defaults parameters for pdepe.

@Ternary_model/create_mesh.m

@@ -41,4 +41,4 @@ elseif strcmp(T.mode, 'Client') && T.v ~= 0
            linspace(T.x(ind_minus)+0.0001, T.x(ind_plus)-0.0001, 1000),...
            T.x(ind_plus:end)];
 end
-end
+end
\ No newline at end of file

@Ternary_model/solve_tern_frap.m

@@ -7,13 +7,12 @@ fh_bc = @(xl, ul, xr, ur, t) flory_bc(xl, ul, xr, ur, t, T.u0, T.e);
 fh_pde = @(x, t, u, dudx) flory_hugg_pde(x, t, u, dudx, T.a, T.b, T.e,...
                                           T.u0, T.e_g0, T.u_g0, T.v,...
                                           T.mode);
-T.sol = pdepe(T.geometry, fh_pde, fh_ic, fh_bc, T.x, T.t);
+opts = odeset('RelTol',1e-3,'AbsTol',1e-6);
+T.sol = pdepe(T.geometry, fh_pde, fh_ic, fh_bc, T.x, T.t, opts);
 toc
 end
 
 %% Function definitions for pde solver
-
-
 function u = flory_ic(x, a, b, u0, e, ic, x0, ic_c)
 if strcmp(ic, 'FRAP')
     % FRAP initial condition