sqwd function

Compute expansion with respect to squared wavelets.