Monte Carlo failure probability with an error bar: does the shaft fit?

Afterwards you can estimate a failure probability by sampling with NumPy, give it a standard error, and say how many samples a target precision needs.