# Integrale di sqrt(1-x^2) tra -1 ed 1 (semicerchio). Puoi # cambiare la funzione e l'intervallo. f <- function(x) sqrt(1-x^2) isto <- function(n) { plot(f,-1,1,col="blue",asp=1) i <- 0:(n-1); k <- -1+1/n+i/(n/2); h <- k-1/n j <- 1:n; kk <- -1+1/n+j/(n/2); hh <- kk-1/n rect(h,0,hh,f(k)); print(sum(2/n*f(k))) } isto(2) # isto(2^2) # isto(2^3) # isto(2^4) # 1.578343 # ... isto(2^15) # 1.570796 pi/2 # 1.570796 # se interessa solo il calcolo usa: integra <- function(n) { i <- 0:(n-1); k <- -1+1/n+i/(n/2); sum(2/n*f(k)) } integra(2^15) # 1.570796