Я пытаюсь построить сумму двух распределений LogNormal, и SymPy просто зависает, не возвращая результат (это прекрасно работает с двумя распределениями Normal) - я подумал, что это может быть связано с «длиной» хвоста, который долженвычислить и попытаться охватить его диапазоном X - это не помогло
Пример кода:
%matplotlib inline
from sympy.stats import LogNormal, Normal, density
from sympy import Symbol, simplify, pprint
from sympy.plotting import plot
x = Symbol("x")
### Works ###
X1 = Normal('n1', 1, 1) # Mean 0, standard deviation 1
X2 = Normal('n2', 2, 1) # Mean 0, standard deviation 1
plot(density(X1+X2)(x))
### Doesn't work ###
X1 = LogNormal('x1', 0, 0.25) # Mean 0, standard deviation 1
X2 = LogNormal('x2', 0, 1) # Mean 0, standard deviation 1
plot(density(X1+X2)(x))