% Thursday 10/7
 %------------Example 1 -------------------------------------
        clf;
        n = 256;
        wplist = [4 0 2; 4 1 2; 4 2 2; 4 4 2; 4 6 5;  4 8 3;  4 10 5; 4 13 6; 4 15 0];
        sz = size(wplist); nr = sz(1);
        for i=1:nr,

                d = wplist(i,1); %<-- max depth is log_2(n)
                b = wplist(i,2); %<-- oscillation index [0 2^d - 1]
                k = wplist(i,3); %<-- shift

                wp = MakeWaveletPacket(d,b,k,'Daubechies',4,n);
                fwp = abs(fft(wp));
                t = (0:(n-1)) ./n;
                wpstr = sprintf('WaveletPacket(%i,%i,%i)',d,b,k)

                figure(i);
                set(0, 'DefaultAxesFontSize', 12);
                set(0, 'DefaultAxesFontName', 'Times')

                subplot(3,1,1);

                plot(fwp(1:(n/2)),(1:(n/2))./(n/2))
                xlabel(['|FFT(' wpstr ')|']);
                ylabel('Frequency');

                subplot(3,1,2);

                plot(t,wp)
                ylabel(wpstr);
                xlabel('Time');

                subplot(3,1,3);
                qmf = MakeONFilter('Daubechies',4);
                ImageAtomicPhase('WP', [1 d b k ], n, wpstr, 256,qmf);

                print(gcf,'-deps',['hand1' int2str(i) '.eps'])
        end