clear all
 close all
 
 %%%%%%%%%%%%% figure setup %%%%%%%%%%%%%
 lw = 2;  
 set(0, 'DefaultAxesFontSize', 15);
 fs = 13;
 msize = 5;
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
figure(1) 
xx = 0:0.01:3;  yy=[]; 
 for i=1:length(xx)
 yy =[yy Phijk(xx(i), 0, 0, MakeONFilter('Daubechies',4), 25)];
  end
  plot(xx,yy,'linewidth',lw)
  axis tight
  print -depsc 'C:\NPbook\Wavelets\daub4f.eps'

 %---------------------------------------------------------------- 
figure(6) 
xx = -1:0.01:2;  yy=[]; 
 for i=1:length(xx)
 yy =[yy Psijk(xx(i), 0, 0, MakeONFilter('Daubechies',4), 25)];
  end
  plot(xx,yy,'linewidth',lw)
    axis tight
  print -depsc 'C:\NPbook\Wavelets\daub4m.eps'

 %---------------------------------------------------------------- 
figure(2) 
xx = 0:0.01:7;  yy=[]; 
 for i=1:length(xx)
 yy =[yy Phijk(xx(i), 0, 0, MakeONFilter('Daubechies',8), 25)];
  end
  plot(xx,yy,'linewidth',lw)
    axis tight
  print -depsc 'C:\NPbook\Wavelets\daub8f.eps'

 %---------------------------------------------------------------- 
figure(3) 
xx = -3:0.01:4;  yy=[]; 
 for i=1:length(xx)
 yy =[yy Psijk(xx(i), 0, 0, MakeONFilter('Daubechies',8), 25)];
  end
  plot(xx,yy,'linewidth',lw)
   axis tight
  print -depsc 'C:\NPbook\Wavelets\daub8m.eps'
 %---------------------------------------------------------------- 
figure(4) 
xx = 0:0.01:15;  yy=[]; 
 for i=1:length(xx)
 yy =[yy Phijk(xx(i), 0, 0, MakeONFilter('Daubechies',16), 25)];
  end
  plot(xx,yy,'linewidth',lw)
    axis tight
  print -depsc 'C:\NPbook\Wavelets\daub16f.eps'

 %---------------------------------------------------------------- 
figure(5) 
xx = -7:0.01:8;  yy=[]; 
 for i=1:length(xx)
 yy =[yy Psijk(xx(i), 0, 0, MakeONFilter('Daubechies',16), 25)];
  end
  plot(xx,yy,'linewidth',lw)
    axis tight
  print  -depsc 'C:\NPbook\Wavelets\daub16f.eps'

 %----------------------------------------------------------------