OFFSET
1,1
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
Primes prime(n) such that (-prime(n) + 2*prime(n+1) - prime(n+2))/((1 - prime(n) + prime(n+1))^(3/2)) < 0.
MATHEMATICA
d2[n_] = Prime[n + 2] - 2*Prime[n + 1] + Prime[n]; d1[n_] = Prime[n + 1] - Prime[n]; k[n_] = -d2[n]/(1 + d1[n])^(3/2); Flatten[Table[If[k[n]<= 0, Prime[n], {}], {n, 1, 100}]]
PROG
(Haskell)
import Data.List (findIndices)
a147813 n = a147813_list !! (n-1)
a147813_list = map (a000040 . (+ 1)) $ findIndices (>= 0) a036263_list
-- Reinhard Zumkeller, Jan 20 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Nov 13 2008
EXTENSIONS
Edited by Alonso del Arte and Joerg Arndt, Nov 01 2013
STATUS
approved