одномерная оптимизация(золотое сечение)

#Целевая функция
def f(x):
    return x**2-4*x

#Начальный отрезок, содержащий минимум
a = 0
b = 3.5

Phi = 1.618#пропорция
acc=0.00000001
y = []

#Алгоритм
while 1:
    x_1=b - ((b - a)/Phi)
    x_2=a + ((b-a)/Phi)
    if f(x_1) >= f(x_2):
        a = x_1
    else:
        b = x_2
    if abs(b-a) < acc:
        print(f((a+b)/2), (a+b)/2)
        break