OFFSET
1,3
COMMENTS
MAPLE
SIGSEQsort := proc(l1::list, l2::list)
if op(1, l1) < op(1, l2) then
return true ;
else
return false ;
end if
end proc:
SIGSEQ := proc(x, vmax)
local TBsrtd, i, j ;
TBsrtd := [] ;
for i from 1 do
if i > vmax then
break ;
end if;
for j from 1 do
if evalf(i+j*x) > vmax then
break ;
end if;
TBsrtd := [op(TBsrtd), [evalf(i+j*x), i]] ;
end do:
end do:
sort(TBsrtd, SIGSEQsort) ;
[seq(op(2, v), v=%)] ;
end proc:
Digits := 100 ;
SIGSEQ(log(2), 50.0) ;
CROSSREFS
KEYWORD
nonn
AUTHOR
R. J. Mathar, May 28 2024
STATUS
approved