Per realizzare l'animazione che illustra l'integrazione di √(1 - x2) si può usare il software R con i comandi seguenti. Volendo si possono 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^14) # isto(2^15) pi/2