summaryrefslogtreecommitdiff
path: root/useful_functions.m
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2010-03-23 18:30:34 +0000
committerEugeniy Mikhailov <evgmik@gmail.com>2010-03-23 18:30:34 +0000
commit964a4b2c8eac531f12511845f2945ca36caba7fd (patch)
treeacd6205dcfa4b6a8c4c54175f68711409f752508 /useful_functions.m
parentf2d5325ab13bc5adfba3c1b0885267d886c34e4b (diff)
downloadmulti_mode_eit-964a4b2c8eac531f12511845f2945ca36caba7fd.tar.gz
multi_mode_eit-964a4b2c8eac531f12511845f2945ca36caba7fd.zip
added number of participating atoms into the calculation
Diffstat (limited to 'useful_functions.m')
-rw-r--r--useful_functions.m9
1 files changed, 6 insertions, 3 deletions
diff --git a/useful_functions.m b/useful_functions.m
index 9127cab..842f0fa 100644
--- a/useful_functions.m
+++ b/useful_functions.m
@@ -303,10 +303,13 @@ function relative_transmission=total_relative_transmission(atom_field_problem)
modulation_freq = atom_field_problem.modulation_freq ;
pos_freq_indexes=(modulation_freq>0);
+ % WARNING INTRODUCED HERE BUT MUST BE DEFINE IN ATOM PROPERTIES FILE
+ %n_atoms - proportional to a number of interacting atoms
+ n_atoms=500;
transmited_intensities_vector = ...
- abs(E_field.linear + (1i)*dEdz_linear).^2 ...
- +abs(E_field.right + (1i)*dEdz_right).^2 ...
- +abs(E_field.left + (1i)*dEdz_left).^2;
+ abs(E_field.linear + n_atoms*(1i)*dEdz_linear).^2 ...
+ +abs(E_field.right + n_atoms*(1i)*dEdz_right).^2 ...
+ +abs(E_field.left + n_atoms*(1i)*dEdz_left).^2;
transmited_intensities_vector = transmited_intensities_vector(pos_freq_indexes);
input_intensities_vector = ...
abs(E_field.linear).^2 ...