From c92334359b5e3d0ad4a4fae6e40dc9d17beb1fcf Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Wed, 7 Oct 2020 11:23:52 +0200
Subject: [PATCH] Comparison with Stefano's special solution works.

---
 prob_laplace.m | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/prob_laplace.m b/prob_laplace.m
index 11c42b1..df8d1bd 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -26,7 +26,7 @@ T1.plot_sim('plot', 1, 'red');
 T2.plot_sim('plot', 1, 'blue');
 shg
 %% Comparison with Stefano's special case Gamma=Dm/Dp analytical solution
-t = linspace(0, 3, 3);
+t = linspace(0, 10  , 3);
 s = csvread('~/Desktop/langevin_droplet/python_3D_Stefano.csv');
 %%
 T_s = Ternary_model(2, 'Gauss', {-2, b(7/3, 10^-18), 0.5, 0.353553,...
@@ -35,7 +35,7 @@ T_s = Ternary_model(2, 'Gauss', {-2, b(7/3, 10^-18), 0.5, 0.353553,...
 T_s.solve_tern_frap();
 %%
 T_s.plot_sim('plot', 1, 'red');
-plot(s(:, 1), 9*s(:, 2), 'b', 'LineWidth', 3);
+plot(s(:, 1), 9*s(:, 2)./s(:, 1).^2, 'b', 'LineWidth', 3);
 axis([0, 4, 0, 1.5])
 %%
 T_t = T1;
-- 
GitLab