Блять, пик перепутал
спасайте, аноны
bump
bump
Вот те гамма, а функции я не знаю.
нахуя?
Есть же готовая реализация
https://www.gnu.org/software/gsl/manual/html_node/Gamma-Functions.html
да лаба ебучая
спасибо, ананас
попробую разобраться
Первая лаба у уже не школьника?
Берешь короче Стирлинга или Ланцоша и аппроксимируешь. Для больших аргументо Спаудж тоже норм.
Реализация есть даже в вики.
Так-то имея готовую формулу для аппроксимации реализовать это — ничего сложного.
Если хочешь понять как получить эту формулу придется разбираться в матане. Например: https://web.viu.ca/pughg/phdThesis/phdThesis.pdf
Чем ты искал?
В педивикии даже реализация на питоне есть
def gamma(z): # great function from Wiki, but maybe could use memorization?
epsilon = 0.0000001
def withinepsilon(x):
return abs(x - abs(x)) <= epsilon
from cmath import sin,sqrt,pi,exp
p = [ 676.5203681218851, -1259.1392167224028, 771.32342877765313,
-176.61502916214059, 12.507343278686905, -0.13857109526572012,
9.9843695780195716e-6, 1.5056327351493116e-7]
z = complex(z)
# Reflection formula
if z.real < 0.5:
result = pi / (sin(piz) gamma(1-z))
else:
z -= 1
x = 0.99999999999980993
for (i, pval) in enumerate(p):
x += pval/(z+i+1)
t = z + len(p) - 0.5
result = sqrt(2pi) t(z+0.5) exp(-t) x
if withinepsilon(result.imag):
return result.real
return result
Анон, я заебался.
Перерыл много интернета, ничего рабочего и дельного не смог найти. Поясни мне, полуёбку, как реализовать расчет этой самой Г(х) на ЯП. Псевлокод, алгоритм - пойдёт всё. Я понять хочу, господа, понять, как считать это.
Да снизойдёт и благословит каждого не прошедшего мимо анона Солнцеликий Abu