search for: innerproduct

...2.500000 2.083333 0.000000 0.000000 [4,] 2.500000 2.083333 0.000000 0.000000 0.000000 [5,] 2.083333 1.785714 1.562500 1.388889 1.250000 Undoubtedly the code is suboptimal. What is wrong with it? a  = 0.08  # alpha M <- 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:...

...500000 2.083333 0.000000 0.000000 [4,] 2.500000 2.083333 0.000000 0.000000 0.000000 [5,] 2.083333 1.785714 1.562500 1.388889 1.250000     Undoubtedly the code is suboptimal. What is wrong with it?   a  = 0.08  # alpha M <- 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...