Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

коды / 9.3

.py
Скачиваний:
12
Добавлен:
20.06.2023
Размер:
503 б
Скачать
import numpy as np

a= -5
b= 5
e= 0.01
def f(x): return x**3+4.88*x**2-2.93*x-14.149;
 
x1= a
x4= b
z= (3-np.sqrt(5))/2
while True:
    x2= x1+z*(x4-x1)
    x3= x4-z*(x4-x1)
    f2= f(x2)
    f3= f(x3)
    if f2 < f3:
        x4= x3
        x3= x2
        f3= f2
        x2= x1+z*(x4-x1)
        f2= f(x2)
    else: 
        x1= x2
        x2= x3
        f2= f3
        x3= x4-z*(x4-x1)
        f3= f(x3)
    if abs(x4-x1) <= 2*e:
        break

x= (x1+x4)/2
print(x,f(x))
Соседние файлы в папке коды