用Mathematica演示三角形外接圆Simson线的包络

用Mathematica演示三角形外接圆Simson线的包络

c=Y Cos[1/2 (-h+u+v+w)]-(Cos[(v-w)/2] Sin[(h-u)/2]+Cos[(h-u)/2] Sin[h-v/2-w/2])-X Sin[1/2 (-h+u+v+w)]==0;
(*Subscript[当\[Theta], 3]增加微小变量t的时候，Simson线也变成了另一条Simson线，这两条Simson线有交点TT*)
d=c/.h->h+t
sol=Solve[{c,d},{X,Y}]//FullSimplify;
TT={X,Y}/.sol;
(*当微小变量t趋向于0的时候，TT的极限状态，就是Simson线与包络线的切点，这个切点的坐标也可以视为包络线的参数方程*)
tt=Limit[TT,t->0]//FullSimplify;
Style[tt//Flatten//TraditionalForm,Blue,20]
Simson=FullSimplify[c/.{u->0,w->(5 \[Pi])/3,v->(2 \[Pi])/3,h->xx}]
ContourPlot[Table[Simson,{xx,0,2 Pi,2 Pi/100}]//Evaluate,{X,-1.1,2},{Y,-1.5,1.1}]
T=tt/.{u->0,w->(5 \[Pi])/3,v->(2 \[Pi])/3}
Points=Point[{
   {Cos[u],Sin[u]},{Cos[v],Sin[v]},{Cos[w],Sin[w]}}/.{u->0,w->(5 \[Pi])/3,v-&g