Symlet wavelet filter computation

The `symaux`

function generates the scaling filter coefficients
for the "least asymmetric" Daubechies wavelets.

is the
order `w`

= symaux(`n`

)`n`

Symlet scaling filter such that ```
sum(w) =
1
```

.

Instability may occur when

`n`

is too large. Starting with values of`n`

in the 30s range, function output will no longer accurately represent scaling filter coefficients.As

`n`

increases, the time required to compute the filter coefficients rapidly grows.For

`n`

= 1, 2, and 3, the order`n`

Symlet filters and order`n`

Daubechies filters are identical. See Extremal Phase.

[1] Daubechies, I. *Ten
Lectures on Wavelets*, CBMS-NSF Regional Conference Series in Applied Mathematics.
Philadelphia, PA: SIAM Ed, 1992.

[2] Oppenheim, Alan V., and Ronald W.
Schafer. *Discrete-Time Signal Processing*. Englewood Cliffs, NJ:
Prentice Hall, 1989.