function Detect(input, cutoff)

if(ischar(input))
    load(input)
else
    A = input;
end

if(nargin < 2)
    cutoff = floor(size(A,1)/2);
end

for rowrange = 1:2
    cutoff1 = cutoff;
    cutoff2 = min(cutoff1+40000, size(A,1));
    if(rowrange == 1)
	    rows = [1:cutoff1]; 
    elseif(rowrange == 2)
	    rows = [cutoff1+1:cutoff2];
    end
	columns = [4:10];
	
	fadj = 1; %0.995; 
    fsample = 1000;
	in = 1:floor(size(rows,2)/20); 
	filter1000 = sin(1000*fadj*(in)*2*3.14159/fsample);
	filter316 = sin(316.2*fadj*(in)*2*3.14159/fsample);
	filter100 = sin(100*fadj*(in)*2*3.14159/fsample);

    for i = columns
        b(:,i) = A(rows,i);
        b1000(:,i) = conv(b(:,i), filter1000); 
        b316(:,i) = conv(b(:,i), filter316); 
        b100(:,i) = conv(b(:,i), filter100); 
	end;
	
	figure; 
	col = length(columns);
	for(i = columns)
        subplot(3,col,i-columns(1)+1); plot(b1000(:,i), '.');
        subplot(3,col,i-columns(1)+1+col); plot(b316(:,i), '.');
        subplot(3,col,i-columns(1)+1+2*col); plot(b100(:,i), '.');
	end;
	clear b;
	clear b1000;
	clear b316;
	clear b100; 
end

return