With the premise that social interactions are described by power-law distributions, we study the stochastic dynamics of SIR (Susceptible-Infected-Removed) compartmental models on static scale-free random networks generated via the configuration model. We compare simulations of our model to analytical results, providing a closed formula and a lower bound for the probability of having a minor epidemic of the disease. We explore the variability in disease spread by stochastic simulations. In particular, we demonstrate how important epidemic indices change as a function of the contagiousness of the disease and the connectivity of the network. Our results quantify the role of the starting node's degree in determining these indices, commonly used to describe epidemic spread. Our results and implementation set a baseline for studying epidemic spread on networks, showing how analytical methods can help in the interpretation of stochastic simulations.