function [Y] = chopchop(fmin,fmax,F,Y,degr)


Y(F>=fmin & F<=fmax)=NaN;
Y2=polyval(polyfit(F(~isnan(Y)),Y(~isnan(Y)),degr),F);
Y(F>=fmin & F<=fmax)=Y2(F>=fmin & F<=fmax);

