from math import * from kandinsky import * fill_rect(0, 0, 320, 222, (0, 0, 0)) def is_prime(n): i = 3 for i in range(i, int(sqrt(n)) + 1): if n % i == 0: return False return True def goldbach(n): if n == 4: return 1 else: sol = 0 for i in range(3, n // 2 + 1, 2): if is_prime(i) and is_prime(n - i): sol += 1 return sol for i in range(4000): x = 4 + 2 * i y = goldbach(x) set_pixel(int(50 + x / 30), 160 - y, (255, )*3)