①計算1と計算2の使い方の違いを教えて下さい。
②キーインを少なくしたいので、piだけでできると助かります。
③(追加)ルートの正解は、sympyで、できるはず?
sympyで入試問題に挑戦!
解説は、理解できてません。方べきの定理ですか?
(参考)入試問題に挑戦! 平成17年宮城県(超難)
http://www.5kaku.net/?p=9299
(参考)WolframAlphaで
28/3 π - 3 sqrt(3)
https://www.wolframalpha.com/input/?i=28%2F3+%C3%97+%CF%80+%E2%80%93+3%E2%88%9A3+&lang=ja
import math
from sympy import *
def my_sanArea(a,b,deg):
return 1/2*a*b*sin(deg*math.pi/180)
def my_yumArea(r,deg):
return 1/2*r**2*(deg*math.pi/180-sin(deg*math.pi/180))
g=8.0
L=2.0
# t=(180-60)*pi/180 #計算1
t=(180-60)*math.pi/180 #計算2
#
P0= Point(0.0, 0.0)
P1= P0+Point(-g/2, 0.0)
P2= P0+Point( g/2, 0.0)
P3= P0+Point(L*cos(t), L*sin(t))
C123 = Circle(P1, P2, P3)
P0P2 = P0.distance(P2)
P1P3 = P1.distance(P3)
AS1 = my_sanArea(g/2,L,60.0)
AR1 = my_yumArea(C123.radius,math.asin(P1P3/2/C123.radius)*180/math.pi*2)
#
result = C123.intersection(Line(P3, P0))
PK = Point(result[1][0], result[1][1])
P0PK = P0.distance(PK)
P2PK = P2.distance(PK)
AS2 = my_sanArea(P0P2,P0PK,60.0)
AR2 = my_yumArea(C123.radius,math.asin(P2PK/2/C123.radius)*180/math.pi*2)
#
print('面積=',AS1+AR1+AS2+AR2)
print('面積=',float(AS1+AR1+AS2+AR2))
# 正解 28/3 × π – 3√3 c㎡
# wolfram 面積= 24.125379010798105011735665886091076485011831966069845778265395256...
# AutoCAD 面積= 2412.53790108 mm2
#
# 計算1の結果 面積= 24.1253790107981
# 計算1の結果 面積= 24.125379010798106
# 計算2の結果 面積= 10.2689725502471 + 5.54256258422039e-13*sqrt(625000000000001464921241321)
# 計算2の結果 面積= 24.12537901079808