function band=spread_hz(hz_array,hz_c)
%hz_array is an array of hz frequencies
%hz_c is the current centre frequency in Hz
barkc=f2bark(hz_c);
band=zeros(size(hz_array));
%---------------
for hi=1:length(hz_array)
    barki=f2bark(hz_array(hi));
    barkd=barki-barkc;
    if barkd >= -2.5 & barkd <=-0.5
        band(hi)=10^(1*(barkd+0.5));
    elseif barkd > -0.5 & barkd <0.5
        band(hi)=1;
    elseif barkd >= 0.5 & barkd <=1.3
        band(hi)=10^-(2.5*(barkd-0.5));
    end
    %---    
   % fprintf('index=%d, bark=%.2f, freq=%.2f, band=%d, barkd=%f, r=%d\n',hi,barki,bark2f(barki),band(hi),barkd,region);
end
