LLMs can use Python to simulate from probability distributions. Though, admittedly they have to code and use their own MCMC samplers (and can’t yet utilize Stan and PyMC directly).