Hmmm... non mi è ancora del tutto chiaro il problema... (es. cos'è nella formula la sommatoria va da 0 a hu, cos'è hu?)
In ogni caso, l'implementazione della formula in R è triviale, ti faccio un esempio di come ottenere l'istogramma. Consiglio di usare la funzione replicate, che semplifica questo tipo di procedure.
Esempio banale: vogliamo calcolare l'istogramma dei risultati del lancio di n dadi
somma.dadi <- function(num.dadi)
{
# "Tiriamo il dado" con la funzione sample che, in questo caso sceglie
# n numeri fra 1 e 6.
return (sum(sample(1:6, num.dadi, replace=T)))
}
n.dadi <- 2
# Chiamiamo la funzione 10000 volte e mettiamo il risultato nel vettore res
res <- replicate(10000, somma.dadi(n.dadi))
# Creiamo l'istogramma
hist(res, 1:(6*n.dadi), col="black", las=1, freq=F)
Come atteso la distribuzione del lancio di 2 dadi sarà una gaussiana centrata sul 6