Change the input "delta" parameter, so that it is a float, without a ten to the power...

diff --git a/SimuPBI_circle_func.m b/SimuPBI_circle_func.m
index 9b7f016d6317a28de1f7fc34fc25a1c64c4c6cca..86fe61b6df66a3b0f9e1fbf30a90681c767b8031 100644 (file)
--- a/SimuPBI_circle_func.m
+++ b/SimuPBI_circle_func.m
@@ -20,7 +20,7 @@
 ## Created: 2016-02-11
 
 ## Example of use :
-## SimuPBI_circle_func('Analytical', 2, 'test', [0 0.1 0.5], 19, 2.7, 300, 180, 1, 1, 0.89, 8.265*10^-7, 100, 500, [200 250], '', 'gaussian', 35)
+## SimuPBI_circle_func('Analytical', 2, 'test', [0 0.1 0.5], 19, 2.7, 300, 180, 1, 1, 0.89, 8.265, 100, 500, [200 250], '', 'gaussian', 35)
 
 ## Please note that dir_out, noise_type, noise_amount are optional arguments.
 
@@ -76,7 +76,7 @@ addpath('utilities_LW')
 
 ## material of the circle (here, PET at 19 keV)
 %mu_mat=0.89889 % (in cm-1); 
-%delta_mat=8.265*10^-7;
+%delta_mat=8.265; (in 10^-7)
 
 ## size and coordinates
 %R_circle=50; % radius of the circle
@@ -160,6 +160,7 @@ end
 
 % absorption and phase parameter 
 mu_mat=mu_mat*100; % in m-1 , to be consistant with lambda!!
+delta_mat=delta_mat*10^-7;
 beta_mat=lambda*mu_mat/(4*pi);
 
 %% Cartes 2D d'attenuation et d/b