Текст подпрограммы и версий
qsf2r_c.zip  qsf2d_c.zip  		Тексты тестовых примеров
tqsf2r_c.zip  tqsf2d_c.zip 

Подпрограмма:  qsf2r_c

Назначение

Вычисление интегралов 

                  B
(1)             ∫  f (x) sin( ρ x + φ ) dx ,
                A
                  B
(2)             ∫  f (x) cos( ρ x + φ ) dx
                A

по формулам интерполяционного типа 5 - ой степени точности с автоматическим выбором шага.

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

qsf2r_c вычисляет интегралы (1), (2) или один из них с заданной абсолютной погрешностью E по формулам интерполяционного типа для осциллирующих функций с автоматическим выбором шага на конечных отрезках интегрирования.

О.В.Брушлинская, Л.Г.Васильева. Набор стандартных программ приближенного вычисления однократных интегралов с автоматическим выбором шага. Сб. "Численный анализ на ФОРТРАНе", вып.8, Изд-во МГУ, 1974.

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

    int qsf2r_c (real *rint1, real *rint2, real *a, real *b, R_fp f,
            real *ro, real *fi, real *e, integer *l1, integer *l2)

   Параметры
rint1 -
rint2   		вещественные переменные, содержащие значения интегралов (1) и (2) соответственно;
a, b - заданные нижний и верхний пределы интегрирования (тип: вещественный);
f - имя вещественной подпрограммы - функции вычисления подинтегральной функции f (x);
ro - заданное значение параметра  ρ (тип: вещественный);
fi - заданное значение параметра  φ (тип: вещественный);
e - заданная абсолютная погрешность вычисления интегралов (тип: вещественный);
l1, l2 - задают режим работы подпрограммы (тип: целый); при этом, если
 

l1 = 1 и l2 = 0, то вычисляется интеграл (1);

l1 = 0 и l2 = 1, то вычисляется интеграл (2);

l1 = 1 и l2 = 1, то вычисляются интегралы (1) и (2); 
   Версии
qsf2d_c - вычисление с удвоенной точностью интегралов (1) и (2) по формулам интерполяционного типа 5 - ой степени точности с автоматическим выбором шага. 
   Вызываемые подпрограммы
utqs10_c - подпрограмма выдачи диагностических сообщений при работе подпрограммы qsf2r_c;
utqs12_c - подпрограмма выдачи диагностических сообщений при работе подпрограммы qsf2d_c. 
   Замечания по использованию
 

Подпрограмма не гарантирует задаваемую погрешность, хотя в большинстве случаев интегралы будут вычислены с требуемой точностью.

В подпрограмме qsf2d_c параметры rint1, rint2, a, b, f, ro, fi, e имеют тип double.

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

int main(void)
{
    /* Local variables */
    extern int qsf2r_c(float *, float *, float *, float *, R_fp,
                       float *, float *, float *, int *, int *);
    static float a, b, e;
    extern float f_c();
    static int l1, l2;
    static float r0, fi, int1, int2;

    a = 0.f;
    b = 1.f;
    r0 = 10.f;
    fi = 0.f;
    e = -1e-7f;
    l1 = 1;
    l2 = 1;
    qsf2r_c(&int1, &int2, &a, &b, (R_fp)f_c, &r0, &fi, &e, &l1, &l2);

    printf("\n %16.7e %16.7e \n",int1,int2);
    return 0;
} /* main */

float f_c(float *x)
{
    /* System generated locals */
    float ret_val, r__1;

/* Computing 4th power */
    r__1 = *x, r__1 *= r__1;
    ret_val = r__1 * r__1;
    return ret_val;
} /* f_c */
      

Результаты:

      int1  =  0.05382478
      int2  =  -0.079553512