Программа для вычисления энтропии

import math
массив_симв = ["а", "6", "в", "г", "д", "е", "ж", "з", "и", "й", "к", "л", "м", "н", "о", "п", "р", "с", "т", "у", "ф", "х", "ц", "ч", "ш", "щ", "ъ", "ы", "э", "ю", "я", " " ]
массив_сочет = [] 
кол_сочет = [] 
сумма_сочет = 0
массив_Р = [] 
сумма_Р = 0
# 1) Отрывок из произведения "Мастер и Маргарита" 
file = open("1.txt") 
text = file.read() 
text = text.lower()
file.close()
 
# 2) Построить матрицу (32x32) для 2x буквенных сочетаний
# Матрица совместных вероятностей P(xi,yj)
for i in массив_симв:
    mas = []
for j in массив_симв: 
    mas.append(i+j) 
мaccив_coчет.append(mas)
 
# 3) Посчитать количество каждых сочетаний в тексте 
# # создадим пустой массив частот сочетаний 
for i in range(32): 
    mas = [] 
for j in range(32):
    mas.append(0) 
кол_сочет.append(mas)
#удалим лишние символы 
for i in '!-,._?...:;""\n':
    text = text.replace(i,'')
    #заменим символы ё на e и ь на ъ
text = text.replace('ё','е')
text = text.replace('ь','ъ')
# количество сочетаний должно быть четное количество 
kol_simv=len(text)
if kol_simv/2!=0:
    kol_simv-=1
# посчитаем все сочетания
for i in range(0, kol_simv, 2):
    ind1 = массив_симв.index(text[i])
    ind2 = массив_симв.index(text[i+1])
    кoл_coчeт[ind1][ind2] += 1 
    сумма_сочет += 1
 
# 4) Вычислим вероятности появления каждого сочетания 
# # Матрица совместных вероятностей Pi j
for i in кoл_coчeт:
    mas = [] 
    for j in i:
        p = j / сумма_сочет
        p = round(p, 4)
        mas.append(p)
        сумма_P += p
    массив_P.append(mas)
сумма_P = round(сумма_P)
 
#5) Вычислим полную условную энтропию H(Y/X) 
# H(Y/X) = -SUM(P(xi)) * SUM( P(yj/xi)*log2P(yj/xi))
 
# Вычислим P(xi)
Pxi=[]
for i in массив_P:
    px = 0
    for j in i:
        px += j
    Pxi.append(px)
 
    # Построить канальную матрицу’ P(yj/xi)
PP_yx = []
for i in range(32):
    px = []
    for j in range(32):
        x = массив_P[i][j] / Pxi[i]
        px.append(x)
    PP_yx.append(px)
 
# Вычислим частные условные энтропии
H_Yxi = []
for i in PP_yx:
    H = 0
    for j in i:
        try:
            h = j * math.log(j, 2)
        except:
            h = 0
        H += h
        H = round( H, 4)
 
    H_Yxi.append(H*(-1))
 
# Вычислим полную условную энтропию H(Y/X)
H_YX = 0
for i in range(32):
    H_YX += Pxi[i] *H_Yxi[i]
H_YX = round(H_YX, 4)
 
# 6) Вычислим энтропию двухбуквенного сочетания 
# # HfX, Y) = -СУММ( СУММ( Pij * logs(Pij)))
 
H_xy = 0
for i in массив_Р:
    for 0 in i:
        try:
            H_xy += p*math.log(p, 2)
        except:
            continue
 
H_xy = round(H_xy, 4) * (-1)
 
# 7) Вывод результатов 
print( "\n\n", "_"*26, "РЕЗУЛЬТАТЫ","_"*26)
 
print("\n\tУсловная энтропия H(Y/X) = ", H_YX)
 
print("\n\tЭнтропия H(X, Y) двухбуквенного сочетания = ", H_xy)
input('\nEnter,please')

for p in i: