Курсовая работа: Отыскание корня уравнения методом половинного деления
Условием окончания поиска корня является (b-a)/2n <E или |f(x)|<E, где Е – точность, [a;b] – начальный отрезок неопределённости, n – число итераций
4. Проверка условий сходимости методов
Интегрирование по методу Симпсона
Для вычисления по методу Симпсона требуется, чтобы функция была непрерывной на отрезке интегрирования.
sin(t)/t=1 при t=0 по первому замечательному пределу, однако при вычислении в QBasic будет выдавать ошибку деления на ноль, поэтом в точке t=0 приравняем 1 искусственно.
Условие для вычисления по методу Симпсона выполняется.
Отыскание корня нелинейного уравнения методом половинного деления:
Условие f(a)*f(b)<0 выполняется при x=[0;5],
Условие единственности корня (sign (f `(x))= const при x=[a;b]) выполняется при xe [0;3]. На отрезке x=[3.3; 2*р] есть ещё один корень, поэтому сократим отрезок неопределённости до [0; 3] (с учетом условия метода Симпсона).
5. Тестирование программных модулей
5.1 Тестирование модуля численного интегрирования
Для тестирования модуля, вызовем его для отыскания интеграла
5.1.1 Схема алгоритма тестирующей программы:
Схема алгоритма управляющей программы:
Схема алгоритма f(x) :
Схема алгоритма модуля численного интегрирования при тестировании:
5.1.2 Код тестирующей программы :
DECLARE function integr (afix,x,E)
DECLARE FUNCTION fint (x)
CLS
PRINT "Itog"; E; integr(0,1,0.001)
END
FUNCTION fint (t)