«Метод половинного деления на Паскале (Pascal)»

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

Создать программный продукт, который находит искомый корень уравнения в отрезке при помощи метода половинного деления.

sin(x-0.5)-x+1=0

Метод половинного деления.

Для этого метода существенно, чтобы функция f(x) была непрерывна и ограничена в заданном интервале [a, b], внутри которого находится корень. Предполагается также, что значения функции на концах интервала f(a) и f(b) имеют разные знаки, т.е. выполняется условие f(a)f(b) .

Обозначим исходный интервал [a, b] как [a0, b0]. Для нахождения корня уравнения f(x) = 0 отрезок [a0, b0] делится пополам, т.е. вычисляется начальное приближение x0 = (a0 + b0)/2. Если f(x0) = 0, то значение x0 = x* является корнем уравнения. В противном случае выбирается один из отрезков [a0, x0] или [x0, b0], на концах которого функция f(x) имеет разные знаки, так как корень лежит в этой половине. Далее выбранный отрезок обозначается как [a1, b1], вновь делится пополам точкой x1 = (a1 + b1)/2 и т.д. В результате на некоторой итерации получается точный корень x* уравнения f(x) = 0, либо бесконечная последовательность вложенных отрезков [a0, b0], [a1, b1], ., [ai, bi], ., таких, что f(ai)f(bi)  (i =1, 2, .), сходящихся к корню x*.

4. Текст программы на Паскале

program mdp;

function f(x: real): real;

begin

.

end;

var

a, b, e, c, x: real;

begin

write('a=');

read(a);

write('b=');

read(b);

write ('e=');

read(e);

c:=(a+b)/2;

while(b-a)>e do

begin

if(a)*f(c)<0 then

b:=c

else

a:=c;

.

readln;

end.

5. Результаты расчёта

Результаты требуемого расчёта:

a=1

b=3

e=0.01

a=1.0000b=2.0000f(a)=0.479425539f(b)=-0.002505013

a=1.5000b=2.0000f(a)=0.341470985f(b)=-0.002505013

a=1.7500b=2.0000f(a)=0.198984619f(b)=-0.002505013

a=1.8750b=2.0000f(a)=0.105893057f(b)=-0.002505013

a=1.9375b=2.0000f(a)=0.053629191f(b)=-0.002505013

a=1.9688b=2.0000f(a)=0.026047790f(b)=-0.002505013

a=1.9844b=2.0000f(a)=0.011893001f(b)=-0.002505013

a=1.9922b=2.0000f(a)=0.004724417f(b)=-0.002505013

x=1.996 f(x)=0.0011

Pascal

X=. на интервале [1; 3]

6. Вывод

1. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы – М.: Лаборатория Базовых Знаний, 2002.

2. Численные методы. Автор: Лапчик М.П., Рагулина М.И., Хеннер Е.К.; под ред. Лапчика М.П.