|
Текст подпрограммы и версий sf12r_c.zip sf12d_c.zip |
Тексты тестовых примеров tsf12r_c.zip tsf12d_c.zip |
Вычисление интеграла Доусона.
Функция sf12r_c производит вычисление интеграла
x
y = e-x*x ∫ et*tdt
0
для всех вещественных x, используя рациональные чебышевские аппроксимации.
W.Cody, K.A.Paciorek, H.C.Thacher, Chebyshev Approximations for Dawson's Integral, Math. Comp., 24, 1970.
real sf12r_c (real *x)
Параметры
| x - | заданное значение аргумента x (тип: вещественный). |
Версии
| sf12d_c - | вычисление интеграла Доусона с повышенной точностью. |
Вызываемые подпрограммы: нет
Замечания по использованию
| Функции sf12r_c и sf12d_c производят вычисления для всех x, представимых на машине с обычной и удвоенной точностью соответственно. |
int main(void)
{
/* Local variables */
extern float sf12r_c(float *);
static float x, y;
x = -1.6f;
y = (float)sf12r_c(&x);
printf("\n %16.7e %16.7e \n",x,y);
return 0;
} /* main */
Результат:
y = -0.399939894324