Create an rng object with np.random.default_rng(), you can seed it for reproducible results. You can draw samples from probability distributions, including from the binomial and normal distributions.

In a recent write-up, [David Delony] explains how he built a Wolfram Mathematica-like engine with Python. Core to the system is SymPy for symbolic math support. [David] said being able to work with ...