function F=daubfdwt(c0,c1,c2,c3,f)
%fprintf('here\n')
%fprintf('N=%d\n',size(f,1))
N=size(f,1);
%k=log10(N)/log10(2);
if N>4 
%lpf(c0,c1,c2,c3,N)*f
%subplot(4,2,10-k)
%plot(lpf(c0,c1,c2,c3,N)*f)  
%t = sprintf('Level %d',10-k);
%title(t)
F=[daubfdwt(c0,c1,c2,c3,lpf(c0,c1,c2,c3,N)*f);hpf(c0,c1,c2,c3,N)*f];

else
%lpf(c0,c1,c2,c3,N)*f
%subplot(4,2,10-k)
%plot(lpf(c0,c1,c2,c3,N)*f)
%t = sprintf('Level %d',10-k);
%title(t)
F=[lpf(c0,c1,c2,c3,N)*f;hpf(c0,c1,c2,c3,N)*f;F];
end

