Я пытаюсь получить асимптотическое расширение c для кумулятивного нормального распределения от sympy. Я пытался
from sympy import *
x = Symbol('x')
N = erf(2**(1/2)*x/2)/2 + 1/2 #Normal(0, 1)
print(N.series(x,oo,3))
Но я получаю:
sympy.core.function.PoleError: Asymptoti c расширение erf вокруг [oo] не реализовано.
Это сообщение об ошибке довольно ясно. Мне интересно, хотя, есть ли способ обойти это, например, выражая cumnorm без erf или подобного.
Спасибо!