# Вторая пачка часть 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. ВЫЧИСЛЕНИЕ ИНТЕГРАЛА ПО РОМБЕРГУ При решении проблемы построения оп¬ти¬маль¬ной квад¬ратуры часто используют различные ме¬то¬ды уточ¬нения приближений к интегралу, по¬лу¬ченные при по¬мощи про¬стей¬ших ква¬дра¬тур¬ных фор¬мул

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