diff options
author | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-03-23 18:30:34 +0000 |
---|---|---|
committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-03-23 18:30:34 +0000 |
commit | 964a4b2c8eac531f12511845f2945ca36caba7fd (patch) | |
tree | acd6205dcfa4b6a8c4c54175f68711409f752508 /useful_functions.m | |
parent | f2d5325ab13bc5adfba3c1b0885267d886c34e4b (diff) | |
download | multi_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.m | 9 |
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 ... |