% daubbasis.m -- Plot basis functions for Daubechies wavelet.

N = 128;	% Number of points -- must be a power of 2.

daubinit;

for n=1:N
  X = [zeros(n-1,1); 1; zeros(N-n,1)]';
  x = daubdwt(c0, c1, c2, c3, X);
  plot(x)
  fprintf('Function %d of %d -- ENTER to continue ', n, N);
  pause
  fprintf('\n')
end
