>>> prop_int_ksigma(1000,500,8,1)
0.084
>>> prop_int_ksigma(1000,500,8,2)
0.118
>>> prop_int_ksigma(1000,500,8,3)
0.192