N=128;T=1;t=linspace(0,T,N);y =cos(100*t);dt=t(2)-t(1);f=1/dt;clear p_saved hn_saved E_saveda=3;%for a=0:0.05:2    fy=Disfrft(y,a);    %fys=cdpei(y,a);    fyss=frft(y,a);    fysa=fft(y);    fyre=frft(fyss,-3);    fyif=ifft(y);    % blue,green,red,cyan    figure(1)    subplot(611);plot(t,real([fyss]));    title(['a = ',num2str(a)]);    subplot(612);plot(t,real([fysa]));    subplot(613);plot(t,real([y]));    subplot(614);plot(t,real([fyre]));    subplot(615);plot(t,real([fyif]));    %subplot(312);plot(x,imag([fyss]));    %subplot(313);plot(x,abs([fyss]));    %pause(0.7);  %end