[LO_D,HI_D,LO_R,HI_R] = WFILTERS('db2')
h=LO_R
g=HI_R
%h(1)=h(1)+0.5
%g(2)=g(2)-0.5
iterations=8
initsignal=[0 1 0]

for a=1:iterations
   phi=upcoef('a',initsignal,h,g,a);
   psi=upcoef('d',initsignal,h,g,a);
   subplot(2,1,1)
   plot(phi)
   subplot(2,1,2)
   plot(psi);
 pause
end