Подправьте код: половинное деление и золотое сечение

@Matan! · Регистрация: 31.05.2013

Author24 — интернет-сервис помощи студентам

Доброго всем времени суток!
Подправьте, пожалуйста,код: результат первого и вторых способов должны совпадать.

Matlab M

restart;
f := proc (x) options operator, arrow; abs(x^5)+(x-1)*cos(x)+2 end proc; a := -1; b := 1; c := (a+b)*(1/2); eps2 := 10^(-4); plot(f(x), x = a .. b, color = red, thickness = 1); with(Optimization); LSSolve([f(x)], {a <= x, x <= b});
 
eps := (1/5)*eps2; xx := (a+b)*(1/2); del := b-a; ax := a; bx := b; while del > 2*eps2 do if evalf(f(xx-eps)) > evalf(f(xx+eps)) then ax := xx-eps; bx := bx else ax := ax; bx := xx+eps end if; xx := evalf((ax+bx)*(1/2)); del := bx-ax end do; xx := (ax+bx)*(1/2); f(xx);
 
NULL;
;
tau := evalf((-1+sqrt(5))*(1/2)); a1 := a; b1 := b; c := a1+(1-tau)*(b1-a1); d := a1+tau*(b1-a1); del := b1-a1; xx := (c+d)*(1/2);
while del > eps2 do if evalf(f(c)) > evalf(f(d)) then a1 := c; b1 := b1; c := a1+(1-tau)*(b1-a1); d := a1+tau*(b1-a1) else a1 := a1; b1 := d end if; xx := (c+d)*(1/2); del := b1-a1; c := a1+(1-tau)*(b1-a1); d := a1+tau*(b1-a1) end do;
koren := (c+d)*(1/2);
f(koren);

Привожу файл,чтобы понятней было:
Дихотомия и золотое сечение.zip

@VSI · 14.05.2014, 19:31

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

@Matan! · 14.05.2014, 22:47 **[ТС]**

Сообщение от VSI

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

Спасибо, но мне интересно,почему же в данном коде неправильно работает прога...

@VSI Модератор $Эксперт по математике/физике$ 5240 / 4027 / 1385 Регистрация: 30.07.2012 Сообщений: 12,288
	14.05.2014, 19:31	2
	Метод половинного деления... Метод золотого сечения... 0

@Matan! 1437 / 1014 / 228 Регистрация: 31.05.2013 Сообщений: 6,645 Записей в блоге: 6
	14.05.2014, 22:47 [ТС]	3
	Сообщение от VSI Метод половинного деления... Метод золотого сечения... Спасибо, но мне интересно,почему же в данном коде неправильно работает прога... 0