Heatmap of the prime numbers. And of the prime numbers except the least significant digit http://liveatthewitchtrials.blogspot.ie/2017/10/yet-more-prime-pictures.html

PrimeHeatmap.py

#!/usr/bin/python3
import numpy as np


#the number of digits to use
numDigits = 60000000

y=450;

a = np.zeros( (10,10) )
previous=0


currentPrime=0
    

def getDigits():
	#0.12345678910111213…
	current = 1;
	print(current)
	while True:
		current_str = str(current);
		for digit in current_str:
			yield int(digit)
		current = current + 1


def getPrimeDigits():
	#0.12345678910111213…
	current = 1;
	global currentPrime
	newP =next_prime(current)
#	print(newP)
	while True:
		newP =next_prime(current)
		#print(newP)
		current_str = str(newP);
		current_str = current_str[:-1]#strip last digit to see if lack of 2,4,6,8,0 causes all issue
		currentPrime =newP
		for digit in current_str:
			yield int(digit)
		current = current + 1
        
gen = getPrimeDigits()


for i in range(numDigits):
	previous = digit-1
	digit = next(gen)
	a[digit-1,previous]+=1

	i=i+1    

print(currentPrime)
a


import matplotlib.pyplot as plt
import numpy as np


plt.imshow(a, cmap='bwr', interpolation='nearest')
plt.show()