Вычисление интеграла фукции f(x) методом Симпсона

Вычисление интеграла фукции f(x) методом Симпсона

С О Д Е Р Ж А Н И Е

Введение . 2

1. Постановка задачи 3

2. Математическая часть . 4

3. Описание метода решения задачи 9

4. Описание алгоритма решения задачи . 10

5. Текст программы . 11

6. Результаты работы программы 15

Заключение 16

Список использованных источников: 17

Введение

История появления и развития персональных компьютеров является одним из наиболее впечатляющих явлений нашего века. С момента появления первых образцов персональных компьютеров прошло меньше 25 лет, но сейчас без них уже немыслимо огромное количество областей человеческой деятельности - экономика, управление, наука, инженерное дело, издательское дело, образование, культура и т.д. Интерес к персональным компьютерам постоянно растет, а круг их пользователей непрерывно расширяется. В число пользователей ПЭВМ вовлекаются как новички в компьютерном деле, так и специалисты по другим классам ЭВМ.

Язык Паскаль - это один из наиболее распространённых языков программирования 80-90х годов , поддерживающий самые современные методологии проектирования программ (нисходящее, модульное проектирование, структурное программирование) имеют свою достаточно богатую историю развития.

Новую жизнь языку дала фирма Борланд, разработавшая на его базе семейство Паскаль – систем, называемых Турбо Паскалем. Интегрированная среда, обеспечивающая многооконную разработку программной системы, обширный набор встроенный в неё средств компиляции и отладки , доступный для работы через легко осваиваемое меню, - всё это обеспечивает высокую производительность труда программиста, недостижимую при работе со старыми средами.

Язык Турбо Паскаль хорошо подходит для обучения программированию.

1. Постановка задачи

Заданием на курсовую работу является создание программы на языке программирования Турбо Паскаль, которая должна осуществлять решение следующей задачи :

Вычислить приближённое значение интеграла функции f(x) на интервале с точностью до 0.01 методами Симпсона и трапеции с целью сравнения.

Интегрируемая функция: .

Определить метод, который решает поставленную задачу за минимальное число повторений.

Построить график функции f(x) на заданном интервале. Решить поставленную задачу с использованием функций и процедур алгоритмического языка Турбо Паскаль.

2. Математическая часть

Для приближённого вычисления интеграла функции f(x) используются методы приближённого интегрирования, наиболее употребительные из них основаны на замене интеграла конечной суммой. Для вычисления промежуток от a(x0) до b(xn) разбивается на n равных частей, и для точек деления x0 , x1 , x2 , x3 , . . . , xn-1 , xn вычисляются значения интегрируемой функции y. Затем необходимо воспользоваться формулой приближённого интегрирования:

1) Формула трапеций (рис.1) :

.(1)

Рис.1.

2) Формула Cимпсона (парабол) (рис.2) :

(2)

Рис.2.

В моей курсовой работе рассматривается приближенное вычисление интеграла (1)

При его аппроксимации заменим функцию f(x) параболой, проходящей через точки т.е представим приближенно f(x) в виде

где - интерполяционный многочлен Лагранжа второй степени,

. (2)

Проводя интегрирование получим

Таким образом приходим к приближенному равенству

(3)

Котрое называется формулой Симпсона или формулой парабол.

На всем отрезке [a,b] формула Симпсона имеет вид

Чтобы не использовать дробных индексов можно обозначить

xi=a+0,5hi, fi=f(xi), i=1,2,…,2N, hN=b-a

и записать формулу Симпсона в виде

(4)

Прежде чем переходить к оценке погрешности формулы (3) заметим, что она является точной для любого многочлена третьей степени, т.е. имеет место точное равенство

если f(x)=a0+a1x+a2x2+a3x3. Это утверждение нетрудно проверить непосредственно.

Для оценки погрешности формулы Симпсона воспользуемся интерполяционным многочленом Эрмита. Построим многочлен третьей степени H3(x) такой, что

.

Такой многочлен существует и единствен.

Однако нам даже не потребуется явный вид многочлена H3(x). Вспоминая, что формула Симпсона точна для любого многочлена третьей степени, получим

(5)

Представим теперь f(x) в виде

f(x)=H3(x)+ri(x), xÎ[xi-1,xi], (6)

где ri(x) – погрешность интерполирования многочленом Эрмита H3(x). Интегрируя (6) и учитывая (5), получим

(7)

Далее имеем

поэтому из (7) для погрешности формулы (3) получаем оценку

где

Вычисляя интеграл приходим к окончательной оценке

(8)

Погрешность составной формулы Симпсона оценивается так

(9)

Отсюда видно, что формула Симпсона существенно точнее, чем формулы прямоугольников и трапеций. На частичном отрезке она имеет точность О(h5), а на всем отрезке – O(h4)

3. Описание метода решения задачи

Для решения поставленной задачи необходимо выполнить следующие действия:

1) Ввести значения границ отрезков;

2) Вывести график функции на экран с учётом масштаба;

3) Вычислить интеграл методом трапеций;

4) Вычислить интеграл методом Симпсона;

Для успешной реализации этих действий программа должна состоять из следующих функциональных модулей: