Предел частичных сумм - PullRequest
Новая
       64

Предел частичных сумм

0 голосов
/ 20 декабря 2018

Итак, я хотел бы рассчитать предел частичных сумм, где a_n равно sp.Lambda(x, 1/x * sp.factorial2(2*x-1)/sp.factorial2(2*x)).Я попробовал это: 

sp.limit(sp.summation(a(i), (i,1,n)), n, sp.oo)

Но когда я запускаю его, он начинает вычисляться, но никогда не заканчивается, поэтому я не получаю ответ.Интересен тот факт, что если я попытаюсь вычислить бесконечную сумму, я получу ответ: 

sp.summation(a(i), (i,1,sp.oo))
>>>2log(2)

Я даже могу создать последовательность частичных сумм: 

s = sp.SeqFormula(sp.summation(a(i), (i,1,n)), (n, 1, sp.oo))

Что, дляпример в Maple, оценивает ограничение своего рода.Так есть ли способ рассчитать такой предел в симпы?На самом деле я хотел бы получить ответ, аналогичный суммированию: 2log(2), но с использованием предела частичных сумм.

1 Ответ

0 голосов
/ 20 декабря 2018 
def limsum(it, epsilon):
    ret = 0
    x = next(it)
    while x<epsilon:
        ret = ret + x
        x = next(it)
    return ret
...