Create Codedef dataforhist(n):     a = []     for i in range(2,1000001):         if is_prime(i):             a.append(int(Mod(i,n)))     return a  b = dataforhist(15) stats.TimeSeries(b).plot_histogram(normalize=False,aspect_ratio=0.001)