function muscwt(file,mystart,docwt)
%logistic1 48000 pd=3 46000 intermittency
%muscwt('logistic2.wav',floor((3.328-2.5)/(4-2.5)*2^16),1)
gencos=0;
%docwt=1;
playsnd=1;
%mystart=1;
factor=10;
epoch=1024;
harm=3;
lps=1;
samprate=8000;
if gencos==0
    %y=wavread('soundw.wav');
    y=wavread(file);
else
    mysiz=40*samprate;
    y=zeros(1,mysiz);
    for j=1:lps
        for i=1:mysiz
            y(i)=y(i)+1/lps*cos(2*pi*i/factor);           
        end
    factor=factor/harm;
    end
end
%y=rand(1,length(y));
%myvec=linspace(1,8*64,64);
myvec=logspace(0.1,3.2,64);
if docwt==1
    c=cwt2sq(y(mystart:mystart+epoch),myvec,'morl','absglb',samprate);
else
    c=fft(y(mystart:mystart+epoch));
    x=linspace(1,epoch+1,epoch+1);
    plot(x,max(imag(c))/2*y(mystart:mystart+epoch)/max(y(mystart:mystart+epoch)),'y');
    hold on
    plot(x,imag(c));
    plot(x,real(c),'r');
    hold off
end
    if playsnd
    player=audioplayer(y(mystart:mystart+epoch),samprate);
    %player=audioplayer(y,40100);
    play(player);
end