For a large sample, a confidence interval for the mean is with sample std (ddof=1).
Implement:
normal_quantile(p): the standard normal quantile, found by bisection on the CDF over to within 1e-12mean_ci(x, conf=0.95): returns (low, high) using normal_quantile(1 - (1 - conf) / 2)