clear all
haar1d([1,2,3,4])
x=[1:2^5]*2*pi/2^5;
y=sin(x);
hold off
plot(x,y)
hold on
yh=haar1d(y);
plot(x,yh,'r')
