在 R 中硬编码双黎曼和
Hard Coding a Double Reimann Sum in R
我正在尝试在 R 中为函数 f(x,y) = x^5 * y^2 编写双 Reimann 和代码，其中 0 < x < 2 和 4 < y < 6。这就是到目前为止，我很难找到正确的答案，因为这段代码输出的值不正确。
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
f <- outer(x, y, function(x, y) x^5 * y^2)
ans <- sum(f)*xsize*ysize
ans
#> [1] 540.4438
尽管显然在这种情况下您可以独立地计算 x 和 y 的总和而不必计算外积。
我正在尝试在 R 中为函数 f(x,y) = x^5 * y^2 编写双 Reimann 和代码，其中 0 < x < 2 和 4 < y < 6。这就是到目前为止，我很难找到正确的答案，因为这段代码输出的值不正确。
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
f <- outer(x, y, function(x, y) x^5 * y^2)
ans <- sum(f)*xsize*ysize
ans
#> [1] 540.4438
尽管显然在这种情况下您可以独立地计算 x 和 y 的总和而不必计算外积。