算額(その5)
四十一 群馬県高崎市下小鳥町 幸宮神社 文政7年(1824)
群馬県和算研究会:群馬の算額,上武印刷株式会社,高崎市,1987年3月31日.
群馬県高崎市 幸宮神社 文政7年(1824)
https://blog.goo.ne.jp/r-de-r/e/95e0d587bc8ce1d15512697ad8ad6102
キーワード:円4個,外円,正六角形3個
#Julia, #SymPy, #算額, #和算
東京都府中市の大国魂神社にも同じ問題(上下が反対)の算額があるが,それは明治 18 年(1885) のものだ。
大円に接する 3 個の正六角形と,大円および正六角形に接する 3 個の円がある。
大円,小円の半径をそれぞれ r, x とする。x を求めよ。
AB = AC + CB = √3/4 \* r + x,OB = r - x とすると,AB / OB = cos(π/6) である。
using SymPy
@syms r::positive, x::positive;
eq1 = (r*sqrt(Sym(3))/4+x) / (r - x) - cos(PI/6);
res = solve(eq1, x);
res[1] |> factor |> println
r*(-3 + 2*sqrt(3))/2
小円の半径は x = r*(2*√3 - 3)/2 である。
include("julia-source.txt");
function hexagon(ox, oy, ol; color=:blue)
θ = range(pi/6, pi*13/6, length=7)
sinθ = ol.*sin.(θ).+oy
cosθ = ol.*cos.(θ).+ox
plot!(cosθ, sinθ, color=color, lw=0.5)
end;
function draw(r=1, more=true)
pyplot(size=(500, 500), grid=false, aspectratio=1, label="", fontfamily="IPAMincho")
pi6 = pi/6
plot()
circle2(0, 0, r, :black)
hexagon(0.5r*cos(pi6), 0.5r*sin(pi6), 0.5r)
hexagon(0.5r*cos(pi6+4pi6), 0.5r*sin(pi6+4pi6), 0.5r)
hexagon(0.5r*cos(pi6+8pi6), 0.5r*sin(pi6+8pi6), 0.5r)
x = (2*sqrt(3) - 3)r/2
circle2(0, r-x, x)
circle2((r-x)*cos(7*pi6), (r-x)*sin(7*pi6), x)
circle2((r-x)*cos(11*pi6), (r-x)*sin(11*pi6), x)
if more
delta = (fontheight = (ylims()[2]- ylims()[1]) / 500 * 10 * 2) /3 # size[2] * fontsize * 2
hline!([0], color=:gray80, lw=0.5)
vline!([0], color=:gray80, lw=0.5)
plot!([0, r*cos(pi6)], [0, r*sin(pi6)])
plot!([0, 0], [0, r-x])
plot!([0, (r-x)/2*cos(pi6)], [r-x, (r-x)/2*sin(pi6)])
annotate!(0, 0, text("O ", 10, :right, :bottom))
annotate!((r-x)/2*cos(pi6), (r-x)/2*sin(pi6), text("A", 10, :left, :top))
annotate!(0, r-x, text(" B", 10, :left, :bottom))
annotate!(x*sin(pi6), r-(x+x*cos(pi6)), text(" C", 10, :left, :top))
end
end
draw(1, true)