f(x)=\dfrac{1}{\sqrt{2\pi}}e^{-\frac{x^2}{2}}

n = 5 (slider from 5 to 5000, step 1)
m: (n * 0.4)
σ: sqrt(((n * 0.4) * 0.6))
H = Histogram(Sequence[(i - m) / σ - 0.5 / σ, i, 0, n + 1], Sequence[(σ * BinomialDist[n, 0.4, k, false]), k, 0, n])
f(x) = (1 / sqrt((2 * π)) * exp(((-x^(2))) / 2))
Gauss = Integral(f, -1, 2)
a: m - σ
b: m + (2 * σ)
P: BinomialDist[n, 0.4, floor(b), true] - BinomialDist[n, 0.4, floor(a), true]
d: Gauss - P