Binomial distribution with `trials` (n ≥ 0, ≤ 1e6) and success `probability` (0..1). `function` selects the mass (pmf), cumulative distribution (cdf) or quantile / inverse-CDF.