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)