Miller-Rabin error probabilities

mr_iterations.py

#!/usr/bin/python

from math import log, pow, sqrt
import sys

# Estimate Miller-Rabin error probability based on bitsize (k)
# and number of MR tests (t)

# https://www.math.dartmouth.edu//~carlp/PDF/paper88.pdf

k = int(sys.argv[1])
t = int(sys.argv[2])

if t * 9 > k:
    print("invalid k,t")

err = log(1 / (pow(k, 3.0/2) * pow(2,t) * pow(t, -.5) * pow(4, 2-sqrt(t*k)))) / log(2)

print("p(%d,%d) = %.02f" % (k, t, err))

You need to provide k and t as command line arguments