Implementation of Trapezoidal Rule (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.33336533334409
# ? 2018 TheFlyingKeyboard and released under MIT License

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.5 * (f(intervalBegin) + f(intervalEnd))

for(i in 0:iterations){
  integral = integral + f(intervalBegin + step * i)

integral = integral * step

print(paste("Integral is equal to: ", integral))


R Trapezoidal Rule
Tagged on:             

Leave a Reply

Your email address will not be published. Required fields are marked *

By continuing to use the site, you agree to the use of cookies. You can read more about it the Cookies&Privacy Policy Section Above. more information

The cookie settings on this website are set to "allow cookies" to give you the best browsing experience possible. If you continue to use this website without changing your cookie settings or you click "Accept" below then you are consenting to this. You can read more about it the Cookies&Privacy Policy Section.