% Daubechies Lagarias for Pollen2 wavelets
% ----------------------------------------
xxs = 0:0.005:6;
xxw = -2:0.005:4;
 yys =[];
 yyw =[];
 poll2 = MakePollen2(0.67890, 0.12345);
for i = 1:length(xxs)
    yys = [yys Phijk(xxs(i), 0, 0, poll2, 25)];
    yyw = [yyw Psijk(xxw(i), 0, 0, poll2, 25)];
end
figure(1)
plot(xxs, yys)
title('Pollen2 scaling function for phi1=0.6789 and phi2=0.12345')
print -deps pollens.eps 
figure(2)
plot(xxw, yyw)
title('Pollen2 wavelet function for phi1=0.6789 and phi2=0.12345')
print -deps pollenf.eps