# Вторая пачка часть 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

Индекс
Элементарные функции Линейные уравнения Нелинейные уравнения Случайные числа
adobe acrobat скачать