(phixonline)-->
constant lim = 1000
sequence van_eck = repeat(0,lim),
pos_before = repeat(0,lim)
for n=1 to lim-1 do
integer vn = van_eck[n]+1,
prev = pos_before[vn]
if prev!=0 then
van_eck[n+1] = n - prev
end if
pos_before[vn] = n
end for
printf(1,"The first ten terms of the Van Eck sequence are:%v\n",{van_eck[1..10]})
printf(1,"Terms 991 to 1000 of the sequence are:%V\n",{van_eck[991..1000]})