"""distribusi_lengkap.py — tiga distribusi inti dalam satu skrip."""

from scipy.stats import binom, poisson, norm

print(f"binomial : P(X = 2) = {binom.pmf(2, 5, 0.5):.4f} (puncak bersama k = 3)")  # 5 lemparan koin
print(f"poisson  : P(X <= 2) = {poisson.cdf(2, 3):.4f}")  # call center: 3 panggilan per jam

tinggi = norm(loc=168, scale=7)  # tinggi badan: mu = 168 cm, sigma = 7 cm
print(f"normal   : P(X <= 175) = {tinggi.cdf(175):.4f}")
print(f"normal   : persentil ke-90 = {tinggi.ppf(0.90):.2f} cm")