Как стать автором
Обновить

Комментарии 8

подходящий пакет с zeta я нашел только в mpmath

Риманова Зета функция, не оно?

from scipy.special import zeta

Я знаю, но судя по описанию, это сумма ряда, а не аналитическое расширение.

>>>zeta(3)
1.2020569031595942
>>> zeta(complex(0.5,14.1))
Traceback (most recent call last):
File "", line 1, in
File "C:\...\site-packages\scipy\special_basic.py", line 2572, in zeta
return ufuncs._riemann_zeta(x, out)
TypeError: ufunc '_riemann_zeta' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''

Существует универсальное интегральное представление для дзета-функции, оно задается формулой Абеля-Плана
image
и работает для всех s, не равных 1

С его помощью можно построить график для вещественных значений аргумента
График дзета-функции для вещественных значений аргумента
image

Было бы здорово посмотреть исходники расчетов и генерации изображений

Я вам пришлю чуть позже

а каким тулом генерятся такие красивые картинки?

Я бы попробовал X[n+1]=Zeta(X[n]+C)**2-Zeta(0)**2. Задумка в том, чтобы аналогично обычному Мандельброту расквадрачивать вокруг нуля, вот только нулей бесконечно много с нерегулярным шаблоном. Если в C подставить ноль Зета-функции, то X[0]=0, X[1]=0, X[2]=0 и так далее, а если что-то около, то должно расфракталиться, как в Мандельброте, вот только вокруг каждого нуля должны быть следы и всех остальных нулей.

Зарегистрируйтесь на Хабре, чтобы оставить комментарий

Публикации