diff --git a/prob_laplace.m b/prob_laplace.m
index d38125cf2b1ead35659f7e80b2abec97c6f831f4..58e27746b6aad5f1410f5928b2c45fdf5f8c546e 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -186,11 +186,29 @@ for i = 1:100%length(T_mov.t)
     shg; pause();
 end
 
-
+%% Different ratios of Dm/Dp
+chi = 7/3; nu = 10^-6; u0 = 0.5; a = -5;
+b1 = b(chi, nu); P = (u0+e(chi))/(u0-e(chi));
+params = {a, b1, u0, e(chi), -0.15432893, 0, 10, 7, 0, 'Const_mob', 0};
+%% Stefano's way of having mobility with square. works, 0.5% diff
+chi = 7/3; nu = 10^-6; u0 = 0.5; a = -5;
+b1 = b(chi, nu); P = (u0+e(chi))/(u0-e(chi));
+params = {a, b1, u0, e(chi), -0.12314694, 0, 10, 7, 0, 'Const_mob', 0};
+%%
+D_i = 0.05672749; D_o = 0.85091231; chi = 7/3; nu = 10^-6; u0 = 0.5; a = -5;
+b1 = b(chi, nu); P = (u0+e(chi))/(u0-e(chi));
+[~, e_g, u_g] = calc_tanh_params(P, D_i, D_o, a, b1, u0);
+%%
+D_i = 0.03; D_o = 0.83; chi = 7/3; nu = 10^-6; u0 = 0.5; a = -5;
+b1 = b(chi, nu); P = (u0+e(chi))/(u0-e(chi));
+[~, e_g, u_g] = calc_tanh_params(P, D_i, D_o, a, b1, u0);
+%%
+params = {a, b1, u0, e(chi), e_g, u_g, 10, 7, 0, 'Constituent', 0};
 %% %%%%%%%%%%%%%%%%% STEADY STATE JUMP LENGTH DISTRIBUTION %%%%%%%%%%%%%%%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 params = {-5, b(7.7/3, 10^-6), 0.5, e(7.7/3), 0, 1, 10, 7, 0,...
           'Constituent', 0};
+%%
 t = [0, 0.05, 0.1, 1];
 direc = 1;
 x0 = sort(5-direc*(0:0.001:4.01));