Histogramme[Séquence[(i - n p) / s - 0.5 / s, i, 0, n + 1], Séquence[s Combinaison[n, k] p^k (1 - p)^(n - k), k, 0, n]]

b = 2.1 (slider from -5 to 5, step 0.01)
a = -2 (slider from -5 to 5, step 0.01)
n = 550 (slider from 10 to 1000, step 10)
p = 0.71 (slider from 0 to 1, step 0.01)
s: sqrt(n * p * (1 - p))
e = 2.718281828459045
f(x) = 1 / sqrt(2 * 3.141592653589793) * e^((-x^(2)) / 2)
a_1 = Histogram(Sequence[(i - n * p) / s - 0.5 / s, i, 0, n + 1], Sequence[s * Binomial[n, k] * p^(k) * (1 - p)^(n - k), k, 0, n])
B: (a, f(a))
D: (b, f(b))
A: (a, 0)
b_1 = Segment(A, B)
C: (b, 0)
c = Segment(C, D)
c_1 = 0
c_2 = 0
aire = Integral(f, a, b)
c_3 = 0
d = BinomialDist(n, p, floor(b * sqrt(n * p * (1 - p)) + n * p), true)
g = BinomialDist(n, p, floor(a * sqrt(n * p * (1 - p)) + n * p), true)
bin: d - g