From 31b53e3d42353fdf1be8167cbb2b4ad9a057f8b6 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Wed, 9 Dec 2009 06:02:42 +0000 Subject: proper zeroing of temporal variables --- liouville.m | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) 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)) -- cgit v1.2.3