Per realizzare l'animazione che illustra l'integrazione di √(1 - x²) 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