算額，和算，数学，その他，Julia ときどき R, Python によるコンピュータプログラム，コンピュータ・サイエンス，統計学<br>
【R前処理講座12】{tibble}データフレームの進化版【tidyverse】https://datasciencemore.com/tibble/#tribbleを使用して，Markdownチックな表記でtibbleが作成できます．usingRCallR"""library(dplyr)trb=tribble(~No.,~language,#----- --- ----1,"R",2,"Python",3,"Julia")"""trbJuliaだと，usingCSV,DataFramesdata="""No.,language1,"R"2,"Python"3,"Julia""""df=CSV.read(IOBuffer(data),DataFrame)tribbleの効用?
base-R の forward pipe operator
pipeOp{base}RDocumentationBackgroundTheforwardpipeoperatorismotivatedbythepipeintroducedinthemagrittrpackage,butismorestreamlined.Itissimilartothepipeorpipelineoperatorsintroducedinotherlanguages,includingF#,Julia,andJavaScript.base-Rのforwardpipeoperator
Re. 【R前処理講座19】{dplyr} group_by：グルーピング【tidyverse】
【R前処理講座19】{dplyr}group_by：グルーピング【tidyverse】https://datasciencemore.com/dplyr-group-by/以下のようなデータフレームにおいて，グループごとの基礎統計量を求める。まずは1変数でグループ化する場合。usingRCallR"""library(dplyr)#definedataframedf=tibble(class=c("a","b","c","c","a","c","b","a","c","b","a","x"),gender=c("M","F","F","M","F","M","M","F","M","M","F","F"),height=c(162,150,168,173,162,198,182,154,175,160,172,...Re.【R前処理講座19】{dplyr}group_by：グルーピング【tidyverse】
Re: dplyrでカラム名や値を変数で指定してもうまく認識されない
dplyrでカラム名や値を変数で指定してもうまく認識されないhttps://trunk28.com/dplyr_for/引用元の記事は2021/07/11のものであるが，「dplyrでカラム名や値を変数で指定してもうまく認識されない」ということで，以下のプログラムはエラーになると...library(dplyr)hoge<-names(iris)[1]iris%>%select(hoge)エラー:Can'tsubsetcolumnsthatdon'texist.xColumn`hoge`doesn'texist.しかし，R4.2.0では，以下の"Note"が出るが，一応正しい結果が得られる。Note:Usinganexternalvectorinselectionsisambiguous.ℹUse`all_of...Re:dplyrでカラム名や値を変数で指定してもうまく認識されない
Juliaをもっと速くしよう！https://zenn.dev/ohno/articles/0ba7970d419898「多次元配列のメモリ上の位置はJulia,R,Fortranでは（Pythonと違い）列優先なので，badよりは，内側のforループでは左側の添字を変化させるgoodのほうが速い」と書かれている。筆者はWindowsのようであり，確かにgoodのほうが速いという結果が提示されている。私の環境はMacminiM1チップなのだが，両者の違いは殆どなかった。functionbad(arr)sum=0foriin1:size(arr)[1]forjin1:size(arr)[2]sum+=arr[i,j]endendreturnsumendfunctiongood(arr)sum=0forjin1:s...Re:Juliaをもっと速くしよう！
今更ながらではあるが，Rの>rep(1:3,each=3)[1]111222333>rep(1:3,3)[1]123123123Juliaではjulia>repeat([1,2,3],inner=3)9-elementVector{Int64}:111222333julia>repeat([1,2,3],outer=3)9-elementVector{Int64}:123123123julia>repeat([1,2,3],3)9-elementVector{Int64}:123123123ついでにjulia>repeat(1:3,inner=2,outer=3)18-elementVector{Int64}:112233112233112233Juliaの小ネタ--044Rのrep()
y=x^2とy=2^xの3つ目の共有点をニュートン法で探してみる
https://p-suugaku.blogspot.com/2022/03/yx2y2x3.htmlですが，煎じ詰めれば，x^2=2^xの解を求めるということ。元記事では第3の解をニュートン法を使って（しかもGoogleスプレッドシートを用いて）求めている。SymPyを使うと以下のようになる。ulia>usingSymPyjulia>@symsx(x,)julia>eq=Eq(x^2,2^x)2xx=2julia>a=solve(eq)3-elementVector{Sym}:24-2*LambertW(log(2)/2)/log(2)julia>a[3].evalf()-0.7666646959621233番目の解は，-2*LambertW(log(2)/2)/log(2)ということではあるが，a[3].e...y=x^2とy=2^xの3つ目の共有点をニュートン法で探してみる
Julia で，散布図と回帰直線の信頼限界帯を描画--その 2
色々あって，最終的にはが成果物として挙げられている。所要のプログラムは51行とのこと#Firstdrawabasicplotwithmainargumentsplot<-ggplot(data=data_iris,#setthedatamapping=aes(x=Sepal.Length,#setxaxisy=Petal.Length,#setyaxiscolour=Species)#setcolorbyspecies)+geom_point()+#setascatterplotgeom_smooth(method="lm")#setalinearline#Thenmakeitnicerplot<-plot+theme(axis.line=element_line(colour="black",#setcolo...Juliaで，散布図と回帰直線の信頼限界帯を描画--その2
ggplot2の入門ということで定番の図であるが，https://data-viz-lab.com/ggplot2にも以下の図がある。背景がネズミイロ，グリッド線がある，凡例が図の外にある，α=1なのだろう，データの重なり具合が分からない。書くのは簡単だよと>plot<-ggplot(data=data_iris,+mapping=aes(x=Sepal.Length,+y=Petal.Length))++geom_point(aes(colour=Species))++geom_smooth(method="lm",+colour="pink")>plotが挙げられているが，なんか知らないが，データフレームを専用に変換しないといけないとかで，事前にいかが必要とのこと。>#データセットのリストを作成>l<-da...Juliaで，散布図と回帰直線の信頼限界帯を描画
>だからCMSは嫌いだ。自分でhtmlやcgiを書いてWinSCPでアップロードする方が何がどうなっているかちゃんとわかるし、ずっと楽だ。なぜCMSがこんなに流行るのか理解できない。Rのグラフィクスでbaseが好きな自分にはggplot2がもてはやされるのが理解できないのだが、それと同じ感覚。たぶん自分がズレているのだが。いえいえ，いまの潮流が本筋からずれているだけだと思いましょう。ggplot2に限らず，tidyverseとかtidydataとかなんとか，わけも分からずグラフックの王道もわからずいいかんげんなグラフを再生産させ続けている御本尊にも相当な責任があると思います。新興宗教的なところがあると思います。世の中の流れが悪いだけです☺
Juliaに限らないのだけれどifmyrank==0a1=ones(n)elsea1=2*ones(n)endなどというのがあると，いやーな気分になるa1=ones(n)myrank!=0&&(a1*=2)とか書きたくなる。病気だなあ。プログラミングの常識?
Julia の MultivariateStats.fit(PCA, x)
method=:covの場合とmethod=:svdの場合の結果が違うという指摘に答えてくれて，MultivariateStatsv0.9.1へバージョンアップされた。それに伴い，以下を修正した。https://blog.goo.ne.jp/r-de-r/e/5e04aee16df93634f81a13d0dd3d5f45https://r-de-r.github.io/stats/Julia-stats8.htmlJuliaのMultivariateStats.fit(PCA,x)
なにはともあれ，実際に見てもらいましょう。ノートブック（Jupyterlabなど）に似ているけど，プログラミング中の変数値を変えると，関連する計算式の結果も変わる。ちょうど，Excelでセルの値を変えると他のセルの結果も変わるように。間違えたら間違った結果が出るので，どこをどう直したらどうなるかがすぐわかる。プログラミングの学習に役立つかな。また，一度作っておけば簡易電卓みたいに入力データだけ変えて結果を得るとかにも使えそう。コーディング部分は見えないようにもできる。JuliaのPluto（プログラミングの学習に最適？）
ケンドールの順位相関係数の信頼区間 kendall.ci() について
中澤先生が言及されていたNSM3:::kendall.ci()であるが，>bootstrap=TRUEと指定すれば（B=でシミュレーション回数）ブートストラップ推定もできるというすぐれものだが、>欠損値対応していないこと（欠損値を含むデータを与えるとエラーになる）、>信頼区間だけ表示されて点推定量が表示されないこと以外にも，プログラミング的にはダメダメなプログラムであった。最後に書き換えたプログラムを置いておくが，欠損値対応はok=complete.cases(x,y)x=x[ok]y=y[ok]の3行でできるし，点推定量は求めているのに表示していないだけなので，cat()で書き出すだけ。問題は，ブートストラップの部分。私の環境でexample=TRUE,B=60000でブートストラップをすると，31秒かかった...ケンドールの順位相関係数の信頼区間kendall.ci()について
Juliaでは，グラフ作成のためのパッケージが色々あるので，少しずつつまみ食いしている。GRMalie（Makie）を試してみたが，その中では以下のものが目に止まった。irisデータを使って書いてみた。ただし，Makieでは，gridを消すことができないそうだ。gridも含めてtickなどもまとめて消すことはできるようだが，それでは副作用が大きすぎる。また，標榜しているほどきれいでもない。なので，今後も使う予定はない。usingGLMakieusingFileIOusingColorTypesusingRDatasetsnoto_sans=assetpath("fonts","NotoSans-Regular.ttf")noto_sans_bold=assetpath("fonts","NotoSans-Bol...JuliaのMakie
「ブログリーダー」を活用して、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）
算額（その1128）四十四岩手県一関市真滝熊野白山滝神社文久元年(1861)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円5個，直線互いに接している大円と中円があり，両者を通る直線と各円の隙間に3個の小円を容れる。大円と小円の直径が10寸，1寸のとき，中円の直径はいかほどか。大円の半径と中心座標をr1,(x1,y1)中円の半径と中心座標をr2,(0,r2)小円の半径と中心座標をr3,(0,2r2-r3),(x3,2r2-3r3),(x1,y1+r1-r3)とおき，以下の連立方程式を解く。注：大円と小円の直径が10寸，1寸のとき，y1=0となるが，これは特殊な場合で常にy1=0が成り立つわけ...算額（その1128）
算額（その1127）四十四岩手県一関市熊野白山滝神社文久元年(1861)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円7個，外円全円の中に大円，中円，小円をそれぞれ2個ずつ容れる。大円，小円の直径がそれぞれ7寸，2寸のとき，中円の直径はいかほどか。山村の図は，意地悪としか思えないほど実際の姿と異なっている。山村の図が正しいものとして解こうとしても解は得られない。全円の半径と中心座標をR,(0,0);R=2r1大円の半径と中心座標をr1,(r1,0)中円の半径と中心座標をr2,(x2,y2)小円の半径と中心座標をr3,(0,y3)とおき，以下の連立方程式を解く。include("julia-...算額（その1127）
算額（その1126）四十一岩手県一関市牧沢牧沢八幡神社明治8年(1875)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円4個，二等辺三角形二等辺三角形の中に，甲円，乙円，丙円を容れる。甲円と丙円の直径がわかっているときに，乙円の直径を求めよ。与えられるものと求めるものが違うが，図形としては算額（その164）と同じである。二等辺三角形の底辺の長さを2x，高さをy甲円の半径と中心座標をr1,(0,r1)乙円の半径と中心座標をr2,(0,2r1+r2)丙円の半径と中心座標をr3,(x3,y3)とおき，以下の連立方程式を解く。一度に解くと有限の時間内に解が求まらないので，逐次解いていく。inclu...算額（その1126）
算額（その1125）三十四岩手県一関市舞川相川菅原神社後額嘉永3年(1850)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円1個，正三角形，正方形，斜線，黒積正方形の中に円，正三角形，斜線を容れる。円の直径が1寸のとき，図で示した灰色の部分の面積（黒積）を求めよ。正方形（正三角形）の一辺の長さをa円の半径と中心座標をr,(r,a-r)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsr::positive,a::positive,A::positive,B::positive,C::positiveA=2(a-√Sym(3)...算額（その1125）
算額（その1124）三十一岩手県一関市舞川観福寺内地蔵堂前額明治34年(1901)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：等球6個，円錐，3次元円錐の側面に6個の等級が接している。等球はの中心は全て円錐の底面に平行な水平面にあり，中は互いに接し合っている。円錐の底面の（円の）直径と高さが与えられているとき等球の直径はいかほどか。なお，「問」には書かれていないが（術には書かれている），等球の中心がある水平面と円錐の底面との距離も「答」に必要である。左上は，真横から見た図。6個の球の中心は同じ平面上にある。右上は，やや上から見た図左下は，円錐の頂点Hと底面の円の縁が重なるように見える位置か...算額（その1124）
算額（その1123）二十八一関市萩荘赤萩観音寺天保2年(1831)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円6個，長方形長方形の中に甲円1個，乙円1個，丙円2個，丁円2個を容れる。長方形の短辺の長さが1寸のとき，長辺の長さはいかほどか。「岩手の算額」も，解けるものは解いてきたが，解く前から「こりゃー，やってみても解けるかどうかわからない」というのが残っている。この算額の問題も，図を見ただけで，「なんだこりゃ」と後回しにしてきた。こんなの算額の図じゃない。案の定，以下のような図になるものであった。しかも，これもこの本に載っている問題の常で，答えが怪しい。算額の元々がおかしいのか，山村の解...算額（その1123）
算額（その1122）二十九一関市山ノ目配志和神社推定文政8年頃山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円8個，外円，二等辺三角形外円の中に二等辺三角形と甲円4個，乙円1個，丙円2個を容れる。丙円の直径が1寸のとき，乙円の直径はいかほどか。二等辺三角形の底辺とy軸の交点座標を(0,y)外円の半径と中心座標をR,(0,0)甲円の半径と中心座標をr1,(r1,y+r1),(r1,y-r1)乙円の半径と中心座標をr2,(0,y2)丙円の半径と中心座標をr3,(x3,y3)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsR::p...算額（その1122）
算額（その1120）十四前沢町赤生津長根前沢月山神社明治11年(1878)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：球3個，四角柱，3次元四角柱（底面が正方形）の中に等球を3個容れる。底面の正方形の一辺の長さが20寸，高さが30寸のとき，等球の直径を求めよ。正方形の一辺の長さをa，高さをh等球の半径と中心座標を下から順にr,(r,r,r),(a-r,a-r,z),(r,r,h-r);上下対称なのでz=h/2とおき，以下の方程式を解く。include("julia-source.txt")usingSymPy@symsa::positive,h::positive,r::positive,...算額（その1120）
算額（その1119）十岩手県胆沢町若柳市野々個人宅安政2年(1855)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円3個，正方形，斜線2本正方形の中に斜線2本を引き，甲円2個，乙円1個を容れる。乙円の直径が1寸のとき，短い方の斜線の長さはいかほどか。正方形の一辺の長さをa，斜線と正方形の斜線との交点を(0,b),(c,a)甲円の半径と中心座標をr1,(x1,r1)乙円の半径と中心座標をr2,(a-r2,y2)とおき，以下の連立方程式を解く。短い方の斜線の長さはsqrt((a-c)^2+a^2)である。include("julia-source.txt")usingSymPy@symsa::...算額（その1119）
算額（その1119）十岩手県胆沢町若柳市野々個人宅安政2年(1855)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円3個，正方形，斜線2本正方形の中に斜線2本を引き，甲円2個，乙円1個を容れる。乙円の直径が1寸のとき，短い方の斜線の長さはいかほどか。正方形の一辺の長さをa，斜線と正方形の斜線との交点を(0,b),(c,a)甲円の半径と中心座標をr1,(x1,r1)乙円の半径と中心座標をr2,(a-r2,y2)とおき，以下の連立方程式を解く。短い方の斜線の長さはsqrt((a-c)^2+a^2)である。include("julia-source.txt")usingSymPy@symsa::...算額（その1119）
算額（その1118）八岩手県水沢市南津田字化粧坂明治18年(1885)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円，直角三角形，正方形，矢，鈎，股，弦算額では基本的な三題である。1.直角三角形内の正方形鈎8.4寸，股14寸のとき，内接する正方形の一辺の長さ（方面）を求めよ。include("julia-source.txt")usingSymPy@syms鈎,股,方面eq1=(鈎-方面)/方面-鈎/股res1=solve(eq1,方面)[1]res1 >printlnres1(鈎=>8.4,股=>11.2).evalf() >println股*鈎/(股+鈎)4.8000000000000...算額（その1118）
算額（その1117）四岩手県花巻市南笹間東光寺慶応2年(1866)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：累円，外円大円の中に小円と初円を容れる。大円，小円と外接する円を初円と呼ぶ。初円が決まれば次の円（二円），また次の円（三円）と云うように，順次円（累円）が決まる。大円，小円，初円の直径がそれぞれ168寸，88寸，14寸のとき，次々に円を決めていくとき，円の直径が3寸になるのは初円を1番目として数えると何番目の円か。小円の中心がx軸上に来るようにして図を描く。大円の半径と中心座標をR,(0,0)小円の半径と中心座標をr1,(R-r1,0)初円の半径と中心座標をr2,(x2,y2)二円...算額（その1117）
算額（その1116）四岩手県花巻市南笹間東光寺慶応2年(1866)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：球，3次元盤の上に大球を載せ，その周りに小球を数個互いに接するように並べる。大球の直径と小球の個数が与えられたとき，小球の直径を求めよ。注：単に大球と書いているが，図にあるように，半球である大球の半径と中心座標をR,(0,0,0)小球の半径と中心座標をr,(x,0,r),(x*cos(2π/n),x*sin(2π/n)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsn::positiveinteger,R::pos...算額（その1116）
算額（その1115）四岩手県花巻市南笹間東光寺慶応2年(1866)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：球4個，円錐台，3次元円錐台の中に直径4寸の大球2個と直径2寸の小球2個が入っている。円錐台の高さを求めよ。注：大球，小球は円錐台の側面にそれぞれ2点で接し，小球は円錐台の上面に接している。2個の小球が重なって見える方向から見た図2個の大球が重なって見える方向から見た図円錐の高さ，底面の円の半径をh，a大球の半径と中心座標をr1,(r1,0,r1)小球の半径と中心座標をr2,(0,r2,z2)とおき，以下の連立方程式を解く。円錐台の高さは，小球の中心のz座標値に小球の半径を加えたも...算額（その1115）
算額（その1114）百二岩手県大船渡市猪川町気仙猪川村安養寺安政年間か山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円5個，楕円，正方形正方形の中に楕円1個と等円を5個容れる。等円の直径が与えられたとき，正方形の一辺の長さを求めよ。正方形の一辺の長さを2a楕円の長半径と短半径，中心座標をa,b,(0,b)等円の半径と中心座標をr,(0,2a-r),(r,b+r),(r,b-r)右上の等円と楕円の接点座標を(x0,y0)とおき，異関連列方程式を解く。include("julia-source.txt")usingSymPy@symsa::positice,b::positive,x0::pos...算額（その1114）
算額（その1113）百一岩手県大船渡市猪川町田茂山町神明杜田茂山神社（奉納年不明）山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円6個，円弧，正三角形円弧の中に正三角形1個，甲円2個，乙円4個を容れる。甲円の直径が1寸のとき，弦の長さを求めよ。正三角形の一辺の長さを2a円弧の半径と中心座標をR,(0,0)甲円の半径と中心座標をr1,(x1,y0+r1)乙円の半径と中心座標をr2,(x21,y0+r2),(x22,y22)弦とy軸の交点座標をy0とおき，以下の連立方程式を解く。弦の長さは2sqrt(R^2-y0^2)である。甲円の中心と正三角形の底辺の右側の頂点を結ぶ直線は円弧の中心を通り，x...算額（その1113）
算額（その1111）百大船渡市猪川町雨宝堂現雨宝山竜宝院文政10年(1827)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円5個，正三角形大円2個と正三角形が交差している。隙間に小円を3個容れる。小円の直径が1寸のとき，正三角形の一辺の長さはいかほどか。正三角形の一辺の長さを2a大円の半径と中心座標をr1,(x1,r1);x1=r1-r2小円の半径と中心座標をr2,(0,r1),(x2,y2);x2=r1とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsa::positice,r1::positive,r2::positiv...算額（その1111）
算額（その1111）百大船渡市猪川町雨宝堂現雨宝山竜宝院文政10年(1827)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円7個，外円大円の中に中円5個，小円1個を容れる。小円の直径が1寸のとき，大円の直径はいかほどか。山村の図には小円が描かれていないが，小円は上の図の位置にある。大円の半径と中心座標をR,(0,0)中円の半径と中心座標をr1,(0,R-r1),(r1,y1),(r1,y1-2r1)小円の半径と中心座標をr2,(0,r2-R)とおき，以下の連立方程式を解く。include("julia-source.txt")usingSymPy@symsR::positive,r1::p...算額（その1111）
算額（その1110）四十一岩手県一関市牧沢牧沢八幡神社明治8年(1875)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円5個，正方形，斜線2本長方形の中に，斜線を2本引き，甲円2個，乙円3個を容れる。乙円の直径が1寸のとき，平（直平;長方形の短辺）はいかほどか。下図に示すように，乙円が同じ大きさでも，長方形の短辺はいかようにもなる。つまり，解はない。「答」，「術」，山村ともに，この点には一切触れていない。要するに条件不足なので，追加条件として，たとえば長方形の長辺の長さを与えるようにすれば解は一意に求まる。長方形の長辺，短辺を2a,b斜線と長辺の交点座標を(c,0)甲円の半径と中心座標をr...算額（その1110）
算額（その1109）四十一岩手県一関市牧沢牧沢八幡神社明治8年(1875)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.htmlキーワード：円5個，外円，弦全円の中に弦を引き，その上下に4個の等円を容れる。等円の直径が1寸のとき，全円の直径はいかほどか。注：弦の傾斜をどのようにしようとも，弦の上下（両側）は対称，すなわち弦は全円の直径であり，弦の傾斜はどのようにしても同じ図形になる（図を回転すれば弦は水平にできる）。弦を水平にすれば，4個の等円はそれぞれがx軸，y軸に接する。全円の半径と中心座標をR,(0,0)第一象限にある等円の半径と中心座標をr,(r1,r1)とおき，以下の方程式を解く。include("j...算額（その1109）
