Implementation of Midpoint Rule (Rectangle Method) (Integral Approximation) written in R.
Enter beginning of interval: -2 Enter end of interval: 2 Enter number of iterations: 1000000 Integral is equal to: 5.33333333334422
# ? 2018 TheFlyingKeyboard and released under MIT License # theflyingkeyboard.net f <- function(x) { return(x * x) } intervalBegin <- as.double(readline(prompt = "Enter beginning of interval: ")) intervalEnd <- as.double(readline(prompt = "Enter end of interval: ")) iterations <- as.integer(readline(prompt = "Enter number of iterations: ")) step = (intervalEnd - intervalBegin) / iterations integral = 0 for (i in 1:iterations) { integral = integral + step * f(intervalBegin + (i - 1) * step) } print(paste("Integral is equal to: ", integral))
R Midpoint Rule (Rectangle Method)