Yellowstone sequence: Difference between revisions

<lang J>

GCD=: +.

relatively_prime=: 1 = GCD

yellowstone=: monad define

start=. #\ i. 4 + y NB. prepare minimal starting values

s=. 3 {. start NB. the sequence vector

start=. 3 }. start

while. y > # s do.

z=. {. start NB. z is the lowest number not in the sequence

while.do.

if. 0 1 -: (_2 {. s) relatively_prime z do.

if. z -.@e. s do.

break.

end.

end.

z =. >: z

end.

start=. start -. z NB. remove z from the list of starting values

s=. s , z

end.

s

)

</lang>

<pre>

yellowstone 30

1 2 3 4 9 8 15 14 5 6 25 12 35 16 7 10 21 20 27 22 39 11 13 33 26 45 28 51 32 17

load'plot'

'marker'plot yellowstone 100

</pre>