Georgi Guninski on Fri, 20 Jan 2023 08:21:42 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Computing integer sequences via integrals


I am trying to compute integer sequences via closed
forms of definite integrals, but probably this is impossible.

The question is: which integer sequences can be computed
by integrals?

Partial results:

Trying to "discretize" the reals, we are using floor()
and the area of the definite integrals.

For a start consider the plot of floor(x).
? ploth(x=1,10,floor(x))

The area is the union of rectangles and for integer
bounds the area is integer.

sagemath computes the following indefinite integrals:
sage: integrate(floor(x),x)
1/2*(2*x - floor(x) - 1)*floor(x)
sage: integrate(floor(x)^2,x)
x*floor(x)^2

The area for floor(x) and the sum(i=1,N-1,x) appear to be computed correctly:
sage: I1=integrate(floor(x),x)
sage: N=13;[I1(x=N)-I1(x=1)-N*(N-1)/2]
[0]

The area for floor(x)^2 doesn't work for me, why?

floor() might not be elementary function, but there is closed
form for it using exp() and log(), taking the principal
branch of the logarithm:

? floor1(x)=x - 1/2*I*log(-exp(-2*I*Pi*x))/Pi - 1/2
%12 = (x)->x-1/2*I*log(-exp(-2*I*Pi*x))/Pi-1/2
? floor1(13.1)
%13 = 13.000000000000000000000000000000000000 + 0.E-39*I