Computing sums of floor() or frac() is related to non-trivial results
in number theory, e.g. `sum(n//i for i in range(1,n+1))` is
https://oeis.org/A006218
sage: var('x,y,n')
(x, y, n)
sage: ig=integrate(floor(n/x),x);ig
x*floor(n/x)
sage: n=13;a=1;b=n
sage: f=floor(n/x)
sage: i1=integrate(f,x);i1a=i1(x=b)-i1(x=a);i2=integrate(f,x,a,b);CC(i2-i1a)
Floor definite integration: can only handle linear < or > condition
28.3417388167000
sage: CC(i2),CC(i1a)
(28.3417388167000, 0.000000000000000)
sage: sum(n//i for i in range(1,n+1))
37
sage: i1
x*floor(13/x)
sage: pl1=plot(i1,a,b);pl1
--
You received this message because you are subscribed to the Google Groups
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-devel+unsubscr...@googlegroups.com.
To view this discussion visit
https://groups.google.com/d/msgid/sage-devel/CAGUWgD_fNcYjjaf%3DS5OAqvRd7wR2e7HbJzipZEhFUBB8KBuw4A%40mail.gmail.com.
