F <- function(x) 1/sqrt(2*pi)*exp(-x^2/2) G <- function(x) F((x-media)/sigma)/sigma # definisci media e sigma, quindi scegli a e b e calcola integrate(G,a,b) # Esempio: media <- 7; sigma <- 3; integrate(G, -Inf, 10)