# potremmo definire la gaussiana esplicitamente:
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, 7)
# ovvero usando il sottoprogramma già incorporato in R:
# dnorm(x, mean = ..., sd = ...)
G <- function(x) dnorm(x, mean = media, sd = sigma)
integrate(G,-Inf,10)
# ottenendo gli stessi valori