import math
n = 40
C1 = 110
C2 = 10
T = [0.05, 0.07, 0.10, 0.13,0.18]

OB = []
OBD = []

for i, Pt in enumerate(T):
    OB.append(n*(C1* sum(T[:i]) + C2)/(i+1))
min_T = min(OB)
print("Очікуване значення:", min_T)
print("Індекс мінімального елементу у OB:", OB.index(min_T) + 1)

def disper2power (iteration):
  sum = 0
  for Pt in T[:iteration]:
    sum += pow(Pt,2)
  return sum


for i, Pt in enumerate(T):
    # OBD.append(C1*(n*sum(T[:i]) + 3 * math.sqrt((n-1/n)*(sum(T[:i]) - disper2power(i))))) + n*C2
    OBD.append( (C1*(n*sum(T[:i]) + (3*math.sqrt(((n-1)/n)*(sum(T[:i])-disper2power(i)) ))) + n*C2) /(i+1) )
min_D = min(OBD)
OBD.index(min_D) + 1

print("Очікуване значення дисперсії:", min_D )
print("Індекс мінімального елементу у OB:", OBD.index(min_D) + 1)