# Example:
#> Daub2
#[1]  0.4829629  0.8365163  0.2241439 -0.1294095
#> T0(Daub2)
#          [,1]       [,2]      [,3] 
#[1,] 0.6830127  0.0000000 0.0000000
#[2,] 0.3169873  1.1830127 0.6830127
#[3,] 0.0000000 -0.1830127 0.3169873
#
 T0  <- function(filter)
{
        N <- length(filter)
        domain <- N - 1
        T0 <- matrix(0, domain, domain)
        #       cat("\n Filter is", filter, "\n Domain and size of T0 is", domain, "\n"
#               )
        for(i in 1:domain) {
                for(j in 1:domain) {
                        if(2 * i - j > N || 2 * i - j < 1)
                                next
                        T0[i, j] <- sqrt(2) * filter[2 * i - j]
                }
        }
        return(T0)
}