Реализация методов оптимизации (метод касательных и метод парабол)

import math
from numpy import *
def func(x):
     return math.cos(x)-math.sqrt(x);
print ("Введите 4 действительных числа, a c b Eps, после каждого введенного числа нажимайте клавишу Enter:")
a = int(input())
c = int(input())
b = int(input())
Eps = float(input())
 
print (" Вы ввели",a,c,b,Eps)
ya=func(a)
yc=func(c)
yb=func(b)
t=c+1/2*(((b-c)*(b-c))*(ya-yc)-((c-a)*(c-a))*(yb-yc))/((b-c)*(ya-yc)+(c-a)*(yb-yc))
if t!=c:
    x=t
else:
    x=(a+c)/2        
y=func(x)
k=0
print(t,x,y)
while (b-a)>(2*Eps): 
    if x<c:
        if y<yc:
            b=c
            c=x
            yb=yc
            yc=y
        elif y>yc:
            a=x
            ya=y
        else :
            a=x
            b=c
            c=(x+c)/2
            ya=y
            yb=yc
            yc=func(c)
    elif x>c:
        if y<yc:
            a=c
            c=x
            ya=yc
            yc=y
        elif y>yc:
            b=x
            yb=y
        else :
            a=c
            b=x
            c=(x+c)/2
            ya=yc
            yb=y
            yc=func(c)
    else:
        break
    k=k+1
Xsr=(a+b)/2
Ysr=func(Xsr)
print("Среднее значение корня равно {0:.6f}, значении функции {1:.6f}, достигнуто за {2:} шагов".format(Xsr,Ysr,k))

from sympy import *
x, y, z, a, b = symbols('x y z a b')
import math
from numpy import *
def func(x):
     return math.cos(x)-math.sqrt(x);
print ("Введите 3 действительных числа, a b Eps после каждого введенного числа нажимайте клавишу Enter:")
x = Symbol('x')
a = int(input())
b = int(input())
Eps = float(input())
print (" Вы ввели",a,b,Eps)
 
y1=func(a)
y2=func(b)
 
str_expr = "diff(cos(x)-sqrt(x),x)" #переделываем строку в формулу
expr = sympify(str_expr)
z1=(expr.subs(x, a)).evalf()
z2=(expr.subs(x, b)).evalf()
print (z1,z2)
k=0
T=0
while b-a>(2*Eps):
    c=((b*z2-a*z1)-(y2-y1))/(z2-z1)
    c=round(c, 6)#округ
    y=func(c)
    y=round(y, 6)
    z=(expr.subs(x, c)).evalf()
    z=round(z, 6)
    print(c,y,z)
    if z==0:
        Xsr=c
        Ysr=y
        T=1
        break
    else:
        if z<0:
            a=c
            y1=y
            z2=z
        else:
            b=c
            y2=y
            z2=z
    k=k+1
if T==0:
    Xsr=(a+b)/2
    Ysr=func(Xsr)
print("Среднее значение корня равно {0:.6f}, значении функции {1:.6f}, достигнуто за {2:} шагов".format(Xsr,Ysr,k))