Displaying 2 results from an estimated 2 matches for "dotproductmatrix_m".
2012 May 23
1
numerical integration
...- function(j,s){ return(exp(-j*a*s)) }
# Inner product in $L^2([0,+\oo))$
#
innerProduct <- function(f,g){
integrand <- function(s){ return(f(s)*g(s)) }
return(integrate(integrand,lower=0,upper=Inf)$value)
}
# The 5x5 matrix of the L^2-inner products of $M_1,\dots,M_5$.
#
numericalDotProductMatrix_M <- function(){
A <- matrix(0,5,5)
for(i in seq(1:5))
for(j in seq(i:5)){
f <- function(s){ return(M(i,s)) }
g <- function(s){ return(M(j,s)) }
A[i,j] <- innerProduct(f,g)
...
2012 May 23
0
numerical integrals
...function(j,s){ return(exp(-j*a*s)) }
# Inner product in $L^2([0,+\oo))$
#
innerProduct <- function(f,g){
integrand <- function(s){ return(f(s)*g(s)) }
return(integrate(integrand,lower=0,upper=Inf)$value)
}
# The 5x5 matrix of the L^2-inner products of $M_1,\dots,M_5$.
#
numericalDotProductMatrix_M <- function(){
A <- matrix(0,5,5)
for(i in seq(1:5))
for(j in seq(i:5)){
f <- function(s){ return(M_j(i,s)) }
g <- function(s){ return(M_j(j,s)) }
A[i,j] <- innerProduct(f,g)
...