Уведомления

Группа в Telegram: @pythonsu

#1 Дек. 18, 2018 15:37:38

ADEIIT
Зарегистрирован: 2018-10-09
Сообщения: 4
Репутация: +  0  -
Профиль   Отправить e-mail  

Расчет дзета функции Римана

нужно посчитать функцию фримана
e можно принять за 2.7 решение можно закодить с помощью mpmath, но я не в курсе как поставить свои аргументы

Прикреплённый файлы:
attachment формула.jpg (70,7 KБ)

Офлайн

#2 Дек. 18, 2018 15:43:56

JOHN_16
От: Россия, Петропавловск-Камчатск
Зарегистрирован: 2010-03-22
Сообщения: 3292
Репутация: +  221  -
Профиль   Отправить e-mail  

Расчет дзета функции Римана

а “к” чему равно?



_________________________________________________________________________________
полезный блог о python john16blog.blogspot.com

Офлайн

#3 Дек. 18, 2018 15:45:43

ADEIIT
Зарегистрирован: 2018-10-09
Сообщения: 4
Репутация: +  0  -
Профиль   Отправить e-mail  

Расчет дзета функции Римана

k берутся от 1 до бесконечности, это сумма должна светись к значению определеному

Офлайн

#4 Дек. 18, 2018 16:59:59

JOHN_16
От: Россия, Петропавловск-Камчатск
Зарегистрирован: 2010-03-22
Сообщения: 3292
Репутация: +  221  -
Профиль   Отправить e-mail  

Расчет дзета функции Римана

а, ну если

real and complex floating-point arithmetic with arbitrary precision
то примерно посчитать можно



_________________________________________________________________________________
полезный блог о python john16blog.blogspot.com

Офлайн

#5 Дек. 19, 2018 06:56:41

scidam
Зарегистрирован: 2016-06-15
Сообщения: 288
Репутация: +  35  -
Профиль   Отправить e-mail  

Расчет дзета функции Римана

Фактически, вам необходимо вычислить значение полилогарифма.
При определенных значениях параметра `s` этот ряд можно просуммировать аналитически и использовать соответствующую формулу.
В общем случае нужно задать какой-нибудь параметр, например, abstol=10e-8, и
считать по рекуррентной формуле, S(n+1) = S(n) + eta**n*theta**(n+1)/(n+1)**s, если
добавляемое значение будет меньше abstol, то прекращать итерации и считать полученное значение S(n) результатом (приближенным, конечно).
Код может быть приблизительно таким:

 def get_multilog(s, eta, theta, abstol=1.0e-8, nmax=1.0e+5):
    result = 0.0
    n = 0
    tol  = abstol
    while (tol > abstol) and (n < nmax):
        tol = eta ** n * theta ** (n + 1) / (n + 1) ** s
        result  += tol
        n += 1






Офлайн

Board footer

Модераторировать

Powered by DjangoBB

Lo-Fi Version