Displaying 2 results from an estimated 2 matches for "testnumericalintegration".
2012 May 23
1
numerical integration
...A[j,i] <- A[i,j]
}
return(A)
}
# The 5x5 matrix of the L^2-inner products of $M_1,\dots,M_5$.
#
dotProductMatrix_M <- function(){
A <- matrix(0,5,5)
for(i in seq(1:5))
for(j in seq(1:5)) A[i,j] <- 1/(a*(i+j))
return(A)
}
testNumericalIntegration <- function(){
A <- dotProductMatrix_M()
print("Dot products, analytic:")
print(A)
# why doesn't the numerical integration work
B <- numericalDotProductMatrix_M()
print("Dot products, numeric:")
print(B)
}
testNumericalIntegratio...
2012 May 23
0
numerical integrals
...A[j,i] <- A[i,j]
}
return(A)
}
# The 5x5 matrix of the L^2-inner products of $M_1,\dots,M_5$.
#
dotProductMatrix_M <- function(){
A <- matrix(0,5,5)
for(i in seq(1:5))
for(j in seq(1:5)) A[i,j] <- 1/(a*(i+j))
return(A)
}
testNumericalIntegration <- function(){
A <- dotProductMatrix_M()
print("Dot products, analytic:")
print(A)
# why doesn't the numerical integration work
B <- numericalDotProductMatrix_M()
print("Dot products, numeric:")
print(B)
}
testNumericalIntegrat...