使用高斯/正态分布钟形曲线投掷100个骰子的概率分布

Question

我有一个程序可以模拟骰子100次。 使用此程序，我需要在直方图上显示一个钟形曲线，以显示正态/高斯分布。 在我的代码中，代表实验结果的变量是z。 钟形曲线的定义包含两个参数：ave确定曲线在z轴上必须具有最大值的位置，而sig确定峰的宽度。

 set.seed(123)

 x <- sample(1:6, size=100, replace = TRUE)

 hist(x,
 main="100 fair rolls",
 xlab = "Dice Result",
 ylab = "Probability",
 xlim=c(0.5,6.5), 
 breaks=-1:100+.5,
 prob=TRUE )

到目前为止我有高斯分布方程：

    bell<-function(z,ave,sig){
exp(-(z-ave)^2/(2*sig))/sqrt(2*pi*sig)
    }

    points<-0:600

    lines(points,bell(points,350,292))