diff options
author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-09 06:01:39 +0000 |
---|---|---|
committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-09 06:01:39 +0000 |
commit | 8aa7f5a5ec78704fa9a86c0cddbda53f36457e03 (patch) | |
tree | a698b5858d629ca1c8d929e21093a064f32256c7 /liouville.m | |
parent | 35af744d95df0d687bccdaaa4dabda6fe830bed7 (diff) | |
download | multi_mode_eit-8aa7f5a5ec78704fa9a86c0cddbda53f36457e03.tar.gz multi_mode_eit-8aa7f5a5ec78704fa9a86c0cddbda53f36457e03.zip |
decay part moved to the end result
Diffstat (limited to 'liouville.m')
-rw-r--r-- | liouville.m | 4 |
1 files changed, 1 insertions, 3 deletions
diff --git a/liouville.m b/liouville.m index 67dde6b..1dc12b0 100644 --- a/liouville.m +++ b/liouville.m @@ -75,8 +75,6 @@ for p=1:N decay_part=\ ( decay_total(g_decay,k)/2 + decay_total(g_decay,j)/2 + g_dephasing(j,k) )* kron_delta(j,m)*kron_delta(k,n) \ - kron_delta(m,n)*kron_delta(j,k)*g_decay(m,j) ; - decay_part=decay_part*hbar/im_one; - L0+=decay_part; Lt=L0; else Li= ( dipole_elements(j,m)*kron_delta(k,n)-dipole_elements(n,k)*kron_delta(j,m) )*E_field(w3i); @@ -84,7 +82,7 @@ for p=1:N endif %Lt=-im_one/hbar*Lt*kron_delta(w_jk-w_iner,w_mn); % above if should be done only if kron_delta is not zero % no need for above kron_delta since the same conditon checked in the outer if statement - Lt=-im_one/hbar*Lt; + Lt=-im_one/hbar*Lt - decay_part; endif endfor if ((p == s)) |