Выбор максимального значения среди всех возможных реализаций формулы

def e_j(сatеgory_stab):
    global e
    if сatеgory_stab == 'A':
        e = 0.05
    elif сatеgory_stab == 'B':
        e = 0.06
    elif сatеgory_stab == 'C':
        e = 0.06
    elif сatеgory_stab == 'D':
        e = 0.12
    elif сatеgory_stab == 'E':
        e = 0.34
    elif сatеgory_stab == 'F':
        e = 0.53
    elif сatеgory_stab == 'G':
        e = 0.53
    return e

def met_factor():
    stability = ['A', 'B', 'C', 'D', 'E', 'F', 'G']
    periodicity = [0.08, 0.09, 0.10, 0.10, 0.12, 0.21, 0.17, 0.13]
    N = 8
    h = 120
    factor_G = []
    sch1 = 0
    while sch1 < len(stability):
        сatеgory_stab = stability[sch1]
        sch2 = 0
        while sch2 < len(periodicity):
            povtor = periodicity[sch2]
            G = (1 - k_b()) * ((2 * N * povtor) / (((2 * math.pi) ** (3.0 / 2.0)) * x)) * (
                f_istosh() / (dispersion() * modul_skor())) * math.exp(
                - ((h + podyem_strui()) ** 2) / (2 * dispersion() ** 2)) + k_b() * (
                (2 * N * povtor) / (math.pow((2 * math.pi), 3.0 / 2.0))) * (f_ist_rad() / (
                dispersion() * modul_skor() * x)) * math.exp(- (h / (2 * dispersion())))
            factor_G.append(G)
            sch2 += 1
        sch1 += 1
    return my_max(factor_G)