③ Why it works

For a normal curve, probability is area. Converting each boundary to a z-score standardizes the problem, and the cumulative function Φ(z) — the total area left of z — is computed with a classic polynomial approximation (Abramowitz & Stegun 7.1.26) accurate to about 7 decimal places. Subtracting two cumulative areas gives the area trapped between two values.