Um dos problemas do código é não estares a alterar (incrementar) o valor de `j`, por isso acabas sempre com os mesmos valores. Uma implementação pode ser feita assim: def e_primo(num): # verificar se e primo return all(num%i!=0 for i in range(2,num)) qtd_tuples = 10 # input results = [] # armazenar resultados curr_prime = 1 while len(results) < qtd_tuples: # enquanto a quantidade de tuplas for inferior ao nosso input curr_prime += 2 if(e_primo(curr_prime) and e_primo(curr_prime+2)): results.append((curr_prime, curr_prime+2)) print(results) Output: [(3, 5), (5, 7), (11, 13), (17, 19), (29, 31), (41, 43), (59, 61), (71, 73), (101, 103), (107, 109)] [DEMONSTRAÇÃO][1] [1]: https://ideone.com/MSdeVF