旧いです。以下へ移動をお願いします。
(その2/2)
①もっと短くなるような気もします。
②数値を代入せず、方程式で解けそうな気もしてきました。
③実数3.16227766でなく、ルート√の使えるCADはありますか。
④作図(文字を含めて)のおすすめのページを教えて下さい。
(参考)
from sympy import *
var('R TA CD')
R=10/2
TA=mpmath.radians(0)
CD=6
CI=Circle(Point(0,0),R)
A=Point(R*cos(TA), R*sin(TA))
C=Point(-R*cos(TA), -R*sin(TA))
CJ=Circle(C,CD)
result=CI.intersection(CJ)
D=Point(result[0].x,result[0].y)
AD=A.distance(D)
D2=Point(R*cos(TA)-AD, 0)
H2=D.midpoint(D2)
H =Segment(A, H2).projection(D)
result=CI.intersection(Line(D, H))
B=Point(result[1].x,result[1].y)
E= Line(A, C).intersection(Line(D, B))
print("#BE=",B.distance(E[0]))
#BE= sqrt(10)