Метод деления отрезка пополам

from math import pi, tan, cos, degrees
 
eps = .0001
al, d, s, v, h, g = 0, 2*eps, 10, 12, 2, 9.81
f = lambda a: s*tan(a) - (g*s**2)/(2*(v*cos(a))**2) - h
while al<pi/2:
    if f(al) * f(al+d) <= 0:
        print(degrees(al+eps))
    al += d

x = 7.3394
y = 3.6697
71.57
45.0

from math import degrees, cos, pi, tan
 
vo = 12
s = float(input('x = '))
h = float(input('y = '))
g = 9.81
eps = 1e-6
 
left = pi / 4
right = pi / 2
alpha = (right + left) / 2
y = lambda x: s * tan(x) - g * (s / (vo * cos(x))) ** 2 / 2
res = []
while abs(y(alpha) - h) > eps:
    if y(alpha) > h:
        res.append(True)
        left = alpha
    else:
        res.append(False)
        right = alpha
    if len(res) == 2 and res[0] == res[1]:
        res = None
        print('Не добросить')
        break
    alpha = (right + left) / 2
if res is not None:
    print(round(degrees(alpha), 2))
 
if res is not None:
    left = -pi / 2
    right = pi / 3
    alpha = (right + left) / 2
    res = []
    while abs(y(alpha) - h) > eps:
        if y(alpha) > h:
            res.append(True)
            right = alpha
        else:
            res.append(False)
            left = alpha
        alpha = (right + left) / 2
    print(round(degrees(alpha), 2))

from math import pi, tan, cos, degrees
 
eps,stp = 1e-6,.1
al,d,s,v,h,g,arr,out = 0,2*eps,10,12,2,9.81,[],[]
sign = lambda x : (x > 0) - (x < 0)
f = lambda x: s*tan(x) - (g*s**2)/(2*(v*cos(x))**2) - h
# Отделяем отрезки с корнями (грубо, с шагом stp):
while al<pi/2:
    if f(al) * f(al+stp) <= 0:
        arr.append((al,al+stp))
    al += stp
#print(arr)
 
# Уточняем корни в выбранных отрезках (бисекции):
for a,b in arr:
    while b-a>d:
        c = (a+b)/2
        if sign(f(a)) != sign(f(c)):
            b = c
        else:
            a = c
    out.append(degrees((a+b)/2))
print(*out,sep='\n')