diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-09 06:02:42 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-09 06:02:42 +0000 |
| commit | 31b53e3d42353fdf1be8167cbb2b4ad9a057f8b6 (patch) | |
| tree | 7ab7986019f4b77029843bef6e9df61f0e9f43e6 | |
| parent | 8aa7f5a5ec78704fa9a86c0cddbda53f36457e03 (diff) | |
| download | multi_mode_eit-31b53e3d42353fdf1be8167cbb2b4ad9a057f8b6.tar.gz multi_mode_eit-31b53e3d42353fdf1be8167cbb2b4ad9a057f8b6.zip | |
proper zeroing of temporal variables
| -rw-r--r-- | liouville.m | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/liouville.m b/liouville.m index 1dc12b0..a67263e 100644 --- a/liouville.m +++ b/liouville.m @@ -66,8 +66,11 @@ for p=1:N % 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; - for w3i=1:length(modulation_freq) + decay_part=0; + Lt=0; + for w3i=1:Nfreq w_iner=modulation_freq(w3i); + decay_part=0; if ((w_iner == w_l)) %such frequency exist in the list of modulation frequencies if ((w_iner == 0)) |