Matstat

import math

f = open("input.txt")
s = f.read()
s = s.replace("\n", "\t")
a = s.split("\t")
for i in range(len(a)):
    a[i] = int(a[i])
a.sort()

n = len(a)
frq = [1]
digit = [a[0]]
for i in range(1, len(a)):
    if a[i] != a[i-1]:
        digit.append(a[i])
        frq.append(1)
    else:
        frq[len(frq) - 1] += 1

print(digit)
print(frq)

N = len(digit)
p = []
for i in range(N):
    p.append(frq[i]/n)

b = [[0], [0], [0], [0]]
for i in range(N):
    b[0][0] += digit[i] * frq[i]
    b[1][0] += digit[i]**2 * frq[i]
    b[0].append(digit[i] * frq[i])
    b[1].append(digit[i]**2 * frq[i])



x = b[0][0] / n
S = math.sqrt(b[1][0] * x**2 / n)

for i in range(N):
    b[2][0] += math.pow(digit[i] - x, 3) * frq[i]
    b[3][0] += math.pow(digit[i] - x, 4) * frq[i]
    b[2].append(math.pow(digit[i] - x, 3) * frq[i])
    b[3].append(math.pow(digit[i] - x, 4) * frq[i])

A1 = b[2][0]/n/S**3
A2 = b[3][0]/n/S**4 - 3

print("x = {}".format(x))
print("S = {}".format(S))
print("A1 = {}".format(A1))
print("A2 = {}".format(A2))
print("Вероятности: {}".format(p))

print("xi*ni = {}".format(b[0]))
print("xi^2*ni = {}".format(b[1]))
print("третий столбец = {}".format(b[2]))
print("четвертый столбец = {}".format(b[3]))