裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

算額(その253)

2023年05月31日 | Julia

算額(その253)

中村信弥「改訂増補 長野県の算額」
http://www.wasan.jp/zoho/zoho.html
県内の算額4(199)
長野県北佐久郡軽井沢町峠 熊野神社 安政4年(1857)

鉤股弦において,鉤と股の積が 52440 平方寸,股と弦の積が 130065 平方寸のとき,鉤,股,弦それぞれの長さを求めよ。

以下の連立方程式を解く。

include("julia-source.txt");

using SymPy

@syms 鉤::positive, 股::positive, 弦::positive;

eq1 = 鉤 * 股 - 52440
eq2 = 股 * 弦 - 130065
eq3 = 鉤^2 + 股^2 - 弦^2
res = solve([eq1, eq2, eq3], (鉤, 股, 弦))

   1-element Vector{Tuple{Sym, Sym, Sym}}:
    (152, 345, 377)

鉤 = 152, 股 = 345, 弦 = 377 である。

using Plots
using Printf

function draw(more)
    pyplot(size=(500, 500), grid=false, aspectratio=1, label="", fontfamily="IPAMincho")
   (鉤, 股, 弦) = (152, 345, 377)
   plot([0, 股, 0, 0], [0, 0, 鉤, 0], linecolor=:black, linewidth=0.5)
   if more == true
       point(0, 鉤/2, " 鉤 = 152寸", :green, :left, :bottom, mark=false)
       point(股/2, 0, " 股 = 345寸", :blue, :center, :bottom, mark=false)
       point(股/2, 鉤/2, " 弦 = 377寸", :red, :left, :bottom, mark=false)
   else
      plot!(showaxis=false)
   end
end;


コメント    この記事についてブログを書く
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする
« 算額(その252) | トップ | 算額(その254) »
最新の画像もっと見る

コメントを投稿

Julia」カテゴリの最新記事