from math import pi
from math import log10

L_wi = [
               [68.0,  70.0,  73.0,  79.0,  81.0,  82.0,  80.0,  73.0],
               [101.0, 102.0, 100.0, 101.0, 99.0,  99.0,  97.0,  95.0],
               [90.0,  91.0,  98.0,  99.0,  97.0,  93.0,  91.0,  86.0]
]

X = 1.0; F = 1.0

Omega = [ 2*pi, pi/2, 2*pi ]

R = [2.0, 8.0, 8.0]

n = 3; f_count = len(L_wi[0])

M = [0.65, 0.62, 0.64, 0.75, 1.0, 1.5, 2.4, 4.2]

B = [32.5, 30.1, 32.0, 37.5, 50.0, 75.0, 120.0, 210.0]

K = [0.72, 0.70, 0.99, 0.99, 0.97, 0.9, 0.83, 0.76]

UDZ_dop = [91, 83, 77, 73, 70, 68, 66, 64]

UZD = []; DELTA_UZD = []


for j in range(f_count):
    sum1 = 0.0
    sum2 = 0.0
    for i in range(n):
        sum1 += (10.0 ** (0.1 * L_wi[i][j]) * X * F) / (Omega[i] * R[i] ** 2.0)
        sum2 += 10.0 ** (0.1 * L_wi[i][j])
    UZD_j = round(10.0 * log10(sum1 + 4.0 / (K[j] * B[j]) * sum2), 2)
    UZD.append(UZD_j)
    DELTA_UZD.append(round(UZD_j - UDZ_dop[j], 2))

print("UZD:     " + str(UZD))
print("UZD_dop: " + str(UDZ_dop))
print("-" * 80)
print("DELTA:   " + str(DELTA_UZD))