By (Crandall, Projects in Scientific Computation - p 202) the
wavelet transform of a function f is :
We can approximate this for values of
A = h/3 * [ y 0 + 4y 1 + 2 y 2 + 4 y 3 + ... + y n ]
This should have n even.