Вторая пачка часть 217
.., bМ; IER (тип in¬¬te¬ger) - признак ошиᬬки во входных параметрах. Если IER = 0 - нет ошибки; IER = 1 при N < М; IER = 2 при М<0.
PROCEDURE FORINT (FNT:MAS21;N,M:INTEGER;
VAR A,B:MAS; VAR IER : INTEGER);
LABEL 70,100;
VAR I, AN, J : INTEGER;
Q,CONS,COEF,S1,S,C1,C,FNTZ,U0,U1,U2:REAL;
BEGIN
IER := 0;
IF M<0 THEN
BEGIN
IER := 2;
EXIT;
END;
IF M-N>0 THEN
BEGIN
IER := 1;
EXIT;
END;
AN := N;
COEF := 2.0 / (2.0*AN+1.0);
CONS := 3.14159265*COEF;
S1 := SIN (CONS);
C1 := COS (CONS);
C:=1.0;
S :=0.0;
J := 1;
FNTZ := FNT[1];
70:
U2 := 0.0;
U1 := 0.0;
I := 2*N + 1;
REPEAT
U0 := FNT [I] + 2.0 *C*U1 - U2;
U2 := U1;
U1 := U0;
DEC (I);
UNTIL I<=1;
A[J] := COEF*(FNTZ+C*U1-U2);
B[J] := COEF*S*U1;
IF J>=(M+1) THEN GOTO 100;
Q := C1*C - S1*S;
S := C1*S + S1*C;
C := Q; INC (J);
GOTO 70;
100:
A[1] := A[1]/2.0;
END.
Для проверки работы процедуры и ее тес¬ти¬ро¬ва¬ния приведены резуль¬таты (табл. 7.6), полученные на отрезке [0, 2???при разложении функ¬ции
f(х) = sin3(х) - 2cos(х) - 3
в ряд Фурье. Для расчетов М = 6; N = 10; 2N + 1 = 21.
Таблица 7.6
Номер Коэффициент Фурье
п/п А i В i
1 -3.0000 0.0000
2 -2.0000 0.75000
3 0.00000 0.00000
4 0.00000 -0.2500
5 0.00000 0.00000
6 0.00000 0.00000
§ 3. НЕПОЛНЫЕ ЭЛЛИПТИЧЕСКИЕ ИНТЕГРАЛЫ I И II РОДА
Неполный эллиптический интеграл I рода в нор¬маль¬ной форме Лежандра определяется выражением [Корн, Корн, 1984; Янке и др., 1968]
,
(k - модуль), которое мо¬жет быть аппроксимировано мнгочленом
(5.20)
где n = 1 - k2, |?(k)| ? 2?10-8; ai, bi - значения ко¬эф-фи¬ци¬ен¬тов, которые при i = 0, 1, 2, 3, 4 приведены в табл. 5.5.
Таблица 5.5
№ Значения коэффициентов
i ai bi
0 1.38629436112 0.5
1 0.09666344259 0.12498593597
2 0.03590092383 0.06880248576
3 0.03742563713 0.03328355346
4 0.01451196212 0.00441787012
Процедура-функция ELLIPT1 вычисления K(k) по¬стре¬на на основании разложения (5
Индекс
Элементарные функции
Линейные уравнения
Нелинейные уравнения
Случайные числа
