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