Метод половинного деления (дихотомии)

@Ильдар Басков · Регистрация: 22.11.2014

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

Функция f(x)	Интервал [a; b]	Начальное значение шага
(1-x)⁴	[0.2; 1.5]	0.25

Написать программу, метод равномерного поиска.
Я накидал тут, но у меня не выводит результат вообще. Писал на ABC Pascal, написал потом на Turbo Pascal, программа запускается но при вводе начального шага зависает и все, и больше никаких движений, в чем ошибка помогите пожалуйста. Метод половинного деления(дихотомии).
Вот программа написанная мною:

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
Program pr_7;
uses crt;
function f(x:real):real;
begin
f:=EXP(4*LN((1-x)));
end;
var a,b,e,c,x:real;
begin
a:=0.2;
b:=1.5;
write ('e=');
read(e);
c:=(a+b)/2;
while abs(b-a)>e do
begin
if f(a)*f(c)<0 then
b:=c
else
a:=c;
c:=(a+b)/2
end;
x:=(a+b)/2;
writeln ('При X=0.2',' Y=',a);
writeln ('При X=1.5',' Y=',b);
end.

@JokeR.BY · 24.11.2014, 14:41

потому что прежде чем копировать будь с какого сайта (код ведь не ваш, не так ли?), нужно внимательно посмотреть на свою формулу и вспомнить на каком отрезке определены функции ln(x) и ln(1-x). для второго: при x < 1. что в ходе расчёта явно не выполняется. поэтому:

Pascal
1
f:=EXP(4*LN((1-x))); <-> f:=sqr(sqr(1-x));

Добавлено через 2 минуты
и да: зачем ~~усы сбрил~~ поменяли последние строки?

@Ильдар Басков · 24.11.2014, 17:07 **[ТС]**

Во-первых не откуда не копировал, во-вторых не до конца понял, степень разве не выставляется?) Смотрел через таблицу, там написано, что степень так выставлять, f:=EXP(4*LN((1-x)));, может я что то и напутал. Извиняюсь. Мне же скобочку в ^4 степень надо возвести.
Разве не так? f:=EXP(4*LN((1-x))); Вот прикладываю табличку:

@Puporev · 24.11.2014, 17:10

А если x>=1?

Добавлено через 1 минуту
Эта формула exp(y*ln(x)) напрямую работает только при x>0

@Ильдар Басков · 24.11.2014, 17:21 **[ТС]**

Получается данная формула f:=sqr(sqr(1-x)); в математике выглядит приблизительно 2*(1-x)^2 так?
Это получается что если x будет =0 то функция не получится? и использовать можно только f:=sqr(sqr(1-x)) это?

@JokeR.BY · 24.11.2014, 17:22

Ильдар Басков, я могу привести сайт (первый в выборке гугла по запросу "метод половинного деления паскаль"), откуда вы скопировали строка в строку (поменяв только вывод и, собственно, формулу). не надо по ушам ездить.
а по теме:
дважды возведение в квадрат - это не ^4 разве?

@Ильдар Басков · 24.11.2014, 17:24 **[ТС]**

я уже понял) ((1-x)^2)^2

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11680&d=1772460536 Одним из. . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .
SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .

@Ильдар Басков 0 / 0 / 0 Регистрация: 22.11.2014 Сообщений: 22
	24.11.2014, 17:07 [ТС]
	Во-первых не откуда не копировал, во-вторых не до конца понял, степень разве не выставляется?) Смотрел через таблицу, там написано, что степень так выставлять, f:=EXP(4LN((1-x)));, может я что то и напутал. Извиняюсь. Мне же скобочку в ^4 степень надо возвести. Разве не так? f:=EXP(4LN((1-x))); Вот прикладываю табличку: 0

@Puporev Почетный модератор 64315 / 47611 / 32743 Регистрация: 18.05.2008 Сообщений: 115,167
	24.11.2014, 17:10
	А если x>=1? Добавлено через 1 минуту Эта формула exp(y*ln(x)) напрямую работает только при x>0 1

@Ильдар Басков 0 / 0 / 0 Регистрация: 22.11.2014 Сообщений: 22
	24.11.2014, 17:21 [ТС]
	Получается данная формула f:=sqr(sqr(1-x)); в математике выглядит приблизительно 2*(1-x)^2 так? Это получается что если x будет =0 то функция не получится? и использовать можно только f:=sqr(sqr(1-x)) это? 0

@JokeR.BY CAPITAL OF ROCK! 1281 / 708 / 982 Регистрация: 03.03.2010 Сообщений: 2,286
	24.11.2014, 17:22
	Ильдар Басков, я могу привести сайт (первый в выборке гугла по запросу "метод половинного деления паскаль"), откуда вы скопировали строка в строку (поменяв только вывод и, собственно, формулу). не надо по ушам ездить. а по теме: дважды возведение в квадрат - это не ^4 разве? 1

@Ильдар Басков 0 / 0 / 0 Регистрация: 22.11.2014 Сообщений: 22
	24.11.2014, 17:24 [ТС]
	я уже понял) ((1-x)^2)^2 0