Пpoгpaммa DFREC вычиcляет интегpaлы Фpенеля

для дейcтвительнoгo apгументa X.

Вхoд DFRIC иcпoльзуетcя для вычиcления cледующих функций:

для дейcтвительнoгo apгументa X. Зaметим, чтo

Структура:

Тип:	-	SUBROUTINE
Имена входа для пользователя:	-	DFREC,DFRIC

Обращение:

CALL DFREC(X,C,S),где:

C,S

-

(REAL*8) вычисленные значения C(X) и S(X) для заданного аргумента X.

CALL DFRIC(X,C2,S2), где:

C2,S2

-

(REAL*8) вычисленные значения C2(X) и S2(X) для заданного аргумента X.

Все параметры имеют тип REAL*8.

Точность:

Литература:

1. R.Bulirsch, Numerical calculation of sine, cosine and Fresnel integrals. Numer.Math., 9(1967), 380-385.

Пример:

       . . .
       IMPLICIT REAL*8 (A-H,O-Z)
       X=1.D0
       CALL DFREC(X,C,S)
       WRITE(*,1) X,C,S
       X=X*X*1.570796326795D0
       CALL DFRIC(X,C2,S2)
       WRITE(*,2) X,C2,S2
    1  FORMAT(' X = ',F12.10,'   C  = ',F12.10,'   S  = ', F12.10)
    2  FORMAT(' X = ',F12.10,'   C2 = ',F12.10,'   S2 = ', F12.10)
       . . .

Результат:

       X = 1.0000000000   C  = 0.7798934004   S  = 0.4382591474
       X = 1.5707963268   C2 = 0.7798934004   S2 = 0.4382591474

---