simp function

Numeric Integration Using Simpson's method

Numeric Integration using the Simpson Method

simp(y, a = NULL, b = NULL, x = NULL, n = 200)

Arguments

• y: y values to integrate
• x: x values to integrate over
• a: NULL
• b: NULL
• n: number of divisions defaults to 200

Returns

Numeric Value of the integration

Examples

# 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)

Author(s)

Rolf Turner