Я пытаюсь закодировать это гамма-распределение на Python; однако мне трудно понять, как включить интеграл в функцию.
Распространение:
Будет ли первая половина уравнения перед интегралом определяться как отдельная функция, а затем умножаться к функции, которая определяется интегралом?
Для более подробного контекста: «t» — это моя переменная времени, которая состоит из моих точек x-данных в диапазоне от 0 до 24.< /p>
Это моя попытка написать это:
Код: Выделить всё
#Define Gamma Function
def func(t,g,n):
y = (1/gamma(g)*n*g)
return y
def func1(t,g,s,n,d):
y1 = (s**(g-1))*(np.exp(1)**(-s/n))*d*s
return y1
I = integrate.quad(func1,0,t)
funcdata = func(xdata,2,1)*func1(xdata,1,1,1,1)
plt.plot(xdata,funcdata,label="Model")
Источник: https://stackoverflow.com/questions/689 ... a-function
