Обратная неполная гамма функция
x = gammaincinv(y,a)
x = gammaincinv(y,a,tail)
x = gammaincinv(y,a) выполняет обратную неполную гамма функцию для соответствующих элементов y и a, такого that y = gammainc(x,a). Элементы y должны быть в закрытом интервале [0,1], и те из a должны быть неотрицательными. y и a должны быть действительными и тот же размер (или любой может быть скаляром).
x = gammaincinv(y,a,tail) задает хвост неполной гамма функции. Выбором является 'lower' (значение по умолчанию), чтобы использовать интеграл от 0 до x или 'upper', чтобы использовать интеграл от x до бесконечности.
Эти два варианта связаны как:
gammaincinv(y,a,'upper') = gammaincinv(1-y,a,'lower').
Когда y близко к 0, опция 'upper' обеспечивает способ вычислить x более точно, чем путем вычитания y от 1.
[1] Cody, J., обзор разработки программного обеспечения для специальных функций, примечаний лекции в математике, 506, числовой анализ Данди, Г. А. Уотсон (редактор)., Springer Verlag, Берлин, 1976.
[2] Abramowitz, M. и И.А. Стегун, Руководство Математических функций, Национального бюро стандартов, Прикладной Математики. Серия № 55, Дуврские Публикации, 1965, секунда. 6.5.