diff --git a/matlab/peak_detect.m b/matlab/peak_detect.m
index 8cee650..cab685d 100644
--- a/matlab/peak_detect.m
+++ b/matlab/peak_detect.m
@@ -28,6 +28,8 @@ function retval = peak_detect ()
 frames = load ('-ascii', '/home/jens/frames.dat');
 [numFrames, N] = size(frames)
 
+DO_SIM_PEAKS = 0;
+
 Nh = N;
 kn = 0.1
 ks = 0.5;
@@ -49,11 +51,10 @@ peak_count = zeros(1, Nh);
 close all;
 figure
 
-
 for k=1:numFrames
 	% Peak generation
 	
-	if 0
+	if DO_SIM_PEAKS
 		z = fix(Nh*rand() + 1);
 		
 		if (rand() < pa)
@@ -86,6 +87,7 @@ for k=1:numFrames
 	% Update Xm at no-peak bins
 	Xm = (1-gate_upd) .* ((1-am)*Xm + am*X) + gate_upd.*Xm*0.7;
 	
+	% Fill peak gaps with guessed values
 	last = Xm(1);
 	for n=2:Nh,
 		if gate_upd(n) == 0