Numeric Integration Using Simpson's method
Numeric Integration using the Simpson Method
simp(y, a = NULL, b = NULL, x = NULL, n = 200)
y
: y values to integratex
: x values to integrate overa
: NULLb
: NULLn
: number of divisions defaults to 200Numeric Value of the integration
# 4-x^2-y^2 fun <- function(x, y){ a <- 4 b <- x^2 d <- y^2 z <- a-b-d return(z) } a <- fun(seq(-1000,1000,1), seq(-1000,1000,1)) simp(a, x=-1000:1000, n=1000)
Rolf Turner