laba_4
.docxАвтономное образовательное учреждение
высшего профессионального образования Ленинградской области
Государственный институт экономики, финансов, права и технологий
Кафедра информационных технологий и высшей математики
Дисциплина «Использование Mathcad для решения экономических задач »
Лабораторная работа №4
на тему
«Численное интегрирование»
Выполнила студентка
3 курса
123 группы
Янченко В.А.
Проверил:
Алексеев Г.В.
Гатчина
2014
Цель работы: ознакомиться с численными методами вычисления определённых интегралов, научиться решать задачи с использованием формулы Симпсона, трапеций, правых и левых прямоугольников, метода Монте-Карло и оценивать погрешность всех перечисленных формул.
Задача №1: Метод прямоугольников
Вычислить приближенное значение интеграла функции
на отрезке [3;5] , используя формулы левых и правых прямоугольников, при n=1000.
Ход работы:
-
Задаем функцию f(x), отрезок [a,b] и функцию нахождения дифференциалов n-го порядка.
В нашем случае: функция , а отрезок [3;5].
-
Находим значение интеграла заданной функции для использования его в дальнейшем решении.
Составим функцию , входными параметрами которой являются : a,b- левая и правая границы интервала ; n- количество разбиений ; char- если имеет значение ”left’’, то идет подсчет по формуле левых прямоугольников , любое другое - по формуле правых прямоугольников.
Погрешность показывает, что полученное значение интеграла верно до третьего знака после запятой.
Integr(a,b,n,’’left’’)=
Integr(a,b,n,’’right’’)=
Ответ: приближенное значение интеграла =6.669 по формуле левых прямоугольников, и - по формуле правых прямоугольников.
Задача №2: Вычислить приближенное значение интеграла функции используя общую формулу Симпсона, при n=1000.
Ход работы:
Составим функцию , входными параметрами которой являются : a,b- левая и правая границы интервала ; n- количество разбиений. Индексы iEven и iUneven обозначают четность и нечетность соответственно.
Simpson(a,b,n)=
Следовательно, решением будет число, равное 6.667. Погрешность показывает, что полученное значение интеграла верно до третьего знака.
Ответ: значение интеграла
Задача №3: Вычислить приближенное значение интеграла функции используя формулу трапеций, при n=1000.
Ход работы:
Введем функцию, реализующую вычисление интеграла методом трапеций.
Trapes(a,b,n)=
Ответ: значение интеграла
Задача №4: Вычислить приближенное значение интеграла функции используя метод Монте-Карло, при n=1000.
Ход работы:
i:=0..n
I=6.665
Ответ: значение интеграла