Make mesh more precise close to zero to have good agreement with Stefano.

@Ternary_model/create_mesh.m

@@ -4,7 +4,7 @@ function create_mesh(T)
 
 if strcmp(T.mode, 'Constituent') || strcmp(T.mode, 'Const_mob') || T.v == 0
 % Start with left side of mesh, at r=0/x=0, with a given step size
-    x = linspace(0, -0.75*T.a, 40);
+    x = linspace(0, -0.75*T.a, 200);
     x = [x, x(end)+0.001/T.precision];
     while x(end)<-T.a
         x_temp = x(end)-x(end-1);