All published worksheets from http://sagenb.org
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