算額（その1123）

算額（その1123）

二十八 一関市萩荘 赤萩観音寺 天保2年(1831)
山村善夫：現存 岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.
http://www.wasan.jp/yamamura/yamamura.html
キーワード：円6個，長方形

長方形の中に甲円 1 個，乙円 1 個，丙円 2 個，丁円 2 個を容れる。長方形の短辺の長さが 1 寸のとき，長辺の長さはいかほどか。

「岩手の算額」も，解けるものは解いてきたが，解く前から「こりゃー，やってみても解けるかどうかわからない」というのが残っている。この算額の問題も，図を見ただけで，「なんだこりゃ」と後回しにしてきた。こんなの算額の図じゃない。案の定，以下のような図になるものであった。しかも，これもこの本に載っている問題の常で，答えが怪しい。算額の元々がおかしいのか，山村の解説（？）も首をひねるものが多い。

このシリーズでは，得られた解で図を描いて提示している。図が正しく描かれているということは，解も正しいという一つの証拠ではあろう。
多くの算額問題では，条件を与えて結果としてパラメータの一つを求めるというのが常道である。しかし，一つのパラメータ以外のパラメータも含めて図を描いて検証するということは殆ど行われてこなかったのではないかと思う次第である。

長方形の長辺，短辺を a, 2r1
甲円の半径と中心座標を r1, (a - r1, r1)
乙円の半径と中心座標を r2, (r2, r1)
丙円の半径と中心座標を r3, (x3, r3)
丁円の半径と中心座標を r4, (r4, r4)
とおき，以下の連立方程式を解く。

include("julia-source.txt")

using SymPy
@syms a::positive, r1::positive, r2::positive,
     r3::positive, x3::positive,
     r4::positive
eq1 = (a - r1 - x3)^2 + (r1 - r3)^2 - (r1 + r3)^2
eq2 = (r2 - x3)^2 + (r1 - r3)^2 - (r2 + r3)^2
eq3 = (r2 - r4)^2 + (r1 -r4)^2 - (r2 + r4)^2
eq4 = (x3 - r4)^2 + (r3 - r4)^2 - (r3 + r4)^2
eq5 = 2r1 + 2r2 - a
res = solve([eq1, eq2, eq3, eq4, eq5], (a, r2, r3, x3, r4))[1]

   (4*r1*(-6*sqrt(239 + 169*sqrt(2)) + 4 + 3*sqrt(2) + 4*sqrt(478 + 338*sqrt(2))), r1*(-12*sqrt(239 + 169*sqrt(2)) + 7 + 6*sqrt(2) + 8*sqrt(478 + 338*sqrt(2))), 2*r1*(-145*sqrt(239 + 169*sqrt(2)) - 100*sqrt(478 + 338*sqrt(2)) + 3131 + 2214*sqrt(2))/(-44*sqrt(239 + 169*sqrt(2)) - 31*sqrt(2)*sqrt(239 + 169*sqrt(2)) + 960 + 679*sqrt(2)), r1*(-2085*sqrt(478 + 338*sqrt(2)) - 2944*sqrt(239 + 169*sqrt(2)) + 64416 + 45549*sqrt(2))/(-44*sqrt(239 + 169*sqrt(2)) - 31*sqrt(2)*sqrt(239 + 169*sqrt(2)) + 960 + 679*sqrt(2)), r1*(-sqrt(1912 + 1352*sqrt(2)) + 31 + 22*sqrt(2)))

res[1] |> println
res[1](r1 => 1/2).evalf() |> println

   4*r1*(-6*sqrt(239 + 169*sqrt(2)) + 4 + 3*sqrt(2) + 4*sqrt(478 + 338*sqrt(2)))
   1.48071257037594

長方形の一辺の長さは甲円の半径の 4(-6*sqrt(239 + 169√2) + 4 + 3√2 + 4sqrt(478 + 338√2)) 倍である。
これはこれ以上簡約化できない。

長方形の短辺が 1 寸のとき，長辺は 1.48071 寸である。

「術」では，1.481136081497563 倍としているが，数値は微妙に違うし，山村は 術の 2 を 2.21 に変更して（理論が不明であるし），「答」に合わせているが，それは本末転倒である。

天 = √2 + 3
(2.21 - √天)*天 + 1

   1.481136081497563

その他のパラメータは以下の通りである。

  r1 = 0.5;  a = 1.48071;  r2 = 0.240356;  r3 = 0.153333;  x3 = 0.426939;  r4 = 0.137032

function draw(r1, more=false)
   pyplot(size=(500, 500), grid=false, aspectratio=1, label="", fontfamily="IPAMincho")
   (a, r2, r3, x3, r4) = (4*r1*(-6*sqrt(239 + 169*sqrt(2)) + 4 + 3*sqrt(2) + 4*sqrt(478 + 338*sqrt(2))), r1*(-12*sqrt(239 + 169*sqrt(2)) + 7 + 6*sqrt(2) + 8*sqrt(478 + 338*sqrt(2))), 2*r1*(-145*sqrt(239 + 169*sqrt(2)) - 100*sqrt(478 + 338*sqrt(2)) + 3131 + 2214*sqrt(2))/(-44*sqrt(239 + 169*sqrt(2)) - 31*sqrt(2)*sqrt(239 + 169*sqrt(2)) + 960 + 679*sqrt(2)), r1*(-2085*sqrt(478 + 338*sqrt(2)) - 2944*sqrt(239 + 169*sqrt(2)) + 64416 + 45549*sqrt(2))/(-44*sqrt(239 + 169*sqrt(2)) - 31*sqrt(2)*sqrt(239 + 169*sqrt(2)) + 960 + 679*sqrt(2)), r1*(-sqrt(1912 + 1352*sqrt(2)) + 31 + 22*sqrt(2)))
   @printf("長方形の短辺が %g のとき，長辺は %g である。\n", 2r1, a)
   @printf("r1 = %g;  a = %g;  r2 = %g;  r3 = %g;  x3 = %g;  r4 = %g\n", r1, a, r2, r3, x3, r4 )
   plot([0, a, a, 0, 0], [0, 0, 2r1, 2r1, 0], color=:black, lw=0.5)
   circle(a - r1, r1, r1)
   circle(r2, r1, r2, :green)
   circle(x3, r3, r3, :orange)
   circle(x3, 2r1 - r3, r3, :orange)
   circle(r4, r4, r4, :magenta)
   circle(r4, 2r1 - r4, r4, :magenta)
   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)
       point(a, 0, " a", :black, :left, :bottom, delta=delta/2)
       point(a - r1, r1, "甲円:r1,(a-r1,r1)", :red, :center, delta=-delta/2)
       point(r2, r1, "乙円:r2,(r2,r1)", :green, :center, delta=-delta/2)
       point(x3, r3, "丙円:r3\n(x3,r3)", :orange, :center, delta=-delta/2)
       point(r4, r4, "丁円:r4\n(r4,r4)", :magenta, :center, delta=-delta/2)
       point(a, 2r1, "(a,2r1)", :black, :right, :bottom, delta=delta/2)
   end
end;