Нахождение вероятности небезопасного моделирования лифта

Question

Я могу найти это с помощью Excel, но я хотел бы знать, как это сделать с помощью Python.

1. Каждый столбец в столбце Person генерируется случайным образом с нормальным распределением стандартного отклонения 15 и среднего значения 80.
2. Количество людей случайным образом генерируется от 1 до 6.
3. Общий вес - это сумма человека с количеством людей. например первый - сумма первых 4-х человек.
4. Небезопасно - логическое значение. Если общий вес больше 500, тогда 1 в противном случае 0.
5. Испытаний - 3997. Сумма небезопасных - 148. Затем вероятность рассчитывается числом, превышающим вес, деленным на испытания.

Я могу создать нормальный дистрибутив, используя numpy. Я создаю случайное число со средним значением 80 и 15 с.д. и случайное число от 0 до 7 (от 1 до 6).

import numpy as np
from decimal import *
total_num = 0
i=1
trial = 100
while i < trial:
    rand_num =np.random.randint(1,7)
    person = np.random.normal(80, 15)
    total_weight = int(rand_num)*Decimal(person)
    # print(total_weight)
    if total_weight > int(500):
        total_num += total_num
    i += 1;
total_num/trial

Но я получаю 0 в моем выводе. Что я тут не так делаю?

Answer 1

Мне просто нужно изменить total_num += 1.

import numpy as np

count = 0
i=int(1)
trial = int(10000)
limit = int(500)
while i < trial:
    rand_num =np.random.randint(1,7)
    person = np.random.normal(80, 15)
    total_weight = int(rand_num)*person
    # print(total_weight)
    if total_weight > limit:
        #print('yes')
        count += 1
        count
    i += 1;
print("Probability is: ")
print(count/trial)