There are several mathematical models that can replicate some of the
statistically random aspects of evolution. Specifically, there are several
theoretical models for mutation of a DNA or protein sequence, and empirical
models of how a particular feature has evolved (e.g., the number of spines
per whorl on some snails shows an increase in variance but no significant
trend over time and the size of forams may increase, decrease, or stay the
same over time). There are many models in population genetics for random
changes in gene frequency. These by no means approach a model that can go
from astrophysics to life, but do model the statistical randomness of
evolution.