Простейшие квадратурные формулы Формулы прямоугольников.Как известно, определенный интеграл в силу своего построения есть предел интегральных сумм: , (22) каждая, из которых соответствует некоторому разбиению : отрезка и произвольному набору точек для каждого разбиения; . Ограничиваясь конечным числом слагаемых в правой части равенства (22) и принимая в качестве набора те или иные значения аргумента из отрезков , можно получить различные формы приближенного интегрирования. Так, принимая в качестве набора значения левых или правых концов отрезков , получим соответственно формулу левых или правых прямоугольников : , (23) . (24) Пример 1: С помощью формул левых и правых прямоугольников вычислить , полагая n = 4. Зная приделы интегрирования a = 1 и b = 9, находим шаг ; тогда точками разбиения служат , , , , , а значения подынтегральной функции в этих точках таковы: Далее найдем числовое значение интеграла, пользуясь формулой (23):  Если вычисление определенного интеграла произвести по формуле (24), то получим:  Наиболее часто используемой формулой, основанной на идее представления определенного интеграла в виде интегральной суммы, является формула прямоугольников, где в качестве берут середины отрезков . Для равномерной сетки эта формула имеет следующий вид: , (25) где ; ; . Выражение для остаточного члена квадратурной формулы (25): (26) Таким образом, оценку погрешности квадратурной формулы (25) можно представить в следующем виде: (27) где  Суммарная вычислительная погрешность составит (28) Эта погрешность не зависит от числа разбиений отрезка интегрирования, а пропорциональна только его длине. Пример 2: Вычислить с помощью формулы прямоугольников интеграл , полагая n = 4. Оценить погрешность полученного приближенного значения. По заданным пределам интегрирования и числу разбиений n определим шаг: . далее на основании формулы (25) имеем  Вычислив необходимые значения функции с тремя верными в узком смысле знаками , получим  Погрешность метода оценим по формуле (27), для чего предварительно найдем максимум абсолютной величины второй производной подынтегральной функции:  Таким образом, погрешность метода есть  Пользуясь формулой (28), найдем вычислительную погрешность:  Следовательно за полную погрешность приближенного вычисления интеграла можно принять , а окончательный ответ записать в виде Для сравнения приведем несколько знаков точного значения вычислительного интеграла: . Формула трапеций.Перейдем теперь к другому способу построения квадратурных формул, связанному с аппроксимацией подынтегральной функции интерполяционным многочленом. Рассмотрим простейший случай аппроксимации многочленом первой степени с узлами в точках a и b: ; . Интегрируя правую и левую части этого равенства и используя вторую теорему о среднем значении функции при интегрировании последнего слагаемого правой части, находим: ; . Таким образом, предполагая, что отрезок интегрирования мал, получаем квадратурную формулу, называемую формулой трапеций: (29) с остаточным членом ; . (30) Используя выражение (27) для остаточного члена, опенку погрешности квадратурной формулы (29) можно представить в виде , (31) где . Полученные выражения для остаточного члена (30) и погрешности (31) показывают, что квадратурная формула (29) является точной для всех линейных функций, поскольку вторая производная таких функций равна нулю, а, следовательно, равны пулю остаточный член и погрешность. Пример 3: Вычислить с помощью формулы трапеций интеграл . Оценить погрешность полученною приближенного значения. На основании формулы (29) имеем . Вычислив необходимые значения функции, получим . Погрешность метода оценим по формуле (31), используя значение М=2: . Вычислительная погрешность, очевидно, равна нулю, так как значения функции и , найдены абсолютно точно. Итак, окончательно имеем . Отметим, что в примере 3 получилось гораздо менее точное решение, чем в примере 2. Однако использование в примере 3 формулы трапеций имеет свои преимущества. Во-первых, если подынтегральная функция задана в виде таблицы ее значений в узлах , то для использования формулы прямоугольников необходимо определить значения этой функции еще и в точках , что вносит дополнительные трудности и дополнительную погрешность. Во-вторых, в примере 3 значения подынтегральной функции были вычислены всего лишь в двух точках, в то время как в примере 2 - в четырех точках, что, естественно, потребовало большего времени. Ценность квадратурной формулы определяется не только формой ее остаточного члена (или погрешности), но и другими факторами, например временем счета. 3. Объекты исследования, оборудование, материалы и наглядные пособия Объектом исследования данной лабораторной работы являются система линейных алгебраических уравнений и методы ее численного решения. Для выполнения работы необходимы ПК и соответствующее программное обеспечение: - MS WINDOWS; - Delphi (любая другая оболочка языка высокого уровня); - MS Office (для оформления отчета). 4. Задание на работу (рабочее задание) Разработать программное обеспечение для численного расчета: 1) Первой и второй производных функции, заданной таблично; 2) с помощью формул левых, правых, центральных прямоугольников и формулы трапеции определенного интеграла. Варианты: № варианта | Задание | № варианта | Задание | № варианта | Задание | 1. |  | 2. |  | 3. |  | 4. |  | 5. |  | 6. |  | 7. |  | 8. |  | 9. |  | 10. |  | 11. |  | 12. |  | 13. |  | 14. |  | 15. |  | 16. |  | 17. |  | 18. |  | 19. |  | 20. |  | 21. |  | 22. |  | 23. |  | 24. |  | 25. |  | 26. |  | 27. |  | 28. |  | 29. |  | 30. |  | |