import sympy def np(a): return(sympy.nextprime(a)) def nnp(a,n): t=[] t.append(a) for i in range(n): a=np(a) t.append(a) return(t) def diffarr(t): a=len(t) s=[] for i in range(a-1): s.append(t[i+1]-t[i]) return(s) def sincarr(t): a=len(t) s=True for i in range(a-1): if t[i+1] <= t[i]: s=False return(s) def sdecarr(t): a=len(t) s=True for i in range(a-1): if t[i+1] >= t[i]: s=False return(s) n1=1 p=2 while p>0: n=n1 m1=nnp(p,n) m2=nnp(m1[-1],n) if sdecarr(diffarr(m1))==True and sincarr(diffarr(m2))==True: print(m1[-1]) n1=n1+1 p=np(p)