БГ
(1/(math.tan(float(angle_arr_1[i_1-1])*math.pi/180))+1/(math.tan(float(angle_arr_2[i_2-1])*math.pi/180)))/(1/(math.tan(float(angle_arr_2[i_2-1])*math.pi/180))-1/(math.tan(float(angle_arr_1[i_1-1])*math.pi/180)))*d/2
Но MATH подтормаживает
ctan1 = 1/(math.tan(float(angle_arr_1[i_1-1])*math.pi/180))
ctan2 = 1/(math.tan(float(angle_arr_2[i_2-1])*math.pi/180
(tan1+tan2)/(tan2-tan1)*d/2
Порядком ускорит код за счёт переиспользования результатов вычисления