ADS_coding.ipynb

Kernel: SageMath 10.6

data=list(map(lambda k:(k,binomial(1024,k)* 0.002^k * 0.998^(1024-k)),[0,1,2,3,4,5]))

data

[(0, 0.128728372417859),
 (1, 0.264164034781337),
 (2, 0.270781370321952),
 (3, 0.184862097841707),
 (4, 0.0945612233949811),
 (5, 0.0386582957366255)]

data2=list(map(lambda k:binomial(1024,k)* 0.002^k * 0.998^(1024-k),[0,1,2,3,4,5,6]))

list_plot(data,plotjoined=True)

bar_chart(data2)

binomial(1024,7)* 0.002^7 * 0.998^(1024-7)

0.00383454893526175

bar=bar_chart(data2)

bar.save('fig-bsc-example.png')