function [axis, cut]=plotcut(time,freq,cutoff,step)
%
% plot the optimal rate vs frequency cutoff
%

cut=zeros(ceil(cutoff/step),1);
axis=cut;
time2=time;

j=1;
for i=1:step:cutoff
time2(find(freq<i))=0;
axis(j)=i;
cut(j)=rate1(time2,findqopt1(time2));
j=j+1;
end;

plot(axis,cut);

