Вычисление интеграла методом Ньютона-Котеса (теория и программа на Паскале)
| Категория реферата: Рефераты по кибернетике
| Теги реферата: решебник 7, реферат география на тему
| Добавил(а) на сайт: Osokin.
Предыдущая страница реферата | 1 2 3 4 5 6 7 | Следующая страница реферата
3) Находим коэффициенты Ньютона –Котеса:
H0=H5=0.065972 ;H1=H4=0.260417 ;H2=H3=0.173611 ;
4)Подставим значения в формулу и получим:
[pic]
Рассмотрим частные случаи формулы Ньйтона-Котеса.
Пусть n=1 тогда
H0=H1=0.5 и конечная формула примет вид:
[pic] Тем самым в качестве частного случая нашей формулы мы получили
формулу трапеций.
Взяв n=3, мы получим
[pic] . Частный случай формулы Ньютона –Котеса – формула Симпсона
Теперь произведем анализ алгоритма и рассмотрим основной принцип работы программы.
Для вычисления интеграла сначала находятся коэффициенты Ньютона-
Котеса. Их нахождение осуществляется в процедуре hkoef.
Основной проблемой вычисления коэффициентов является интеграл от
произведения множителей. Для его расчета необходимо:
А) посчитать коэффициенты при раскрытии скобок при q
(процедура mnogoclen)
Б) домножить их на 1/n , где n –степень при q (процедура koef)
В) подставить вместо q значение n (функция integral)
Далее вычисляем факториалы (функция faktorial) и перемножаем полученные выражения (функция mainint). Для увеличения быстроты работы вводится вычисление половины от количества узлов интерполяции и последующей подстановкой их вместо неподсчитанных.
Процедура koef(w: массив;n:целый;var e:массив);
[pic]
Процедура hkoef(n:целый;var h:массив);
[pic]
[pic]
Процедура mnogochlen(n,i:целые;var c:массив );
[pic]
[pic]
Процедура funktia(n:целая;a,b:вещест.;var y:массив;c:вещест.;f:строка);
[pic]
Функция facktorial(n:целый):двойной;
[pic]
Функция integral(w:массив;n:целый):двойной;
[pic]
Функция mainint(n:целый;a,b:вещест.;y:массив):двойной;
[pic]
Основная программа
[pic]
Программа состоит из 8 файлов:
. K_main.exe – файл загрузки основной программы
. K_unit.tpu – модуль вычислительных процедур и функций
. K_graph.tpu – модуль графических процедур
. Graphic.tpu – модуль процедур для построения графика
. Egavga.bgi – файл графической инициализации
. Sans.chr, litt.chr – файлы шрифтов
. Keyrus.com (не обязательно) – файл установки русского языка.
Для работы программы с русским интерфайсом желательно запускать ее в
режиме DOS.
Рекомендуем скачать другие рефераты по теме: контрольная по русскому, решебник по математике 6 класс.
Категории:
Предыдущая страница реферата | 1 2 3 4 5 6 7 | Следующая страница реферата