Thanks for that Onuralp. That's a good insight. I did get some good results for all the following:

stochastic_motor.visualize_attributes()
print(stochastic_motor.total_impulse)
stochastic_rocket.visualize_attributes()
stochastic_flight.visualize_attributes()

but something that stuck out to me as a potential problem was that when I tried print(stochastic_motor.total_impulse), I got back the following:

(np.float64(2221943.0652824277), np.float64(207658.2304002269), <bound method RandomState.normal of RandomState(MT19937) at 0x1D833D56840>)

Do you think this may have been the cause of the problem since it's not a single float? I had set it up like so in the script:

total_impulse=(1.07*(RedMotor.total_impulse), (RedMotor.total_impulse/10)),

type(stochastic_motor.total_impulse) returns that it is a tuple. Could this be the cause of the problem?