Текст подпрограммы и версий 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