--- a/comp/comp_warpedfoff.m
+++ b/comp/comp_warpedfoff.m
@@ -1,4 +1,4 @@
 function foff=comp_warpedfoff(fc,bw,fs,L,scaletofreq)
 %COMP_WARPEDFOFF  foff for warped filters
 
-foff=floor(scaletofreq(fc-.5*bw)/fs*L);
+foff=floor(scaletofreq(fc-.5*bw)/fs*L)+1;