diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-10 17:20:14 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-10 17:20:14 +0000 |
| commit | dbbb852d00038ac7c7873783c4442724dee1031b (patch) | |
| tree | 496ea15efe7d2564fa0b3526f541030abade9422 | |
| parent | 92651328a20576a67a9ff5b92c71573cb7eb7024 (diff) | |
| download | multi_mode_eit-dbbb852d00038ac7c7873783c4442724dee1031b.tar.gz multi_mode_eit-dbbb852d00038ac7c7873783c4442724dee1031b.zip | |
some small typos fixes
| -rw-r--r-- | useful_functions.m | 13 |
1 files changed, 6 insertions, 7 deletions
diff --git a/useful_functions.m b/useful_functions.m index 2ae4ced..1765dd9 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -59,9 +59,8 @@ function L=Liouville_operator_matrix( modulation_freq, rhoLiouville_w, rhoLiouville_r, rhoLiouville_c ) %------------------------- - useful_constants; L=zeros(N); % NxN matrix - + useful_constants; Nfreq=length(modulation_freq); for p=1:N @@ -70,16 +69,15 @@ function L=Liouville_operator_matrix( k=rhoLiouville_c(p); m=rhoLiouville_r(s); n=rhoLiouville_c(s); - % we garanted to know frequency of final and initial rhoLiouville + % we guaranteed to know frequency of final and initial rhoLiouville w1i=rhoLiouville_w(p); w2i=rhoLiouville_w(s); w_jk=modulation_freq(w1i); w_mn=modulation_freq(w2i); % thus we know L matrix element frequency which we need to match w_l=w_jk-w_mn; - % lets search this wrequency in the list of available frequencyes - % but since we not garanteed to find it lets assign temporary 0 to Liouville matrix element - L(p,s)=0; + % lets search this frequency in the list of available frequencies + % but since we not guaranteed to find it lets assign temporary 0 to Liouville matrix element decay_part=0; Lt=0; for w3i=1:Nfreq @@ -108,7 +106,8 @@ function L=Liouville_operator_matrix( if ((p == s)) Lt+=-im_one*w_jk; endif - L(p,s)=Lt; + Lps=Lt; + L(p,s)=Lps; endfor endfor endfunction |