Compute probability density, cumulative probability, interval probabilities, and inverse CDF (quantiles) for any Gaussian distribution. Visualize the bell curve with dynamic shading.
The normal distribution, also known as the Gaussian distribution, is the most important continuous probability distribution in statistics. Its probability density function (PDF) forms the symmetric bell curve defined by the mean (μ) and standard deviation (σ). The formula is:
$$f(x) = \frac{1}{\sigma\sqrt{2\pi}} e^{-\frac{1}{2}\left(\frac{x-\mu}{\sigma}\right)^2}$$
The cumulative distribution function (CDF) gives the probability that a random variable X is less than or equal to x, while the inverse CDF (quantile function) finds the x-value for a given cumulative probability. This calculator uses high-precision approximations (error < 7.5×10⁻⁸) based on the Abramowitz & Stegun algorithm, widely adopted in scientific computing.
Algorithm Accuracy Comparison — Validated against reference values from R (pnorm/qnorm) and SciPy:
| Test Case | This Calculator | Reference (R / SciPy) | Difference |
|---|---|---|---|
| Φ(1.96) [standard normal] | 0.975002 | 0.975002 | < 5×10⁻⁷ |
| Φ⁻¹(0.975) | 1.959964 | 1.959964 | < 1×10⁻⁷ |
| N(100,15): P(X ≤ 130) | 0.977250 | 0.977250 | < 1×10⁻⁶ |
| N(50,2): 95th percentile | 53.289707 | 53.289707 | < 1×10⁻⁶ |
All values rounded to 6 decimals. Maximum absolute error < 7.5×10⁻⁸ for CDF, < 1×10⁻⁹ for inverse CDF (Acklam method).
The standard normal CDF Φ(z) is computed using the rational approximation from Abramowitz and Stegun (26.2.17). For a general normal N(μ,σ), we transform: z = (x - μ)/σ, then CDF = Φ(z). The inverse CDF (quantile) is obtained via the Acklam algorithm, which is accurate to nearly machine precision. The PDF is computed directly from the closed form. Our implementation ensures numerical stability even for extreme tails.
The empirical rule (68‑95‑99.7) emerges naturally: about 68.27% of values lie within μ±σ, 95.45% within μ±2σ, and 99.73% within μ±3σ. This calculator also shows these ranges as reference.
Important note on assumptions: The normal distribution assumes continuous, symmetric data without heavy tails. For real datasets, always verify normality assumptions (e.g., using Q-Q plots or Shapiro-Wilk test) before applying parametric methods. This calculator is intended for theoretical calculations and educational purposes.
A manufacturing process produces shafts with target diameter 50 mm and σ = 0.2 mm. Using the calculator with μ=50, σ=0.2, the probability that a shaft lies within the tolerance [49.8, 50.2] is P(49.8 ≤ X ≤ 50.2) = ? Enter bounds 49.8 and 50.2 → probability ≈ 0.6827 (68.27%). For stricter limits [49.7, 50.3] (3σ), the probability rises to 0.9973. This directly quantifies process capability.
Assume daily returns of a stock follow a normal distribution with μ = 0.05% and σ = 1.2%. A risk manager wants the 5% Value at Risk (VaR) – the loss level that will not be exceeded with 95% confidence. Using the inverse CDF with p = 0.05, the calculator returns x = μ + σ·Φ⁻¹(0.05) ≈ -1.92%. Thus, there is a 5% chance of a daily loss exceeding 1.92%.
Clinical laboratories often establish reference intervals assuming normality. For a blood test with μ=120 mg/dL and σ=15 mg/dL, the central 95% reference interval is [μ - 1.96σ, μ + 1.96σ] = [90.6, 149.4] mg/dL. This calculator instantly provides these critical values via the quantile function (p=0.025 and p=0.975).