% hybridbasis
p = 7;

rows = floor(sqrt(p)); % how many plot rows
cols = ceil(p/rows); % how many plot columns

% levels and shifts
LS = [([-1:1:p-1]'),zeros(p+1,1)];

for i = 1:(rows*cols),
  lev = LS(i,1); % level
  sh = LS(i,2); % shift
  
  if(lev == -1),
    w = [1, zeros(1, 2^p - 1)];
  else,
    w = putlevel(zeros(1, 2^p), ([0:(2^lev - 1)]==sh), lev);
  end
  
  f3 = idpwt(w(1:2^3),daubcofs(1),0);
  f5 = idpwt([f3, w((2^3 + 1):2^5)], daubcofs(2),3);
  f = idpwt([f5, w((2^5 + 1): 2^p)], daubcofs(8),5);
  m = max(abs(f));
  ff = f/m;
  subplot(rows,cols,i);
  plotperdisc(ff);
  %title(sprintf('level %d',i-2));
  drawnow
end;