function bms=gammatonebank(x,lowcf,highcf,numchans,fs)

% Run gammatone filterbank and output selected parameters
% NB outputs weighted by Fay (1988) hearing threshold curve
% modified from code by martin which was modified from code by me!

cfs=MakeErbCFs(lowcf,highcf,numchans);

bms = zeros(numchans,length(x));

for c=1:numchans
  cf=cfs(c);
  bms(c,:) = gammatone(x,fs,cf)/fay(cf);
end;