Prime spiral from matplotlib import pyplot as plt import numpy as np from sympy import isprime upper_limit = pow(2, 17) theta = [] r = [] ptheta = [] pr = [] def f(x): return (np.sqrt(x),x) for i in range(1, upper_limit): res = f(i) r.appe…