function b = hm(a,f0,fc)
% hm calculates response of Colburn filters to amplitude a
% a represents stimulus amplitude, f is stimulus frequency,
% fc is characteristic frequency

alpha = 4;
temp = size(f0);
numcomp = temp(2);
b = zeros(1,numcomp);

for i = 1:numcomp,
f = f0(i);

if f <= fc
	b(i) = a(i)*(f/fc)^alpha;

elseif f > fc
	b(i) = a(i)*(fc/f)^(alpha*2);

end	
end
