[LO_D,HI_D,LO_R,HI_R] = wfilters('db1')
h=LO_R
g=HI_R
h(1)=h(1)+0.01
%g(2)=g(2)-0.001
iterations=10
%initsignal=[1]
initsignal=[0 0 0 1 0 1 0 0 0 0 1]

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