Shorten probability variable name.

prob_laplace.m

@@ -405,16 +405,16 @@ bp = 5;
 direc = -2; % 1: jump from left to right, -1: jump from right to left.
 
 %% 
-[ls, p_in_in] = frap_distro(t, params, 2, bp);
-[~, p_out_out] = frap_distro(t, params, -2, bp);
-[~, p_in_out] = frap_distro(t, params, 1, bp);
-[~, p_out_in] = frap_distro(t, params, -1, bp);
+[ls, p_i_i] = frap_distro(t, params, 2, bp);
+[~, p_o_o] = frap_distro(t, params, -2, bp);
+[~, p_i_o] = frap_distro(t, params, 1, bp);
+[~, p_o_i] = frap_distro(t, params, -1, bp);
 %%
 figure; hold on;
-plot(ls, sum(p_in_in, 2));
-plot(ls, sum(p_out_out, 2));
-plot(ls, sum(p_in_out, 2));
-plot(ls, sum(p_out_in, 2));
+plot(ls, sum(p_i_i, 2));
+plot(ls, sum(p_o_o, 2));
+plot(ls, sum(p_i_o, 2));
+plot(ls, sum(p_o_i, 2));
 %% Save time points for comparison with python
 csvwrite('lb_in_out.csv', ls)
 csvwrite('prob_in_out.csv', p);