# 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