四十一 群馬県高崎市下小鳥町 幸宮神社 文政7年(1824) 群馬県和算研究会:群馬の算額，上武印刷株式会社，高崎市，1987年3月31日. 幸宮神社の算額 https://gunmawasan.web.fc2.com/files/sangak-corner/satimiya1.htm 東京都府中市 大国魂神社 明治18年(1885)…

大阪府茨木市 総持寺 嘉永7年(1854) 近畿数学史学会：近畿の算額「数学の絵馬を訪ねて」，平成4年5月16日 初版第一刷，大阪教育図書株式会社，大阪市. 奈良県大和郡山市小泉町 耳成山口神社 嘉永7年(1854) 近畿数学史学会：近畿の算額「数学の絵馬を訪ねて」…

三十 群馬県高崎市石原町 清水寺 文化12年(1815) 群馬県和算研究会:群馬の算額，上武印刷株式会社，高崎市，1987年3月31日.キーワード：累円#Julia #SymPy #算額 #和算 #数学 一点で接する半径の等しい弧がある。その弧にはさまれて，それぞれの弧に接し，ま…

macOS Ventura バージョン 13.0 (22A380) になったよ！！！(10月25日） iPhoneをWebカメラ代わりに使える「連係カメラ」 複数のアプリやウィンドウをまとめてグループ化する「ステージマネージャ」 メールアプリでは、一度送信操作を行なった直後の取り消し…

Python 3.11.0 がダウンロード可能になった Python.org からダウンロードしよう > Python を使って開発や実験を行うときは、用途に応じて専用の実行環境を作成し、切り替えて使用するのが一般的です。 などというフェークニュースには惑わされないで... まあ…

運用は自己責任で R : Nightly builds Mac は https://mac.r-project.org/ Windows は https://cran.r-project.org/bin/windows/base/R-devel-win.exe

「算額（その0002）」で，-2*r1^(3//2)*r2^(3//2)/(r1^2 - 2*r1*r2 + r2^2) + r1*r2*(r1 + r2)/(r1 - r2)^2 が単純には simplify できないと書いた。 以下のようにすれば，simlify できた。 a = -2*r1^(3//2)*r2^(3//2)/(r1^2 - 2*r1*r2 + r2^2) + r1*r2*(r1…

新しい Windows マシンを使えるようにしているのだがいろいろトラブルが発生する。 その一つに、R の立ち上げにはなんの問題もないのに、RStudio の立ち上げ時に以下の警告メッセージが出るというのに気づいた。オープニングの最初に出るので、狭い R コンソ…

四十一 群馬県高崎市下小鳥町 幸宮神社 文政7年(1824) 群馬県和算研究会:群馬の算額，上武印刷株式会社，高崎市，1987年3月31日. 算額に見る江戸時代の幾何学 ― 算額の問題に挑戦してみませんか？ ― 問題2 https://www.nikkei-science.com/page/magazine/980…

割当問題 Hungarian method n 人の人に n 個の仕事を割り当てるとき，最も効率のよい（コストの低い）割当方法を探る。 Project Euler の Problemm 345 で使用する。 Julia では Hungarian パッケージが用意されている。 using Hungarian 普通に使用すると，…

【C言語】フィボナッチ数列の10個の和と7番目の数との関係 https://sironekolab.com/archives/3021 1. 命題 フィボナッチ数列において，連続する 10 項の和は，第 7 項の 11 倍である。 2. 確認 フィボナッチ数列の生成関数 Fibonacci() を以下のように定義…

「青チャート式数学II」基本例題141 https://existence-scholar.com/modules/d3diary/index.php?page=detail&bid=2439&req_uid=1 極座標系でのグラフで悩んでいるようであるが，直交座標系で考えるべきである。（その後解決したようである）https://existenc…

ランキング参加中数学 ランキング参加中数学・科学・工学 ランキング参加中高校数学の基本 ランキング参加中ぼっちねこのパーフェクト数学教室☆*:.｡. o(≧▽≦)o .｡.:*☆ ランキング参加中【公式】2025年開設ブログ 「大学入試問題を SymPy で解く」を，R の car…

まずは，SymPy の summation() でやってみる。 using SymPy@syms n, keq = summation(n^2, (n, 1, k)) |> factor string(eq) # 二乗の和の公式 "k*(k + 1)*(2*k + 1)/6" eq(k => 3) 14 eq(k => 3000000000000000000) 90000000000000000045000000000000000005…

「Julia/SymPy: 数値代入法による係数の決定」であるが，一番簡単に解くには以下のようにする。 using SymPy@syms a, b, c, d, xeq = 9/(x-1)^2/(x+2)^2 - a/(x-1) - b/(x-1)^2 - c/(x+2) - d/(x+2)^2solve(eq, a, b, c, d) これで，解が以下のように示され…

三重県三重郡三重郡菰野町 廣幡神社 文化9年(1812) http://www.wasan.jp/mie/hirohata2.html 山形県鶴尾市大山 椙尾神社 文政元年(1818) http://www.wasan.jp/yamagata/sugio.html 一〇一 埼玉県大宮市高鼻町 氷川神社 明治31年(1898) 埼玉県立図書館：埼玉…

整式を整式で割ったときの商と剰余 using SymPy@syms x (x,) 1. 商 quo() quo(f, g, *gens, **args) Compute polynomial quotient of f and g. quo はエクスポートされていないので，sympy.quo() として使用する 例 x^2+1 を 2x - 4 で割ったときの商を求め…

三点を通る円の方程式を求める 以下を参照した 3点を通る円「ますいしいからの挑戦状(16)!!」をwolframAlphaとsympyでやってみた。https://qiita.com/mrrclb48z/items/62d3502529d09b0ef843 注: Julia の SymPy では，Point(), Circle(), Point2D() ともにエ…

Julia Version 1.8.2 (2022-09-29) が公開されていました

UTF-8 で，ASCII 文字と漢字かな混じりの文字列から，m番目からn番目までの文字列を取り出すとき，Julia では直接的な関数（指定法）がないのかなあ？ index = collect(eachindex(s)) で何文字目か？という情報を取り出し，それに基づいて文字列を抽出する。…

0^0 はいくつ？ インターネット上にも多くのページがある。 多くのプログラミング言語では 0^0 = 1 とされている。Julia も同じである。0^0 1 図を描いてみよう。 using Plotsx = 0:0.001:2plot(x, x.^x, label="", size=(300, 200)) 正の方向から 0 に近づ…

Julia で二重根号を外す 二重根号を外すには `sympy.sqrtdenest()` を使う。 using SymPyx = √(5 + 2√Sym(6))x |> N3.1462643699419723423291350657155704455124771291873287012324867174426654953709 y = sympy.sqrtdenest(√(Sym(5) + 2√Sym(6)))y |> N3.14…