在 R 中对双黎曼和进行硬编码

Question

我正在尝试在 R 中为函数 f(x,y) = x^5 * y^2 编写双 Reimann 和，其中 0

a <- 0 # lower bound for x
b <- 2 # upper bound for x
c <- 4 # lower bound for y
d <- 6 # upper bound for y
xsize <- (b-a)/1000 # step size in x-direction
ysize <- (d-c)/1000 # step size in y-direction
x <- seq(a + (xsize/2), b - (xsize/2), xsize)
y <- seq(c + (ysize/2), d - (ysize/2), ysize)

f <- x^5 * y^2

ans <- sum(f*xsize*ysize)
ans

Answer 1

f <- outer(x, y, function(x, y) x^5 * y^2)
ans <- sum(f)*xsize*ysize
ans
#> [1] 540.4438

虽然显然在这种情况下，您可以独立地计算 x 和 y 的总和，而不必计算外积。