P(n,p,k<=r)

n = 50 (slider from 0 to 2000, step 1)
p = 0.67 (slider from 0 to 1, step 0.01)
mu: n * p
sigma: (n * p * (1.0 - p))^(0.5)
A: (mu, 0.0)
B: (mu, 0.4)
b = Segment(A, B)
r = 30 (slider from 0 to 600, step 1)
s = 35 (slider from 0 to 500, step 1)
t = 4000 (slider from 1 to 4000, step 1)
φ(x) = 1.0 / (2.0 * 3.141592653589793)^(0.5) * exp((-x^(2.0)) / 2.0)
c = BarChart((-0.5 - mu) / sigma, (r + 1.0 - 0.5 - mu) / sigma, Sequence[Binomial[n, k] * p^(k) * (1.0 - p)^(n - k) * sigma, k, 0.0, r, 1.0])
a = BarChart((-0.5 - mu) / sigma, (n + 1.0 - 0.5 - mu) / sigma, Sequence[Binomial[n, k] * p^(k) * (1.0 - p)^(n - k) * sigma, k, 0.0, n, 1.0])
d = BarChart((s - 0.5 - mu) / sigma, (n + 0.5 - mu) / sigma, Sequence[Binomial[n, k] * p^(k) * (1.0 - p)^(n - k) * sigma, k, s, n, 1.0])
Φu = LowerSum(φ, -5.0, (r - mu + 0.5) / sigma, t)
Φo = UpperSum(φ, -5.0, (r - mu + 0.5) / sigma, t)