算額，和算，数学，その他，Julia ときどき R, Python によるコンピュータプログラム，コンピュータ・サイエンス，統計学<br>
算額（その306）「三重県に現存する算額の研究」福島完(2007/2/13)https://mie-u.repo.nii.ac.jp/?action=repository_uri&item_id=7216三重県伊賀市永保寺天保15年(1844)問題文2三辺のの和が26.14丈となる鉤股弦（直角三角形）がある。鉤股の差は分からないが，弦は鉤より4.44丈長い。この直角三角形に内接するように正三角形を入れるとき，正三角形の一辺，中句，三角面を求めよ。中句と弦の交点を(x,y)，正三角形の他の2点を(x,y2),(x3,0)とおき，以下の連立方程式を解く。まず，鉤,股,弦,中句を求める。include("julia-source.txt");usingSymPy@syms鉤::positive,股::posit...算額（その306）
算額（その305）「三重県に現存する算額の研究」福島完(2007/2/13)https://mie-u.repo.nii.ac.jp/?action=repository_uri&item_id=7216三重県伊賀市菅原神社嘉永7年(1854)問題文4外大円の中に長方形を起き，菱形を切り取った残りの領域に円を置く。円の直径が3尺，長方形の残りの長さBEが4.5尺である。このとき，菱形の一辺（菱面=BD），長軸（同立=AD），短軸（同横=BC），外大円の円周，矢（r0-y）を求めよ。外大円の半径をr0とおく。r0=同立/2長方形の右上の頂点の座標を(x,y)とする。以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsr1::positive,x::p...算額（その305）
算額（その304）「三重県に現存する算額の研究」福島完(2007/2/13)https://mie-u.repo.nii.ac.jp/?action=repository_uri&item_id=7216三重県伊賀市菅原神社嘉永7年(1854)問題文1正方形に内接する円を描き，円の中に図のように菱形と2種類の円を2つずつ入れる。菱形の一辺と短軸は等しい。小円の直径が3.55尺のとき，正方形の一変，正方形に内接する円の円周，菱形の一変，菱形内の大円の直径，矢（菱形の上の頂点から正方形の上辺までの距離）を求めよ。正方形の中心を原点とする。正方形の一辺の長さを2r0とする（内接する円の半径がr0）。菱形内の小円の半径，中心座標をr1,(r1,y1)，大円の半径，中心座標をr2,(r2,0)とする。以下の連立方程...算額（その304）
算額（その303）「三重県に現存する算額の研究」福島完(2007/2/13)https://mie-u.repo.nii.ac.jp/?action=repository_uri&item_id=7216三重県伊賀市菅原神社嘉永7年(1854)問題文3外側の正三角形の一辺の長さを2a，一番大きい正方形，二番目，三番目の正方形の一辺の長さをそれぞれ2b,2c,2dとし，以下の連立方程式を解く。2d=7.8179が与えられているが，dも未知数として方程式を解く。include("julia-source.txt");usingSymPy@symsa::positive,b::positive,c::positive,d::positived=7.8179/2sqrt3=sqrt(Sym(3))eq1=sqrt...算額（その303）
算額（その302）「三重県に現存する算額の研究」福島完(2007/2/13)https://mie-u.repo.nii.ac.jp/?action=repository_uri&item_id=7216三重県伊賀市菅原神社嘉永7年(1854)問題文5面積が163350歩平方の鈎股弦がある。鈎股弦の長さの差は分からないが，弦が82.5尺であるとき，鈎，股，中勾，方面，小面を求めよ。注:1尺=10歩正方形の一辺の長さをaとして，以下の連立方程式を解く。include("julia-source.txt");usingSymPy@syms鈎::positive,股::positive,弦::positive,a::positiveeq1=鈎*股/2-1633.50eq2=鈎^2+股^2-82.5^2eq3=(...算額（その302）
算額（その301）「三重県に現存する算額の研究」福島完(2007/2/13)https://mie-u.repo.nii.ac.jp/?action=repository_uri&item_id=7216三重県伊賀市永保寺天保15年正三角形の中に図のように円，正方形が入っている。正三角形の面積が2625平方寸のとき，それぞれの大きさを求めよ。正三角形の一辺の長さを2aとする。下の同心円の中心座標は(0,a/2)，大きい円の半径r1はa/2である。上にある円の半径をr2とし，以下の方程式を解き，a,r1,r2を求める。include("julia-source.txt");usingSymPy@symsa::positive,r1::positive,r2::positiveeq1=a*sqrt(Sym(3...算額（その301）
算額（その300）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（80)長野県長野市若槻東條雁田八幡宮文化3年(1806)文面一部欠損県内の算額（274)長野県長野市北尾張部連開庵奉納年不明二等辺三角形に全円が内接し，斜線で区切られた三角形にそれぞれ等円が内接している。全円の直径が3寸，中鉤（高さ）が4寸のとき，等円の直径を求めよ。5元連立方程式はSymPyのsolve()では解けないが，include("julia-source.txt");usingSymPy@symsh::positive,a::positive,w::positive,x0::positive,r0::positive,r1::positive,x1::positiv...算額（その300）
算額（その299）愛媛県伊佐爾波神社明治6年高阪金次郎峻則http://www.wasan.jp/ehime/isaniwa22.htmlhttps://www.city.matsuyama.ehime.jp/kanko/kankoguide/rekishibunka/bunkazai/ken/isaniwa_sangaku.htmlhttps://isaniwa.official.jp/2012/12/28/%E7%AE%97%E9%A1%8D%E3%81%AB%E6%8C%91%E6%88%A6/扇面（中心角120度）に東円1個，西円，南円，北円が2個ずつ入っている。南円の直径を知って北円の直径を求めよ。以下の方程式を解けばよいが，一挙に解くのが難しい。方程式は[eq4,eq7]と残りの方程式は(r0...算額（その299）
算額（その298）愛媛県伊佐爾波神社昭和12年(1937)算額22中村正教http://www.wasan.jp/ehime/isaniwa22.htmlhttps://isaniwa.official.jp/%E5%BE%A1%E5%AE%9D%E7%89%A9/%E7%AE%97%E9%A1%8D/大円の中に中円4個，小円1個が入っている。大円の径を知って，中円，小円の径を求めよ。大円，中円，小円の半径をr0,r1,r2とし，図に示すように記号を定め，以下の連立方程式をr1,r2について解く（r0は未知数として記号のまま残る）。include("julia-source.txt");usingSymPy@symsr1::positive,r2::positive,r0::positiveeq1=2r1...算額（その298）
算額（その297）中村信弥「算額への招待追補2」http://www.wasan.jp/syotai/syotai.html長野県千曲市八幡武水別神社（八幡神社）天保8年(1837)5個の正方形と1個の円が直角三角形に内接している。円の直径が7寸，一番小さい正方形の一辺の長さが6寸のとき，鉤，股，弦の長さを求めよ。(復元）算額の図が小さく（あるいは誤解？）分かりづらいが，円と一番小さい正方形の関係は，以下の拡大図のようになっている（算額では円が正方形に内接しているように描かれている）。また，「算額への招待追補2」は算額の術の解説をしており，上に示す1つの解のみを提示しているが，解は二通りある。もう一つの解の拡大図は以下のとおりである。記号の説明のために全体図も二通り示しておく。この2つの解は，図を裏返して...算額（その297）
算額（その296）中村信弥「算額への招待追補2」http://www.wasan.jp/syotai/syotai.html長野県千曲市八幡武水別神社（八幡神社）天保8年(1837)天円，地円，人円が互いに外接し，三角形に内接している。天円と地円の直径がそれぞれ8寸，4.5寸のとき，三角形の三辺の長さを求めよ。図に示すように記号を定め，以下の連立方程式をnlsolve()により解く。include("julia-source.txt");usingSymPy@symsr1::positive,r2::positive,r3::positive,x1::positive,x2::positive,x3::positive,x32::positive,a::positive,b::positive,c::po...算額（その296）
算額（その295）中村信弥「算額への招待追補2」http://www.wasan.jp/syotai/syotai.html長野県千曲市八幡武水別神社（八幡神社）天保8年(1837)直角三角形（鉤股弦）の中に，天円，地円，人円それぞれ1個，小円2個が入っている。天地人円の直径の和と大円の直径の積，および直角三角形の周の二乗の和が46216歩のとき，鉤および股を求めよ。図に示すように記号を定め，以下の連立方程式をnlsolve()により解く。include("julia-source.txt");usingSymPy@symsr1::positive,r2::positive,r3::positive,r4::positive,r5::positive,x2::positive,x3::positive,x...算額（その295）
算額（その294）中村信弥「算額への招待追補2」http://www.wasan.jp/syotai/syotai.html長野県千曲市八幡武水別神社（八幡神社）天保8年(1837)大円内に互いに外接し大円に内接する3個の小円がある。3個の小円の面積と，小円に囲まれる面積（水色）の和に小円の直径をかけた積に小円の直径の3倍の3乗を加えると794050立方寸になる。このとき，大円の直径を求めよ。大円の半径をr0，小円の半径，中心座標をr1,(r1,y),(0,r1)として，以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsr0::positive,r1::positive,y::negative;eq1=r1^2+(r0-r1-y)^2-4r1^2...算額（その294）
算額（その293）中村信弥「算額への招待追補1」http://www.wasan.jp/syotai/syotai.html長野県飯山市滝澤家所蔵http://www.wasan.jp/nagano/takizawa.html長方形内に3個の円が内接している。円の直径が5寸のとき，長方形の短辺の長さを求めよ。問では3個の円について述べており，右上，左上の2個の円についての記述はない。意図は良くわからないがここではこの円についても半径を求めておく。長方形の短辺の長さをyとする。半径，中心座標を以下のように決める。円r1,(r1,r1),(0,y-r1)小円r2,(r2,y-r2)以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsr1::positi...算額（その293）
算額（その292）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(282)長野県中野市片塩大徳寺天保9年(1838)大円と小円は互いに外接し，また2本の平行線の上下に接している。斜線は大円の共通接線である。大円と小円の直径がそれぞれ100寸，61寸のとき，黒く塗られた部分の面積（歩）を求めよ。原点と平行線の距離をyとする。小円の半径と中心座標r1,(0,r1+y)大円の半径と中心座標r2,(x2,r2+y)共通接線の傾きは±((r2+y)/x2)，共通接線と平行線の交点のx座標は±xである。include("julia-source.txt");usingSymPy@symsr1::positive,r2::positive,x2::posi...算額（その292）
算額（その291）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(281)長野県中野市片塩大徳寺天保9年(1838)上下の甲円は外接しそれぞれ左右の甲円にも外接している。乙円2個は互いに外接し，下の甲円と容円および左右の甲円に外接している。容円は上の甲円に内接している。甲円の直径が125寸のとき，容円の直径を求めよ。半径，中心座標を以下のように決める。甲円r1,(0,r1),(x1,0)乙円r2,(r2,y2)容円r3,(0,2r1-r3)r1=125/2であるが，r1も未知数として以下の連立方程式を解き(x1,r2,y2,r3)を求める（それぞれの式にr1が含まれる）。include("julia-source.txt");usingSym...算額（その291）
Juliaで2**(4+6)を入力すると，エラーなしで2*(4+6)を計算する。julia>2**(4+6)20Juliaの挙動不審
算額（その290）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(249)長野県諏訪市中洲諏訪大社上社明治12年(1897)長方形内に右下の頂点で長辺に接する円弧の一部があり，他に大円，中円，小円がある。それぞれは図のように互いに接している。中円の直径が500寸のとき，大円の直径を求めよ。長方形の長辺の長さをaとする。その一部が円弧である円の半径，中心座標をr0,(0,r0)とする。長方形の短辺の長さはsqrt(r0^2-a^2)である。大円の半径，中心座標をr1,(a-r1,r1)中円の半径，中心座標をr2,(x2,x2)小円の半径，中心座標をr3,(x3,y3)および(a-r3,r3)とする。以下の連立方程式をnlsolve()で解く。i...算額（その290）
算額（その289）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(250)長野県諏訪市中洲諏訪大社上社明治12年(1897)長方形内に斜線を描き，大円，小円が2個ずつ入っている。大円は長方形と斜線に接している。小円は長方形，大円，斜線に接している。長方形の短辺の長さが360寸，小円の直径が40寸のとき，長方形の長辺の長さを求めよ。長方形の長辺の長さをaとする。短辺の長さは大円の直径に等しい。大円の半径と中心座標をr1,(r1,r1)小円の半径と中心座標をr2,(x,2r1-r2)とする。斜線と長方形の長辺の交点の座標を(b,2r1),(a-b,0)とする。以下の方程式を解く。include("julia-source.txt");usingS...算額（その289）
算額（その288）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(250)長野県諏訪市中洲諏訪大社上社明治12年(1897)円内に2本の平行な弦と3個の正方形が入っている。正方形の一辺の長さが1寸のとき，円の径を求めよ。円の半径をRとする。正方形の一辺の長さをaとすれば，対角線の長さの半分はa/√2=√2a/2で，これをbとする。⊿OABにおいて，OB=R,OA=3b-R,AB=2b以下の方程式を解く。include("julia-source.txt");usingSymPy@symsR::positive,a::positive,b::positiveb=a/√Sym(2)eq=(3b-R)^2+(2b)^2-R^2solve(eq,R)...算額（その288）
算額（その287）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(251)長野県諏訪市中洲諏訪大社上社明治12年(1897)鉤股弦で，股×弦-鉤×弦=425平方寸弦-鉤=18寸のとき，鉤，股，弦を求めよ。以下の連立方程式を解く。(鉤^2+股^2)-弦^2を使わないのは，虚数解を除くため。include("julia-source.txt");usingSymPy@syms鉤::positive,股::positive,弦::positive;eq1=(股*弦-鉤*弦)-425eq2=(弦-鉤)-18eq3=(鉤^2+股^2)-弦^2eq3=sqrt(鉤^2+股^2)-弦res=solve([eq1,eq2,eq3],(鉤,股,弦));解は2組...算額（その287）
算額（その286）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(252)長野県諏訪市中洲諏訪大社上社明治12年(1897)等脚台形内に，甲，乙，丙の5個の円がある。それぞれ台形に内接したり，互いに外接しあっている。甲円の直径が1寸のとき，乙円の直径を求めよ。台形の上底と下底の長さを2b,2a，高さをcとする。甲円の半径と中心座標r1,(0,r1)乙円の半径と中心座標r2,(0,r2+2r1)上の丙円の半径と中心座標r3,(0,r3+2r1+2r2)右の丙円の半径と中心座標r3,(x3,r3)c=2(r1+r2+r3)である。以下の連立方程式を解く。なお，円の中心から直線までの平方距離は分数式になるので分子=0を解けば数式は若干短くなる。in...算額（その286）
算額（その285）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(253)長野県諏訪市中洲諏訪大社上社明治12年(1897)外円の中に弦を引き，等円を6個入れる。そのうち3個は弦の下，2個は弦の下に置き，それぞれが弦に接し，互いに外接する。また弦の下の3個の円の左右の等円は外円に内接する。残り1個は弦の上の2個に外接し，外円に内接する。等円の径が1寸のとき，外円の径を求めよ。外円，等円の半径をR,r，弦の位置をyとする（y<0,外円の中心から弦までの距離は-y）。等円の中心座標は図に示すようになる。rを変数のままにして，以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsR::posi...算額（その285）
算額（その284）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(247)長野県諏訪市中洲諏訪大社上社明治12年(1897)菱形の長，平（横の対角線と縦の対角線），面（一辺の長さ）の和に，面と長の再乗冪（3乗）を加えると901430720，菱形の面積は42504になる。このとき，長，平，面を求めよ。「面と長の再乗冪」は「面の再乗冪と長の再乗冪」ではない。連立方程式を解く。include("julia-source.txt");usingSymPy@syms長::positive,平::positive,面::positive;eq1=(長+平+面)+(面+長^3)-901430720eq2=長*平/2-42504eq3=(長^2+平^2)/4...算額（その284）
算額（その283）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(246)長野県諏訪市中洲諏訪大社上社明治12年(1897)大円内に，正三角形，元円，亨円，利円をそれぞれ3個，貞円を6個が入っているそれぞれは，大円，正三角形，隣り合う円に接している。大円の半径，中心座標R,(0,0)元円の半径，中心座標r1,(x1,r1)亨円の半径，中心座標r2,(x2,r2)利円の半径，中心座標r3,(0,R-r3)貞円の半径，中心座標r4,(x4,R-2r3+r4ここで，r2=16//2r3=(1//2-√3/4)Rx1=(R-r1)√3/2x2=(R-2r1-r2)√3/2である。図のように記号を定め，連立方程式を解く。include("julia-s...算額（その283）
算額（その282）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(246)長野県諏訪市中洲諏訪大社上社明治12年(1897)外円の中に菱形，その中に正方形がある。それらに等円が図のように接している。等円の径が8453寸のとき，外円の径を求めよ。正方形の一辺の長さを2a，等円の半径をr，外円の半径，中心座標をR,(0,0)とする。上の等円の半径，中心座標はr,(0,R-r)右の等円の半径，中心座標はr,(a+r,0)である。図のように記号を定め，連立方程式を解く。等円の半径rを未知数のまま解く。include("julia-source.txt");usingSymPy@symsR,r,ar=8453//2eq1=(R-2r)/R-a/(R-a...算額（その282）
算額（その279）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(239)長野県松本市筑摩筑摩神社明治6年(1873)甲円4個，乙円2個が互いに外接し，3個の円は直線に接している。乙円の径が1寸のとき，甲円の径を求めよ。甲円の半径をr1とする。右の甲円の中心は(√3r1,2r1)である。乙円の半径をr2とする。右の乙円の中心を(x2,r2)である。以下の連立方程式をr2を変数のままで解く。include("julia-source.txt");usingSymPy@symsr1::positive,x1::positive,y1::positive,r2::positive,x2::positive;r2=1//2y1=2r1x1=sqrt(...算額（その279）
算額（その281）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(240)長野県松本市筑摩筑摩神社明治6年(1873)等脚台形に対角線を引き，時計回りに90°回転させ上半分だけを考える。対角線で区切られた2つの区域に対角線及び辺に接する等円を入れる。元の台形の高さaを12寸，等円の径rを4寸とするとき，元の台形の底辺の半分bはいかほどか。図のように記号を定め，連立方程式を解く。include("julia-source.txt");usingSymPy@symsr::positive,a::positive,b::positive,c::positive,d::positive,y1::positive,dummy;(r,a)=(2,12)e...算額（その281）
算額（その280）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(240)長野県松本市筑摩筑摩神社明治6年(1873)鉤股弦（直角三角形）に2本の斜線を入れ，区切られた領域に2個の等円を入れる。鉤が5寸のとき，等円の直径を求めよ。図のように記号を定め，連立方程式を解く。include("julia-source.txt");usingSymPy@symsa::positive,b::positive,r::positive,x1::positive,y1::positive,x::positive,y::positive;eq1=distance(0,5,a,0,r,y1)-r^2eq2=distance(0,5,a,0,x1,r)-r^2e...算額（その280）
算額（その278）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(239)長野県松本市筑摩筑摩神社明治6年(1873)等脚台形に2本の斜線を引き，分割された区画に2個の等円が入っている。上底，下底がそれぞれ2寸4分，1尺2寸のとき，等円の径をもとめよ。小数点を避けるため，長さを「分」単位で表す。すなわち，上底，下底をそれぞれ24（分），120（分）とし，y軸を対称軸とする。等円の半径をr，台形の高さをy1,上底と下底の長さを2a,2b，斜線と脚の交点座標を(x2,y2)とし，図のように記号を定め，方程式を解く。include("julia-source.txt");usingSymPy@symsy1::positive,x2::positiv...算額（その278）
マイナ520万人分、総点検年金情報閲覧で政府>政府は13日、年金情報の閲覧問題で、国家公務員や私立学校教職員らの共済組合にもマイナンバー登録状況の総点検を求めると明らかにした。地方公務員の組合と合わせ、対象は約520万人になる。松本剛明総務相は「原因は人為的なミス」と説明し、制度そのものに問題はないと強調。野党側は混乱の拡大に批判を強めた。ヒューマンエラーは「ありうるもの」として，それを防ぐための防波堤を設計するのが普通。いままでの，「消えた年金」とか，「幼児の公金受け取り口座が親の名義の口座」とか，チェックすれば防げるのをやらないから問題が生じる。消えた年金だと，仮名書の姓名と実際の戸籍名を突合しないとか幼児のマイナンバーの名義と口座の名義が一致するかどうか確認しないとかむしろ，システム設計者やお役人は...事務処理のいろはのい
遅ればせながら，今日，市役所にマイナンバーカードを受け取りに行きました。そこで，実態に触れてみて，「これでいいのか?」と思う点が多々あったので記録しておきます。まず，私が窓口に行った段階で，待ち行列の人数は44人。待つこと1時間以上。窓口にたどり着いてから何をやったかというと。1.事前に書かされた証明用電子証明書の暗証番号を，受付のおねいさんが入力。2.同じく，利用者証明用電子証明証暗証番号を，受付のおねいさんが入力。3.その確認。これだと，1分もかからないと思うでしょう?いやいや，測っていないけど少なくとも20分はかかった。おまけに，「証明用電子証明書の暗証番号」は大文字アルファベットと数字の組み合わせで6〜16桁というのだけど，私の汚い文字をみて「これは数字ですか，アルファベットでか?」などというやり...個人番号カード...マイナンバーカードですわ...
算額（その277）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額(238)長野県松本市筑摩筑摩神社明治6年(1873)大小2つの正方形が接し，頂点を結ぶ斜線で区切られた区画内に大小の円が入っている。大円，小円の径がそれぞれ1尺8寸，8寸のとき，大きい方の正方形の一辺の長さを求めよ。大小の正方形の一辺の長さをa,bとおく。斜線が正方形と交わる座標を(a,c)とおく。include("julia-source.txt");usingSymPy@symsa::positive,b::positive,c::positive;eq1=a/(a+b)-c/b#三角形が相似であることからeq2=(a-c)/b-9//4#同じく相似関係からeq3=2a-c...算額（その277）
算額（その276）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（262）長野県上水内郡信州新町日名日置神社明治44年(1911)正三角形に円が内接している。三角面冪（正三角形の一辺の平方）が64平方寸，外積（正三角形と内接円の間の面積）が12.9平方寸のとき，内接円の直径を求めよ。三角率（正三角形の一辺の長さをaとしたときの正三角形の面積=a^2*(√3/4)における√3/4=0.4330127018922193）の近似値を0.433，円積率（直径がbのときの円の面積=(b/2)^2*π=b^2*(π/4)におけるπ/4=0.7853981633974483）の近似値を0.7とするとあるが，この近似値計算は基本的に間違っている。正三角形の一...算額（その276）
算額（その275）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（262）長野県飯綱町牟礼神社明治31年(1898)外球に円柱が内接している。甲球1個，乙球2個，丙球2個，丁球1個，戊球1個がそれぞれが隣の球（甲乙丙丁丙乙甲）の順に外接しており，外球に内接し円柱に外接している。また，外球を含めてすべての球の中心は同一平面上にある。外球の直径が85寸，甲球の直径が17寸のとき，乙球の直径を求めよ。問題は三次元で難しそうであるが，全ての球の中心が同一平面上にあるということから，上方からの透視図を考えると図のような二次元の問題であることがわかる。それぞれの円の半径，中心座標を以下のようにする。外円r0,(0,0)既知r0=85//2甲円r1,(0,...算額（その275）
算額（その274）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（165）長野県飯山市下木島鳥出神社弘化4年(1847)大円と小円は同心円で，半径の差は小円の半径に等しい。3個の小円の中心は大円の演習場にあり，小円は中遠に外接しかつ大円（中円）の弦にも接している。小円の径が4寸，矢が1寸のとき，中円の径はいかほどか。矢の長さをaとする。中円の半径と中心座標をr2,(0,0)とする。中央と右の小円の半径と中心座標をr1,(0,r1+r2),(x1,r2-a+r1)とする。r1,aを変数のまま以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsr1::positive,r2::positi...算額（その274）
今まで，ユニコード文字列のスライスに手こずっていたけど，Unicode.graphemes()でなんとかうまくいくようになった。多バイト文字とASCIIなどが混じっていても文字単位でスライスできる。In[23]:a="私の名前は中野です"Out[23]:"私の名前は中野です"In[24]:a[2:5]#これはエラーになるERROR:StringIndexError:invalidindex[2],validnearbyindices[1]=>'私',[4]=>'の'Stacktrace:[1]string_index_err(s::String,i::Int64)@Base./strings/string.jl:12[2]getindex(s::String,r::UnitRange{Int64})@Ba...Julia/v1.9.1での新機能の二つ目
新しい機能で，入力と出力に番号プレフィックスをつけることができるようになった。しかし，println()を使うとその機能が使えない。バグ？julia>usingREPL#まずはこれをやらないといけないのかな?julia>REPL.numbered_prompt!()#ついで，これ。In[2]:(a,b)=(2,45)Out[2]:(2,45)In[3]:println(a+2b)#println()の出力には[番号]がつかない92In[4]:a+2b#REPLだと，つく。Out[4]:92Julia/v1.9.1での新機能の一つ
Julia / Current stable release: v1.9.1 (June 7, 2023)
Juliaがv1.9.1へバージョンアップされましたCurrentstablerelease:v1.9.1(June7,2023)Julia/Currentstablerelease:v1.9.1(June7,2023)
算額（その273）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（165）長野県飯山市下木島鳥出神社天保14年(1843)大円内に中円1個，小円2個が入っている。外円，中円の径を知って小円の径を求めよ。大円の半径と中心座標をR,(0,0)とする。中円の半径と中心座標をr1,(0,R-r1)とする。右下の小円の半径と中心座標をr2,(r2,y2)とする。以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsR::positive,r1::positive,r2::positive,y2::negative;eq1=r2^2+y2^2-(R-r2)^2eq2=r2^2+(R-r1-y2)^2...算額（その273）
算額（その272）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（257）長野県下高井郡野沢温泉村湯沢神社明治17年(1884)円内に等円3個，三角形1個が入っている。等円の径を知って三角形の高さを求めよ。右下の等円の半径と中心座標をr,(r,R-3r)とする。外円の半径をRとすると，三角形の頂点Aの座標(x0,y0)で，y0=R-2r,x0=sqrt(R^2-(R-2r)^2)である。上の等円の中心から線分ABへの距離がrである。以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsR::positive,r::positive,x::positive,y::negative,x0::...算額（その272）
算額（その271）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（256）長野県下高井郡野沢温泉村湯沢神社明治17年(1884)直線の上に天円2個が互いに外接している。月円4個も直線に外接し両端の月園は天円にも外接している。人円2個は月円2個と，天円および日円に外接している。日円は天円とも外接している。天円の径を知って月円の径を求めよ。y軸が天円の交点を通るようにして，右半分だけを考える。天円の半径，中心座標をr1,(r1,r1)日円の半径，中心座標をr2,(0,y2)人円の半径，中心座標をr3,(x3,y3)月円の半径，中心座標をr4,(r4,r4)および(3r4,r4)天円と右端の月円が外接することから，天円の径が与えられれば月円の径は方...算額（その271）
算額（その270）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（260）長野県下高井郡木島平村天満宮推定明治21年(1888)直線の上に丁円，乙円，甲円，乙円の順に4円が外接している。この4円に外接する丙円がある。甲円，乙円の径がそれぞれ2寸，3寸のとき，丁円の径を求めよ。丁円の半径，中心座標をr0,(0,r0)甲円の半径，中心座標をr1,(x1,r1)乙円の半径，中心座標をr2,(x2,r2),(x3,r2),x3=x1+(x1-x3)=2x1-x2丙円の半径，中心座標をr4,(x1,2r1+r4)include("julia-source.txt");usingSymPy@symsr0::positive,r1::positive,r2...算額（その270）
算額（その269）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（260）長野県下高井郡木島平村天満宮推定明治21年(1888)大円2個，中円5個，小円4個が配置されている。小円は大円と中円に外接している。大円の直径が1尺（10寸）のとき，小円の直径を求めよ。術は「近似値」を求めているが，右上と左下の小円の中心間の距離を「大円の半径+小円の直径」とする無茶なものである。大円，中円の半径をr1,r2とおく（r2=r1/2）。小円の直径をr3,中心座標を(x3,x3)とおき，以下の方程式を解く。条件式が2本で，r3,x3の二変数について解を求める。その解の式にはr1が変数として含まれる。include("julia-source.txt");us...算額（その269）
算額（その268）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（183）長野県下高井郡木島平村西小路満昌院嘉永5年(1852)左下の頂点が原点の正方形と，2つの頂点がその正方形の辺上にあり，最初の正方形の右上の頂点が45度回転した正方形がある。左下の二等辺直角三角形中に甲円，右上の2個の二等辺直角三角形中にそれぞれ乙円が内接している。甲円の直径を乙円の直径で表わせ。甲円と乙円の入っている二等辺直角三角形は相似で，相似比が√2である。よって，「甲円の直径=√2乙円の直径」である。a=10.000000;b=6.666667;r1=1.952621;r2=1.380712;r1/r2=1.4142135623730951include("jul...算額（その268）
算額（その267）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（64）長野県下高井郡木島平村往郷水穂神社寛政12年(1800)甲，乙，丙の3個の正方形がある。甲，乙の面積の和をA，乙と丙の面積の和をBとする。このとき丙の一辺の長さを求めよ。ただし正方形の一辺の差は等しいとする。甲，乙，丙の3個の正方形の一辺の長さを，甲，乙，丙とする。include("julia-source.txt");usingSymPy@syms甲::positive,乙::positive,丙::positive,A::positive,B::positive;eq1=甲^2+乙^2-Aeq2=乙^2+丙^2-Beq3=(甲-乙)-(乙-丙)res=solve([e...算額（その267）
算額（その266）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額（171）長野県下高井郡木島平村上木島天然寺弘化2年(1845)図のように，正三角形の内側に3個の正方形がある。それぞれの正方形の2つの頂点は正三角形の辺の上にあり，他の2個の頂点は他の正方形の頂点と一致している。正三角形の一辺が8寸のとき，正方形の一辺の長さを求めよ。正三角形の一辺の長さをa，正方形の一辺の長さをxとすると，⊿ABCは一辺の長さがxの正三角形，∠OADは30°，OB=a/2ゆえ2(x+√3/2x)=a。xについて解くとx=(2-√3)a。aに8を代入するとx=2.14359353944898。include("julia-source.txt");usingSy...算額（その266）
算額（その265）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額3（91）長野県下高井郡木島平村穂高一川谷大元神社文化8年(1811)直角三角形（鉤股弦）内に円と正方形がある。三角形の面積は294平方寸，鉤の3乗と股の3上の和に円の直径（径）と正方形の一辺の長さ（方面）の差を掛けると62426となる（数字の掛け算で次元は無視）。このとき，円の直径を求めよ。include("julia-source.txt");usingSymPy@syms鉤::positive,股::positive,弦::positive,径::positive,方面::positive;eq1=(鉤+股-弦)-径eq2=(鉤*股/2)-294;eq3=(鉤^3+股^3)...算額（その265）
算額（その264）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額3（135）長野県長野市元善町善光寺天保3年(1832)正方形内に，甲，乙，丙，丁の4円が入っている。甲円か乙円の径が与えられているとき，正方形の一辺の長さを求めよ。正方形の一辺の長さをxとする。甲円の半径，中心座標をr1,(x-r1,x-r1)乙円の半径，中心座標をr2,(r2,r2)丙円の半径，中心座標をr3,(x-r3,r3)丁円の半径，中心座標をr4,(x4,r4)とする。図のように記号を定め，方程式を解く。甲円の径を与えるか乙円の径を与えるかは本質的なものではない。どちらにどのような値を与えても，得られる解（図）は相似である。つまり，各円の径及び正方形の一辺の長さは，ど...算額（その264）
算額（その263）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額2（94）長野県長野市三輪美和神社文化10年(1813)図のような，横4丈2尺，縦2丈4尺の長方形がある。長方形の内側に図のように糸を張ったところ，糸の長さが15丈になった。このとき，等弦の本数（折れ線の数）を求めよ。弦の本数をnとする。各弦は直角三角形の斜辺である。直角を挟む辺の長さは片方は24尺であり，もう一方は42/n尺であるので，弦の長さはsqrt((42/n)^2+24^2)である。これがn本あり，その合計が15丈=150尺なので，以下の方程式を解けばよい。include("julia-source.txt");usingSymPy@symsn::positive;eq...算額（その263）
算額（その262）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額2（71）長野県長野市久保寺観音堂享和3年(1803)鉤股弦（直角三角形）の中に円と正方形が入っている。方の一辺と円の直径の積が672平方寸，鉤股弦の周と中鉤の積が5644.8平方寸であるとき，方面（正方形の一辺）を求めよ。注：この2つの数値を，小数でなくするために10倍するのは，正しくない。以下の連立方程式を解く。include("julia-source.txt");usingSymPy@syms鉤::positive,股::positive,弦::positive,中鉤::positive,長弦::positive,短弦::positive,方面::positive,径::...算額（その262）
算額（その261）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額2（71）長野県長野市久保寺観音堂享和3年(1803)鉤股弦（直角三角形）の中鉤を隔てて，大円，小円が入っている。大円，小円の直径が68.8寸，51.6寸であるとき，鉤，股，弦を求めよ。大円，小円の中心座標を(x1,r1)，(股-r2,y2)とおき，以下の連立方程式を解く。include("julia-source.txt");usingSymPy@syms鉤::positive,股::positive,弦::positive,中鉤::positive,長弦::positive,短弦::positive,x1::positive,y2::positive;r1=688//2r2=...算額（その261）
算額（その260）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額2（70）長野県長野市久保寺観音堂享和3年(1803)圭（二等辺三角形）内に，大円1個，小円2個が入っている。中鉤（圭の高さ）が36寸，大円と小円2個の径の和が44寸のとき，大円の径を求めよ。大円と小円の半径をそれぞれr1,r2，圭の底辺の長さを2xとする。以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsl::positive,x::positive,r1::positive,r2::positive;l=sqrt(36^2+x^2)eq1=2r1+4r2-44#大円の径と小円の径の2倍（2個ぶん）の和が44寸eq2=...算額（その260）
算額（その259）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額2（70）長野県長野市久保寺観音堂享和3年(1803)鉤股弦（直角三角形）において，(1)鉤^3+股^3=223.894125(2)中鉤^3+短弦^3=48.361131のとき，鉤，股，弦を求めよ。「中鉤」，「短弦」は図の如し。以下の連立方程式を解く。include("julia-source.txt");usingSymPy@syms鉤::positive,股::positive,弦::positive,中鉤::positive,短弦::positive,長弦::positive;eq1=鉤^3+股^3-223.894125eq2=中鉤^3+短弦^3-48.361131eq3=...算額（その259）
算額（その258）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額2（70）長野県長野市久保寺観音堂享和3年(1803)菱形の中に大円1個と小円2個が入っている。菱形の辺の長さの和と菱形の横の対角線の長さの積が160平方寸，菱形の辺の長さの和と大円の直径の積が96平方寸のとき，小円の直径を求めよ。菱形の中心を原点とし，長い方と短い方の対角線の長さを2x,2yとおく。大円の半径，中心座標をr1,(0,0)，小円の半径，中心座標をr2,(r2,0)とおく。以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsr1::positive,r2::positive,x::positive,y::po...算額（その258）
算額（その257）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額2（112）長野県長野市信更町田野口清水神社文政11年(1828)正方形内に斜線で区切られた部分に円と正方形がある。円の径と内側の正方形の一辺の長さが等しい。外側の正方形の一辺の長さが100寸のとき，円の径を求めよ。外側の正方形の対角座標を(0,0),(x,x)とする。円の半径と中心座標をr,(x-r,x-r)内側の正方形の対角座標を(0,0),(r,r)とする。以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsx::positive,r::positive;x=100eq=2((x-r)-2r)^2-r^2res=so...算額（その257）
算額（その256）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額4（186)長野県長野市北尾張部連開庵嘉永6年(1853)正方形内に四分円，半円，天円，地円，人円が入っており，互いに内接・外接している。地円，人円の径がそれぞれ1寸4分4厘，1寸のとき天円の径を求めよ。四分円の半径をr0とする。半円の半径はr0/2である。天円の半径と中心座標をr1,(x1,r0-r1)地円の半径と中心座標をr2,(x2,r2)とする。r2=1.44である。人円の半径と中心座標をr3,(r3,y3)とする。r3=1である。以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsr0::positive,r1:...算額（その256）
算額（その255）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県外(285)東京都あきる野市二宮神社寛政6年(1794)勾股弦（直角三角形）に5個の等円（直径が同じ円）が弦（直角に対する変）及び鉤（直角を挟む短い方の辺），股（同じく長い方の辺）に接している。内接する全円の径が44寸のとき，等円の径を求めよ。全円の半径をr1=22寸とする。等円の半径をr2とする。以下の方程式をとき，まず勾，股，弦の長さを求める。include("julia-source.txt");usingSymPy#全円半径r1=22寸、#等円半径r2@syms鉤,股,弦,r1::positive,r2::positive,x::positive,y::positive弦=110...算額（その255）
算額（その249）中村信弥「改訂増補長野県の算額」http://www.wasan.jp/zoho/zoho.html県内の算額3（230)長野県諏訪市下諏訪諏訪大社下社慶応4年(1868)長方形内に甲，乙，丙，丁の円が1個ずつ入っている。丙円の径が5寸0分6厘2毛5糸，丁円の径が1寸4分4厘のとき，甲円の径はいかほどか。注：以下の説明図は算額の図に近くなるように，r3=2,r4=1としたものである。長方形の長辺，短辺の長さをそれぞれx，yとする。左下隅角を原点とする。甲円の半径，中心座標をr1,(x−r1,r1)とする。y=2r1である。乙円の半径，中心座標をr2,(r2,r2)とする。丙円の半径，中心座標をr3,(x3,y-r3)とする。r3=50625/2丁円の半径，中心座標をr4,(x4,r4)と...算額（その249）
「ブログリーダー」を活用して、r-de-rさんをフォローしませんか？
算額（その1694）六十二岩手県花泉町金沢大門神社・大門観世音菩薩慶応2年(1866)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html今有如図03052https://w.atwiki.jp/sangaku/pages/291.htmlキーワード：球4個，3次元#Julia,#SymPy,#算額,#和算長方形の中に4個の等半楕円と5個の等円を容れる。半楕円と2個の等円は1点で接している。楕円の短径が1寸のとき，長径はいかほどか。相変わらず，山村の図はいいかげんだ。右下の楕円の長半径，短半径，中心座標をa,b,(a,0)等円の半径と中心座標をr,(x,0),(0,b)共通接点の座標を(x0,y0)とおき，以...算額（その1694）
算額（その1694）四十七岩手県一関市平沢平沢白山神社慶応2年(1866)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html今有如図03056https://w.atwiki.jp/sangaku/pages/330.htmlキーワード：球4個，3次元#Julia,#SymPy,#算額,#和算盤上に甲球2個，乙球2個が互いに外接して載っている。高さ（盤面から乙球の天辺までの高さ）が69寸，甲球の直径が46寸のとき，乙球の直径はいかほどか。甲球の半径と中心座標をr1,(r1,0,r1),(-r1,0,r1)乙球の半径と中心座標をr2,(0,y2,r2),(0,y3,z2)高さをK=z2+r2とおき，以下の連立...算額（その1694）
算額（その1693）四岩手県花巻市北笹間東光寺慶応2年(1866)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html今有如図03054https://w.atwiki.jp/sangaku/pages/138.htmlキーワード：円，外円，正方形3個，斜線2本#Julia,#SymPy,#算額,#和算円の中に2本の斜線を隔てて，大正方形1個，小正方形2個を容れる。大小の正方形はそれぞれ円周上の2点で接し，大小正方形は一点で接する。明記されていないが2本の斜線は直角に交わる。すなわち斜線の2つの終端は円の直径上にある。山村の図はミスリーディングである。円の半径と中心座標をR,(0,0)大正方形，小正方形の一辺...算額（その1693）
「見たまま編集」のとき，テキスト中に半角の二重括弧「((」と「))」で囲まれる部分は脚注になる。ああだこうだの議論はさておき，解決法は間に空白を入れて二重括弧にしないこと。「((」にする。はてなブログを開設した段階で，何も設定しない場合，「見たまま編集」になっている。「見たまま編集」以外に設定した場合，はそれなりの別の問題が出てくる。一旦記事を作成（引っ越し）したら，別の編集方法に設定することはできない。以下のような場合，意図しないでもそうなってしまう。以下は数式(a+b)/((c+d)-e)再び文章また数式(x-y)*(a-(z-u))+2t再度文章以下のように表示される。人によっては重要：はてなブログへの移行時の問題：テキスト形式の数式の一部分が勝手に脚注になる不具合
引っ越し中です。なかなか時間がかかりますね。今後とも宜しくお願いします。https://sangaku0418.hatenablog.com/引っ越し中です
算額（その1692）新潟県三条市三条本成寺北越蒲原郡三十番神寛政12年(1800)神谷定令門人東都本郷松下清六郎興昌藤田嘉言(1807)：続神壁算法http://www.wasan.jp/jinpeki/jinpeki.html涌田和芳，外川一仁：三条本成寺の紛失算額，長岡工業高等専門学校研究紀要第49巻（2013）.https://kinpoku.nagaoka-ct.ac.jp/lib/kiyo/vol_46-50/vol_49/49_07wakuta2.pdfキーワード：菱形，等円#Julia,#SymPy,#算額,#和算菱形の中に2本の斜線を引き，等円を数個容れる（全部で偶数個）。菱形の対角線の長い方が4寸，短い方が3寸，等円の個数が全部で10個のとき，等円の直径はいかほどか。解法は涌田，外川（『...算額（その1692）
算額（その1691）は，5月1日に解禁します。算額（その1691）
算額（その1690）東都愛宕本地堂文化元年(1818)八木林平質門人總州參谷郷尾高卯之助頼之藤田貞資(1807)：続神壁算法http://www.wasan.jp/jinpeki/zokujinpekisanpo.pdfキーワード：円7個，楕円#Julia,#SymPy,#算額,#和算，#数学外円内に楕円1個と等円6個を容れる。上下の4個の等円は楕円と外接し外円に内接する，左右の2個の等円は楕円と1点で内接する（曲率円である）。外円の直径が19.9寸のとき，等円の直径はいかほどか。外円の半径と中心座標をR,(0,0)楕円の長半径，短半径と中心座標をa,b,(0,0);a=R等円の半径と中心座標をr,(r,y),(a-r,0)右上の等円と楕円の接点座標を(x0,y0)とおき，以下の連立方程式を解く。incl...算額（その1690）
算額（その1689）岐阜県郡上郡八幡町小野郡上八幡神社（小野天満宮）嘉永3年(1850)http://www.wasan.jp/gifu/gujyohatiman.html米山忠興：算額（その11）―郡上八幡神社―https://toyo.repo.nii.ac.jp/record/2533/files/shizenkagakuhen50_089-106_OCR.pdf上村文隆：はまぐりの数学https://hamaguri.sakura.ne.jp/mokuzi.html(153)美濃・飛騨の国の和算の歴史―算額の問題に挑戦しよう―https://hamaguri.sakura.ne.jp/minohidawasan.htmキーワード：回転楕円体3個，円錐台#Julia,#SymPy,#算額,#和算，#...算額（その1689）
算額（その1688）山形県新庄市戸沢神社文政元年(1818)http://www.wasan.jp/yamagata/tozawa.htmlキーワード：正方形3個，長方形#Julia,#SymPy,#算額,#和算，#数学長方形の中に大中小3個の正方形を容れる。長方形の長辺，短辺の長さが与えられたとき，小正方形の一辺の長さはいかほどか。長方形の長辺，短辺をa,b点Aの座標を(0,y2)点Bの座標を(x2,y1)点Cの座標を(x3,y6)点Dの座標を(x1,b)点Eの座標を(x5,0)点Fの座標を(a,y3)点Gの座標を(x4,y4)点Hの座標を(a,y5)点Iの座標を(x6,b)とおき，以下の連立方程式を解く。SymPyの能力上，解析解は求められないので，数値解を求める。include("julia-sou...算額（その1688）
算額（その1686）九十八江刺市雨宝堂現中善観音堂文政10年(1827)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html今有如図03006https://w.atwiki.jp/sangaku/pages/204.htmlキーワード：直角三角形，正三角形，正方形#Julia,#SymPy,#算額,#和算，#数学正方形の中に小さな正方形を2個容れ，甲円と乙円を2個ずつ容れる。乙円の直径が与えられたとき，甲円の直径はいかほどか。甲円の半径をr1乙円の半径をr2正方形の一辺の長さをaとおき，以下の連立方程式を解く。include("julia-source.txt");#julia-source.txtソースh...算額（その1686）
算額（その1685）七十一岩手県一関市川崎町薄衣諏訪前浪分神社(第4面)明治35年(1902)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html今有如図03088https://w.atwiki.jp/sangaku/pages/277.htmlキーワード：円9個，円弧2個#Julia,#SymPy,#算額,#和算，#数学2個の円弧が交わる中に，大円1個，中円4個，小円4個を容れる。小円の直径が与えられたとき，大円の直径を求める術を述べよ。円弧の半径と中心座標をR,(0,R-2r2)大円の半径と中心座標をr1,(0,0)中円の半径と中心座標をr2,(r1+r2,0)小円の半径と中心座標をr3,(x3,y3)...算額（その1685）
算額（その1684）六十四岩手県一関市花泉町金沢大門沢大門神社・大門観世音菩薩明治26年(1893)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html今有如図03077https://w.atwiki.jp/sangaku/pages/293.htmlキーワード：球1個，円錐3個，3次元#Julia,#SymPy,#算額,#和算，#数学盤上に3個の円錐が底面の中心が円周上にあり底面で互いに外接している。その中心に1個の球を容れる。球は盤面に接し，またそれぞれの円錐と1点で外接している。円錐の底面の直径と高さが与えられたとき，球の直径はいかほどか。算額（その1683）の類題である。いつもどおり，山村の図は不適...算額（その1684）
算額（その1683）四十六岩手県一関市真滝熊野白山滝神社明治31年(1898)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html今有如図03081https://w.atwiki.jp/sangaku/pages/329.htmlキーワード：球1個，円錐5個，3次元#Julia,#SymPy,#算額,#和算，#数学盤上に5個の円錐が底面の中心が円周上にあり底面で互いに外接している。その中心に1個の球を容れる。球は盤面に接し，またそれぞれの円錐と1点で外接している。球の高さと円錐の高さは等しい。円錐の底面の直径が与えられたとき，球の直径はいかほどか。いつもどおり，山村の図は不適切なものである。なぜか，円錐とも...算額（その1683）
算額（その1682）四十五岩手県一関市真滝熊野白山滝神社明治8年(1875)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html今有如図03061https://w.atwiki.jp/sangaku/pages/327.htmlキーワード：円5個，斜線2本，直線上#Julia,#SymPy,#算額,#和算，#数学直線上に全円1個，大円2個，小円2個が載っている。全円と大円，全円と小円は外接しており，大円と小円は共通接線を持つ。大円と小円の直径が与えられたとき，全円の直径はいかほどか。いつもどおり，山村の図は不適切なものである。全円の半径と中心座標をr1,(0,r1)大円の半径と中心座標をr2,(c,r2)小...算額（その1682）
算額（その1681）四十五岩手県一関市真滝熊野白山滝神社明治8年(1875)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html今有如図https://w.atwiki.jp/sangaku/pages/327.htmlキーワード：円1個，楕円3個，正三角形#Julia,#SymPy,#算額,#和算，#数学楕円3個が正三角形を囲み，正三角形に内接する円と楕円は外接している。楕円の短径が1寸のとき長径はいかほどか。山村の図はいつもどおり，ひどいものである。内接円の半径と中心座標をr,(0,0)楕円の長半径，短半径と中心座標をa,b,(0,3r/2)2個の楕円の接点座標を(x0,y0)とおき，以下の連立方程式を解...算額（その1681）
算額（その1680）七十岩手県一関市川崎町薄衣浪分神社文久年間山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円4個，楕円#Julia,#SymPy,#算額,#和算，#数学楕円と大円が交差しており，隙間に小円3個を容れる。楕円の長径，短径が与えられたとき，小円の直径はいかほどか。楕円の長半径，短半径，中心座標をa,b,(0,0)大円の半径と中心座標をr1,(0,r1-b)小円の半径と中心座標をr2,(d,0),(0,2r1-b-r2)とおき，以下の連立方程式を解く。include("julia-source.txt");#julia-source.txtソースhttps://blog.goo....算額（その1680）
5月1日に公開算額（その1679）
算額（その1678）宮城県塩竈市森山塩竈神社文化14年(1817)昭和60年復元奉納http://www.wasan.jp/miyagi/siogama4.htmlキーワード：円3個，楕円#Julia,#SymPy,#算額,#和算，#数学楕円の中に，甲円，乙円，丙円を容れる。楕円の長径が13寸，短径が5寸，乙円の直径が4寸のとき，丙円の直径はいかほどか。楕円の長半径，短半径，中心座標をa,b,(0,0)甲円の半径と中心座標をr1,(d1,0)乙円の半径と中心座標をr2,(d2,0)丙円の半径と中心座標をr3,(x3,y3)丙円と楕円の接点座標を(x0,y0)とおき，以下の連立方程式を解く。SymPyの能力的に，解析解を求めることができないので，数値解を求める。include("julia-source.t...算額（その1678）
算額（その1158）十七大里郡岡部村岡稲荷社文化13年(1816)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円10個，外円，直角三角形外円の中に直角三角形，甲円，乙円，丙円を1個ずつ，丁円，戊円，己円を2個ずつ容れる。丁円，戊円の直径がそれぞれ4寸3寸のとき，己円の直径はいかほどか。算額(その589）を一段階複雑にしたものである。直角三角形の斜辺は外円の直径である（中心を通る）。外円の半径と中心座標をR,(0,0)線分AB,CA,CBの長さをAB,CA,CBとする。甲円の半径と中心座標をr1,(x1,y1)乙円の半径と中心座標をr2,(0,r2-R)丙円の半径と中心座標をr3,(r3-R,0)丁円の半径と中心座標をr4,(x41,y41),(x4...算額（その1158）
算額（その1157）十七大里郡岡部村岡稲荷社文化13年(1816)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円11個，半円2個，長方形長方形の中に甲円，乙円，丙円の11個，半円2個を容れる。甲円の直径が5寸5分のとき，乙円の直径はいかほどか。長方形の長辺，短辺を2r1,r1甲円の半径と中心座標をr1,(0,0)乙円の半径と中心座標をr2,(r1+r2,0),(x2,r1-r2)丙円の半径と中心座標をr3,(x3,y3)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsr1::positive,r2::positive,x2::positive,r3::positive,x3::po...算額（その1157）
オクラの花ですハイビスカスと同じ仲間です2024/07/19
算額（その1156）九八鴻巣市三ツ木山王三ツ木神社明治28年(1895)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円7個，外円，半円1個，外円の中に，天円2個，地円2個，等円2個半を入れる。天円の直径が1寸のとき，地円の直径はいかほどか。外円の半径と中心座標をR,(0,0)天円の半径と中心座標をr1,(x1,4r2-R)地円の半径と中心座標をr3,(x3,2r2-R)等円の半径と中心座標をr2,(0,r2-R),(0,3r2-R),(0,5r2-R)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsR::positive,r1::positive,x1::positive,r2::...算額（その1156）
算額（その1155）九八鴻巣市三ツ木山王三ツ木神社明治28年(1895)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円4個，楕円，正方形，菱形正方形の中に楕円，菱形を1個ずつ，等円を4個容れる。楕円の長径，短径が4寸，3寸のとき，等円の直径はいかほどか。正方形の一辺の長さを2s楕円の長半径，短半径をa,b楕円と菱形の一辺との接点座標を(x0,y0)等円の半径と中心座標をr,(s-r,s-r)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symss::positive,r::positive,x0::positive,y0::positive,a::positive,b::positive...算額（その1155）
Juliaで統計解析2024/7/18にバージョンアップhttps://r-de-r.github.io/stats/Julia-stats1.htmlhttps://r-de-r.github.io/stats/Julia-stats2.htmlhttps://r-de-r.github.io/stats/Julia-stats3.htmlhttps://r-de-r.github.io/stats/Julia-stats4.htmlhttps://r-de-r.github.io/stats/Julia-stats5.htmlhttps://r-de-r.github.io/stats/Julia-stats6.htmlhttps://r-de-r.github.io/stats/Julia-stat...Juliaで統計解析2024/7/18版
算額（その1154）一〇一大宮市高鼻町氷川神社明治31年(1898)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円7個，楕円1個外円の中に楕円1個，甲円2個を容れ，隙間に乙円2個，丙円2個を容れる。乙円の直径が1寸のとき，丙円の直径はいかほどか。注：これだけでは条件不足である。「問」には書かれていないが，乙円は楕円の曲率円である。そのように解釈すれば条件は足りる。楕円の長半径，短半径と中心座標をa,b,(0,0);a=r1+r3,b=a-2r3外円の半径と中心座標をR,(0,0);R=a甲円の半径と中心座標をr1,(0,R-r1)乙円の半径と中心座標をr2,(R-r2,0)丙円の半径と中心座標をr3,(0,R-r3)とおき，以下の連立方程式を解く。...算額（その1154）
算額（その1153）九九春日部市小渕観音院明治30年(1897)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：正三角形，正方形，三角法面積が65996035.98832495593坪の正三角形の土地の一辺の長さはいかほどか。注：10間四方=1坪正三角形の一辺の長さがa，高さがhのとき，正三角形の面積Sは，S=(a/2)*(√3a/2)=a^2*√3/4である。√3/4=0.4330127018922193であるが，この近似値0.433が「三角法」と呼ばれていた。これは，辺の長さが同じ正三角形の面積(S)と正方形の面積(a^2)の比である。a=sqrt(S/(√3/4))≒sqrt(S/0.433)演算精度が必要なので，長精度演算を行う。正三角形の面積...算額（その1153）
算額（その1152）埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円7個，最大値面積が262155154.5032769612坪の正五角形の土地の一辺の長さはいかほどか。以下の図はGeoGebraで，一辺を与えた正多角形をかき，その面積を求めたものである。算額の「答」を得るために，「問」の数値を変えている。注：10間四方=1坪正五角形の一辺の長さをa，一辺の両端と正五角形の中心を結ぶ二等辺三角形の高さをhとする。h=a/2/tand(Sym(72)/2)で，正五角形の面積は，5*a*h/2=5*(a*(a/2)/tand(72/2))/2である。include("julia-source.txt")usingSymPy@symsS::positive...算額（その1152）
算額（その1151）七八加須市大字外野棘脱地蔵堂明治9年(1876)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円7個，最大値交差する2個の大円の交差部分に中円1個，小円2個を容れる。大円の直径が3.5寸のとき，小円の最大直径はいかほどか。注：交差する部分が小さければ小円も小さい。交差部分が大きくなるに連れ小円も大きくなるが，ある点をすぎれば逆に小さくなる。なお，図には大円の内部にもう一個（二個）の円があるがそれは交差部分の大きさに連動するので解を得るには無関係である。大円の半径と中心座標をr1,(0,r2)中円の半径と中心座標をr3,(0,0)小円の半径と中心座標をr4,(r3+r4,0)とおく。ただし，今回の方程式eq1は，大円，中円，小円の中...算額（その1151）
算額（その1150）五七羽生市小松小松神社安政6年()埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円2個，楕円，正六角形2個楕円の中に2個の等円と2個の正六角形を容れる。等円の直径が1寸のとき，正六角形の一辺の長さはいかほどか。正六角形の一辺の長さを2a楕円の長半径と短半径を4a,b等円の半径と中心座標をr,(0,b-r)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsa::positive,b::positive,r::positiveeq1=dist2(0,0,a,√Sym(3)a,0,b-r,r)eq2=(3a)^2/(4a)^2+(√Sym(3)a)^2/b^2-1res=...算額（その1150）
算額（その1149）番外九武州慈恩寺埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円2個，正三角形正方形と大円が直線上に載っている正方形の大円寄りにある頂点を通る大円の接線と，正方形の上辺を延長した直線と正方形の辺の3点で接する小円を描く。大円の直径が12寸，正方形の一辺の長さが10寸のとき，小円の直径はいかほどか。正方形の一辺の長さをa大円の半径と中心座標をr1,(r1,r1)小円の半径と中心座標をr2,(-a-r2,a-r2)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsa::positive,b::positive,r1::positive,r2::positiveeq1=d...算額（その1149）
算額（その1148）番外三川口市峯後峯が岡八幡宮（昭和10年代よりは前）埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円2個，直角三角形，正三角形直角三角形内に大円，小円，正三角形を容れる。鈎（直角を挟む二辺の短い方）が5寸のとき，図形の各寸法を得よ。大円の半径と中心座標をr1,(r1,y1)小円の半径と中心座標をr2,(r2,r2)正三角形の頂点の座標を(a,0),(b,0),((a+b)/2,√3(b-a)/2)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@syms鈎::positive,a::positive,b::positive,r1::positive,y1::positive...算額（その1148）
算額（その1147）一〇八加須市騎西町玉敷神社大正4年(1915)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：正三角形，正方形交差する大小の正三角形の中に3個の正方形を容れる。大きい正三角形の一辺の長さが10寸のとき，正方形の一辺の長さはいかほどか。大きい正三角形の一辺の長さを2a正方形の一辺の長さを2b下にある正方形の左右の頂点を通る，大きい正三角形の斜辺と並行な直線が作る正三角形は上にある正方形を内包する正三角形と合同である。下に新たに作った正三角形の底辺の長さを2cとおく。eq1はb,cのみを含む関係式，eq2はa,b,cを含む関係式である。両者を併せてb,cを求めれば，それぞれはaのみを含む式になる。include("julia-sourc...算額（その1147）
算額（その1146）一〇四北足立郡桶川町加納氷川天満神社明治43年(1910)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円2個，直角三角形，正方形直角三角形の中に正方形と2個の円を容れる。鈎，股が3寸，4寸のとき，正方形の一辺の長さと円の直径を求めよ。正方形の一辺の長さをa，円の直径をrとおき，以下の方程式を解く。include("julia-source.txt")usingSymPy@syms鈎::positive,股::positivr,a::positivr,r::positivreq1=(鈎-a)/a-鈎/股eq2=2a-sqrt(2a^2)-2r;まず，正方形の一辺の長さを求める。a=股*鈎/(股+鈎)で，鈎，股が3寸,4寸のとき，正...算額（その1146）
算額（その1145）一〇一大宮市高鼻町氷川神社明治31年(1898)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円2個，楕円2個，長方形，斜線2本長方形の中に2本の斜線を引き，大円2個と楕円2個を容れる。楕円の長径が4寸，大円の直径が8寸のとき，楕円の短径はいかほどか。長方形の長辺の長さを2c,（短辺の長さは2r）斜線と長方形の長辺との交点の座標を(d,0)大円の半径と中心座標をr,(c-r,r)楕円の長半径，短半径，中心座標をa,b,(0,b),(0,2b+a)楕円と斜線の接点座標を(x1,y1),(x2,y2)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsa::positiv...算額（その1145）
算額（その1144）一〇一大宮市高鼻町氷川神社明治31年(1898)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円5個，正三角形2個，正方形，対角線正方形の中に対角線，大円2個，中円1個，小円2個，正三角形2個を容れる。小円の直径が1寸のとき，大円の直径はいかほどか。図形を反時計回りに45°回転させると考えやすい。正方形の対角線の長さを2a大円の半径と中心座標をr1,(x1,-r1)中円の半径と中心座標をr2,(0,-r2-b)小円の半径と中心座標をr3,(x3,r3)小さな正三角形の高さをb=a*(3-√3)/2とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsa::positiv...算額（その1144）
算額（その1143）二九戸田市美女木八幡社天保4年(1833)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：直角三角形，正方形,繰り返し直角三角形の中に甲，乙，丙，丁の正方形を4個容れる。甲と丁の正方形の一辺の長さが与えられたとき，丙の正方形の一辺の長さはいかほどか。直角三角形の直角を挟む二辺の短い方を鈎，長い方を股とする。甲，乙，丙，丁の正方形の一辺の長さをa,b,c,dとおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@syms鈎::positive,股::positive,a::positive,b::positive,c::positive,d::positive,e::positive...算額（その1143）
算額（その1142）二五武州妻沼聖天宮文政11年(1828)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：累円，外円大円の中に等円3個を容れ，甲円から始まり，乙円，丙円，...，己円，...を容れる。等円の直径が3.5寸のとき，己円の直径はいかほどか。等円の半径と中心座標をr,(0,r),(0,0),(0,-r)大円の半径と中心座標を3r,(0,0)中円の半径と中心座標を2r,(0,r)甲円の半径と中心座標をr1,(x1,y1)乙円の半径と中心座標をr2,(x2,y2):とおき以下の連立方程式を解く。まず，甲円について，r1,x1,y1を求める。include("julia-source.txt")usingSymPy@symsr::positive...算額（その1142）
算額（その1141）二五武州妻沼聖天宮文政11年(1828)埼玉県立図書館：埼玉県史料集第二集『埼玉の算額』，昭和44年，誠美堂印刷所，埼玉県与野市.キーワード：円5個，外円，弦外円の中に水平な弦を引き，大円，中円，小円を容れる。大円，中円の直径がそれぞれ10寸，2.5寸のとき，小円の直径はいかほどか。弦とy軸の交点座標を(0,y);y=2r2-R外円の半径と中心座標をR,(0,0)大円の半径と中心座標をr1,(x1,y+r1)中円の半径と中心座標をr2,(x2,y+r2),(0,y-r2)小円の半径と中心座標をr3,(0,y+r3)とおき以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsR::positive,y::positive,r1::po...算額（その1141）
