シムソン線 - DyGeom.m on Mathematica+KETpic

(* 0.40 以降 *)
SetDirectory[NotebookDirectory[]];
<< DyGeom.m
pa = 点[{0, 0}, pointname -> "A", pointposition -> "w"]
pb = 点[{4, 0}, pointname -> "B", pointposition -> "e2"]
pc = 点[{3, 3}, pointname -> "C", pointposition -> "n2"]
lbc = 直線[pb, pc, clipping -> {pb, pc}]
lab = 直線[pa, pb, clipping -> {pa, pb}]
lca = 直線[pc, pa, clipping -> {pc, pa}]
pp = 中点[pb, pc]
lp = 垂線[pp, lbc]
pq = 中点[pa, pb]
lq = 垂線[pq, lab]
po = 交点[lp, lq]
c1 = 円[po, pa]
pr = 円上の点[2.5, -3, c1, pointname -> "R", pointposition -> "s2"]
ld = 垂線[pr, lbc]
pd = 交点[ld, lbc, pointname -> "D", pointposition -> "ne2"]
ld = AddOption[ld, clipping -> {pr, pd}]
le = 垂線[pr, lca]
pe = 交点[le, lca, pointname -> "E", pointposition -> "n2"]
le = AddOption[le, clipping -> {pr, pe}]
lf = 垂線[pr, lab]
pf = 交点[lf, lab, pointname -> "F", pointposition -> "sw2"]
lf = AddOption[lf, clipping -> {pr, pf}]
mk1 = 角度印[pr, pd, pc, markshape -> "rightangle"]
mk2 = 角度印[pr, pe, pc, markshape -> "rightangle"]
mk3 = 角度印[pr, pf, pb, markshape -> "rightangle"]
simsonline = 直線[pd, pe]
lbc = AddOption[lbc, clipping -> {pc, pd}]
作図[pa, pb, pc, lab, lbc, lca, po, c1, pr, ld, pd, le, pe, lf, pf, simsonline, mk1, mk2, mk3]

---