function yyy = bm_haircell(y, nchan, flo)
% bm_haircell(stimulus)
% Basilar-membrane / Haircell model based on Malcolm Slaney's implementation
% of John Holdsworth's gammatone and Ray Meddis's haircell models

global sr;			% sampling rate

% synthesize & apply gammatone filter
[A, B] = MakeERBFilters(sr, nchan, flo);
yy = ERBFilterBank(A, B, y);
yy = flipud(yy);

% apply hair cell model
yyy = MeddisHairCell(yy, sr);