# 6.2.py
# 2016-04-09
# $Id: 6.2.py 1.2 2020/11/22 23:29:49 s Exp $

import datetime

# ex   : E(x)
# varx : var(x)

def keisen():
  print(60 * "-")

# input
n = 6
x = [i * j for i in range(1, n + 1) for j in range(1, n + 1)]

# calculating
m    = len(x)
sumx = sum(x)
ex   = sumx / m
sumv = sum((xi - ex) ** 2 for xi in x)
varx = sumv / m

# output
print(datetime.date.today())
print("6.2.py")
f1 = "{:>5}{:>5}{:>5}{:>10}{:>10}{:>10}{:>15}"
print(f1.format("k", "i", "j", "x=i*j", "E(x)", "x-E(x)", "(x-E(x))^2"))
keisen()
k = 0
for i in range(1,n+1):
  for j in range(1,n+1):
    f2 = "{:5d}{:5d}{:5d}{:10d}{:10.2f}{:10.2f}{:15.2f}"
    print(f2.format(k , i, j, x[k], ex, x[k]-ex, (x[k]-ex)**2))
    k = k + 1
keisen()
print("{:15}{:10d}{:20}{:15.2f}".format("sum", sumx, "", sumv))
print()
print("{:10}{:10d}".format("m= ", m))
print("{:10}{:10.2f}".format("E(x)= ", ex))
print("{:10}{:10.2f}".format("var(x)= ", varx))
# eof