Wavelet-Galerkin method