Лабораторная работа №5 По дисциплине: Вычислительная математика

Лабораторная работа № 5

По дисциплине: Вычислительная математика

2011
Написать программу для нахождения максимального значения функции на отрезке [0, 0.5] методом золотого сечения с точностью 0.0001. Считается, что требуемая точность достигнута, если выполняется условие, (e – заданная точность, ak, bk – границы интервала неопределенности, k = 0,1,2,… ), при этом, ,

N – последняя цифра пароля = 1.
Программа

Максимум найден не верно. Минимум искать не требовалось.

Исходный код программы
program Lab5;
uses

crt;

N=1;

b0=0.5;

{функция из условия}

function f(x:real):real;

begin

end;
var

x1,x2,f1,f2,b,a,x:real;

begin

a:=a0; {начальный интервал}

b:=b0;

x1:=0.618*a+0.382*b; {выбираю точки на интервале}

f1:=f(x1); {нахожу в точках значения}

f2:=f(x2);

while (b-a)>e do {цикл поиска максимума с погрешность е}

begin

if f1>f2 then {функция возрастает}

b:=x1; {передвигаю границы интервала}Не туда

F2:=f1;

x2:=x1;

f1:=f(x1); {ее значение}

end

else

a:=x2; {передвигаю границы интервала} Не туда

f1:=f2;

x1:=x2;

f2:=f(x2); {ее значение}

end;

end;

writeln('Максимум:');

writeln('x=',x:3:3); {вывожу результат}

writeln('y=',f(x):3:3);

{поиск минимума такой же только меняю знак сравнения}

a:=a0;

x1:=0.618*a+0.382*b;

x2:=0.382*a+0.618*b;

f1:=f(x1);

f2:=f(x2);

while (b-a)>e do

begin

if f1

begin

b:=x1;

x2:=x1;

f1:=f(x1);

end

else

a:=x2;

x1:=x2;

f2:=f(x2);

end;

end;

writeln;

writeln('x=',x:3:3);