function wcshow(wc,L)

[N,J] = quadlength(wc);

for j=L:J-1
  low = 1:2^j;
  high = 2^j+1:2^(j+1);
  scale = 3*2^(j-L+1);
  wc(high,high) = scale*wc(high,high);
  wc(high,low)  = scale*wc(high,low);
  wc(low,high)  = scale*wc(low,high); 
end

low = 1:2^L;
wc(low,low) = 2^((L-J))*wc(low,low);

wc=abs(wc);
imshow(wc,[0 1]);
