yt = np.array(known_primes)
n = np.random.rand(20) + 2
i = 0
while True:
    p = np.stack([nth_prime_explicit(n[i], K=50) for i in range(20)])
    err = yt - p
    n += 0.01 * err
    print(p[0], np.mean(err**2))
    i+=1