Processing math: 100%

First time here? Check out the FAQ!

Ask Your Question
0

Sage doesn't show limit of ratio of sums

asked 10 years ago

nespinoza gravatar image

Hi everyone,

I'm new to Sage and I'm trying to obtain the following limit: limNNi=1x5/2iNi=1x2i with xi=(Ni)/(N1). What I tried was:

i = var('i')
assume(i>0,'integer')
N = var('N')
assume(N>0,'integer')
x_i = (N-i)/(N-1)
numerator = sum(x_i^(5./2.),i,1,N)
denominator = sum(x_i^(2.),i,1,N)
limit(numerator/denominator,N=oo)

However, this just gives me:

6*limit((N - 1)*sum(((N - i)/(N - 1))^2.5, i, 1, N)/(2*N^2 - N), N, +Infinity)

From what I can see, Sage has no problems with the denominator, the general formula of which was obtained easily. However, apparently it has problems obtaining the general formula for the numerator due to the fact that the exponent in the sum is not an integer. Is there a way to cope with this? How can I obtain the limit?

Thanks in advance for all the help!

Preview: (hide)

1 Answer

Sort by » oldest newest most voted
0

answered 10 years ago

rws gravatar image

Sage uses Maxima for formal sums and so depends on that package. Maxima can compute

sage: sum(1/i^3,i,1,N)
gen_harmonic_number(3, N)

but unfortunately Sage cannot use this result symbolically. A ticket for this is waiting for review: http://trac.sagemath.org/ticket/16671

However, to compute your sum Maxima would need to know about generalized harmonic numbers with rational index. So this functionality must first be implemented in Maxima. Even then, this won't give you the limit which, I think, could be computable with some holonomic package. Something like this is missing in Sage, at the moment, as well.

Finally, I could well be wrong in the above, and there may be a special shortcut for implementation.

Preview: (hide)
link

Comments

I see...well, I'll have to wait then. Simulations are good enough anyways to actually see the limiting value!

nespinoza gravatar imagenespinoza ( 10 years ago )

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 10 years ago

Seen: 331 times

Last updated: Nov 03 '14