Много интересна задача с това, че ще се опитам да я реша като използвам числени методи:
Ще изпозлваме много отличния чертеж на Hephaestus.
Цялото решение ще стъпи на тоеа, че ще търсим CD като функция на ъгъла при върха А и когато тази финкция $f$ върне 5, ще сме готови. За тази цел ще използваме една косинусова и 2 синусови теореми:
- Код: Избери целия код
In [93]: def f(A):
...: a = sqrt(10**2 + 11**2 - 2*10*11*cos(A))
...: x = asin(11*(sin(A)/a))/3
...: y = pi -(A + x)
...: CD = sin(x)*10/sin(y)
...: return CD
...:
In [94]: float(f(pi/3))
Out[94]: 3.7167888629845
In [95]: float(f(pi/4))
Out[95]: 4.453439699192902
In [96]: float(f(pi/5))
Out[96]: 5.0666706589651564
Добре,виждаме, че отговора е някъде около $\frac{pi}{5}$
Сега ще запалим scipy solvera:
- Код: Избери целия код
In [102]: from scipy.optimize import fsolve
In [107]: fsolve(lambda t: f(t)- 5, float(pi/5))
Out[107]: array([ 0.64350111])
In [108]: float(10*11*sin(0.64350111)/2)
Out[108]: 33.00000005309549
Ок и за да намерим лицето изпозлвахме формулата по две страни и синус на ъгъл межди тях и получихме 33.00000005309549, което дава грешка около 1е-7. Като за числов метод мисля що годе добре се справи, макар, че се надявах на малко по-добър резултат. Явно някъде загубихме малко точност около плаващита запетаи.