Вторая пачка часть 129
0777778;
C6 [1] := 0.0659722;
C6 [2] := 0.2604167;
C6 [3] := 0.1736111;
C6 [4] := 0.1736111;
C6 [5] := 0.2604167;
C6 [6] := 0.0659722;
C7 [1] := 0.0488095;
C7 [2] := 0.2571429;
C7 [3] := 0.0321429;
C7 [4] := 0.3238095;
C7 [5] := 0.0321429;
C7 [6] := 0.2571429;
C7 [7] := 0.0488095;
C8 [1] := 0.0434606;
C8 [2] := 0.2070023;
C8 [3] := 0.0765625;
C8 [4] := 0.1729745;
C8 [5] := 0.1729745;
C8 [6] := 0.0765625;
C8 [7] := 0.2070023;
C8 [8] := 0.0434606;
NA := N - K + 1;
IF (N - K) <= 0 THEN
BEGIN
F := 0.0;
FNK := F;
EXIT;
END;
NJ := KU (NA);
CASE NJ OF
1 : BEGIN
F := 0.0; FNK := F;
EXIT;
END;
2 : FOR I := 1 TO NJ
C [I] := C2 [I];
3 : FOR I := 1 TO NJ
C [I] := C3 [I];
4 : FOR I := 1 TO NJ
C [I] := C4 [I];
5 : FOR I := 1 TO NJ
C [I] := C5 [I];
6 : FOR I := 1 TO NJ
C [I] := C6 [I];
7 : FOR I := 1 TO NJ
C [I] := C7 [I];
8 : FOR I := 1 TO NJ
C [I] := C8 [I];
END; { **** CASE **** }
N1 := NJ - 1;
D := N1*H;
NZ := NA DIV N1;
IF (NJ - 2) < = 0 THEN NZ := NZ - 1;
F := 0.0;
FOR J := 1 TO NZ DO
BEGIN
JJJ := N1*(J - 1) + K - 1;
FOR I := 1 TO NJ DO
BEGIN
III := I + JJJ;
F := F + C[I]*D*P[III];
END;
END;
FNK := F;
END.
Подпрограмма FNK тестировалась на ЭВМ Pen¬ti¬um-166 для массивов от 21 до 5001 чле¬нов. Ре¬зультаты по¬ка¬за¬ли, что подпрограмма обес¬пе¬чи¬ва¬ет за¬данную точность вы¬числений при разумных зна¬ч嬬ни¬ях шага сетки h, кроме то-го, она эффективно ис¬поль¬¬зо¬ва¬лась в системе мо¬де¬ли¬ро¬ва¬ния динамики не¬ли¬ней¬ных волн в плазме, описываемых урав¬нениями типа ура⬬не-ния Кадомцева - Петвиашвили, пред¬ставленной в се¬рии ра¬бот В.Ю.Белашова [Белашов, 1989, 1991a, б, 1997].
§ 7. ВЫЧИСЛЕНИЕ ИНТЕГРАЛА ПО РОМБЕРГУ
При решении проблемы построения оп¬ти¬маль¬ной квад¬ратуры часто используют различные ме¬то¬ды уточ¬нения приближений к интегралу, по¬лу¬ченные при по¬мощи про¬стей¬ших ква¬дра¬тур¬ных фор¬мул
Индекс
Элементарные функции
Линейные уравнения
Нелинейные уравнения
Случайные числа
