-
Notifications
You must be signed in to change notification settings - Fork 0
/
solve1dmonot.m
44 lines (40 loc) · 1.05 KB
/
solve1dmonot.m
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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
# Принимает на вход монотонно возрастающую функцию F
# x0 - начальное приближение
# x_eps - допустимая точность x
# F_eps - допустимая точность F
# R - максимальное значение x
function sol = solve1dmonot (F, x0, x_eps, F_eps, R)
sol = [];
# вычисляем нижнюю границу для корня
left = x0;
x_step = x_eps;
while (F(left) > 0 && abs(left) < 2 * R)
left -= x_step;
x_step *= 2;
endwhile
if (F(left) > 0)
#disp("error")
return
endif
# вычисляем верхнюю границу для корня
right = x0;
x_step = x_eps;
while (F(right) < 0 && abs(right) < 2 * R)
right += x_step;
x_step *= 2;
endwhile
if (F(right) < 0)
#disp("error")
return
endif
# F(left) <= 0, F(right) >= 0
while (right - left > x_eps)
mid = (left + right) / 2;
if (F(mid) > 0)
right = mid;
else
left = mid;
endif
endwhile
sol = left;
endfunction