2036.sagews

kl = (0.3, 0.2, 0.2, 0.13, 0.09, 0.04, 0.025, 0.01, 0.004, 0.001)

M = 0
for k in range(10):
    M += kl[k]*(1-0.5^k)
print M

JM = 0
for k in range(10):
    tmp = (1-2*(0.5^k))
    if k == 0:
        tmp += 1
    JM += kl[k]*tmp
#    print kl[k]*tmp
print JM

print JM/M

0.526388671875000 0.352777343750000 0.670184148327897