Текст подпрограммы и версий ( Фортран )
sf40r.zip  sf40d.zip 
Тексты тестовых примеров ( Фортран )
tsf40r.zip  tsf40d.zip 
Текст подпрограммы и версий ( Си )
sf40r_c.zip  sf40d_c.zip 
Тексты тестовых примеров ( Си )
tsf40r_c.zip  tsf40d_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf40r_p.zip , sf40e_p.zip
Тексты тестовых примеров ( Паскаль )
tsf40r_p.zip , tsf40e_p.zip

Функция:  SF40R

Назначение

Вычисление неполной гамма - функции  Гx (p) вещественного аргумента и вещественного параметра.

Математическое описание

SF40R вычисляет неполную гамма - функцию Гx (p)

                   x
     Гx (p) = ∫ e - t t p - 1 dt    ,
                  0 

где  x  и  p  - принимают вещественные значения, причем x ≥ 0 и p > 0.

M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions New York, Dover Publications, INC., 1965.

Использование

    REAL  FUNCTION  SF40R (X, P, IERR) 

Параметры

X - заданное значение аргумента x (тип: вещественный);
P - заданное значение параметра p (тип: вещественный);
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
IERR=65 - когда значение аргумента x меньше 0; значение функции полагается равным 3.4E38;
IERR=66 - когда значение параметра P неполной гамма - функции не принадлежит полуинтервалу (0, 21.4] ; значение функции полагается равным 3.4E38.

Версии

SF40D - вычисление неполной гамма - функции Гx (p) с повышенной точностью.

Вызываемые подпрограммы

SF39R - подпрограмма вычисления гамма - функции распределения вероятностей I (x, p);
SF39D - подпрограмма вычисления гамма - функции распределения вероятностей I (x, p) с повышенной точностью;
SF23C - подпрограмма вычисления гамма - функции и ее натурального логарифма от комплексного аргумента;
SF23P - подпрограмма вычисления гамма - функции и ее натурального логарифма от комплексного аргумента с повышенной точностью;
      UTSF10 -
      UTSF12  
подпрограммы выдачи диагностических сообщений при работе функций SF23C, SF39R, SF40R;
      UTSF11 -
      UTSF13  
подпрограммы выдачи диагностических сообщений при работе функций SF23P, SF39D, SF40D

Замечания по использованию

  Для функции SF40D параметры X и P (0 < P ≤ 408.909) должны иметь тип DOUBLE PRECISION и при IERR ≠ 0 значение SF40D полагается равным 1.7D308.

Пример использования

         X = 1.
         P = 1.
         Y = SF40R (X, P, IERR)

Результаты:

       Y  =  0.6321205588074 , 
       IERR = 0