和算の1458次方程式をPython、SciPyの非線形ソルバーoptimize.rootで計算...しようと思ったがそのままではたぶん無理なのでa,b,c,d,e,fの6変数のまま計算する。a=10.000005716191469となった。他の変数も切りがいいので、たぶん答えから式を立てたんじゃないだろうか。
今日、この話を見た。
おお、これは以前、tsujimotterさんの動画でも見て、一回SciPyで計算していた。
ただその時は、定数が間違っていたのと、5つしか式がないままに計算していた。実は6つ目があったということで今回やり直し。
プログラムはめちゃくちゃ短い。初期値の選び方が大事ですが、それはのちほど。
答えはこれ
この文献の結果ともだいたい合ってる。
関孝和の問題を解く
古今算法記遺題の数値解について
多分答えからa=10, b=9, c=8, d=7.6666…, e=5, f=4と選んでから式を立てたんじゃなかろうか。
初期値の選び方ですが、aが決まればあと全部決まるので、aをいくつか振って6番目の式が最小になるように決めました。
こんな感じ。
するとa=10くらいで最小になるので、それを使って初期値を選んでいる。
func(sol.x)で検算しても
[0.0,
0.0,
-7.105427357601002e-14,
-5.684341886080802e-14,
1.4210854715202004e-14,
8.731149137020111e-11]
なかなかの精度。
« 数学フィールドワーク＆乱数（ちくま学芸文庫）を買った。数学フィールドワークはゾウリムシ、原子力、二酸化炭素、文学、カオスなど話題が豊富。乱数は学生時代買おうとして買わなかった本の復刻版。乱数の検定ってどうやるんだ？というのが気になっていた。 | トップページ | パンサー尾形さんのNHK笑わない数学、バーチ・スウィンナートン=ダイアー予想（BSD予想）を観てメモ。GeoGebraでいろいろお絵描きしたり、EDSACIIを調べたり（真空管コンピュータ！）、予想の論文見てみたりした。自分でもBSDの計算追いたくなった。 »
「パソコン・インターネット」カテゴリの記事
- 高周波回路シミュレータQucsStudioがuSimmicsに名称変更し、バージョンも4.8.3から5.8にアップデートされた。Qucsと区別するためだそうだ。また、Pythonの高周波用ライブラリscikit-rfもv1.5.0にバージョンアップされていた(2024.12.04)
- MATLAB Onlineで高周波基板設計用のRF PCB Toolboxを使ってみる。Coupled line バンドパスフィルタやratraceカプラが設計できる。モーメント法（MoM）や有限要素法(FEM)でちゃんと計算してくれているようだ。(2024.12.06)
- MATLAB Onlineで高周波用のRF Toolboxを使ってみる。Touchstoneファイルの読み込み、dB表示グラフ、スミスチャートなど簡単にできるし、フィルタ合成やIEEE P370 De-embedding(ZC-2xThru)も使える（MATLABで書かれたものがオリジナル）。(2024.12.05)
- MATLAB OnlineのSimulinkでローレンツ方程式をode8で計算してみる。Interface 2025年1月号でMATLAB Onlineの半年ライセンスがついてきたので。Simulinkを使うのは初めてだったが、わかりやすいSimulink入門コースを修了したのですぐできた。(2024.12.04)
- Interface2025年1月号はMATLABで1ニューロンから手作り 数学＆図解でディープ・ラーニング。初歩からAlexNetの転移学習、CNNまで話題が豊富で、なんとMatlab Onlineの半年ライセンスがついてくる。Simulinkや各種toolboxも使える。早速MATLAB入門オンラインコース修了した。(2024.12.03)
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2024年12月6日 NGMNが無線パフォーマンス評価フレームワーク発行、5GAAがC-V2Xのロードマップ発行、Marvellの3nm 1.6Tbps PAM4インターコネクト、Nokiaの2.4Tbps光伝送、Silicon Labsの低消費電力モジュール、Xiaomi 14T Pro分解動画(2024.12.06)
- 高周波回路シミュレータQucsStudioがuSimmicsに名称変更し、バージョンも4.8.3から5.8にアップデートされた。Qucsと区別するためだそうだ。また、Pythonの高周波用ライブラリscikit-rfもv1.5.0にバージョンアップされていた(2024.12.04)
- 日経サイエンス2025年1月号の特集 和算再発見の佐藤賢一さんの記事「算聖 関孝和の実像」に出てきた矢高に対する円弧の２乗の近似式をカシオの高精度計算サイトkeisan.casio.jpの自作式として作った。ものすごい精度であることがよくわかる。(2024.12.03)
- MATLAB Onlineで高周波基板設計用のRF PCB Toolboxを使ってみる。Coupled line バンドパスフィルタやratraceカプラが設計できる。モーメント法（MoM）や有限要素法(FEM)でちゃんと計算してくれているようだ。(2024.12.06)
- MATLAB Onlineで高周波用のRF Toolboxを使ってみる。Touchstoneファイルの読み込み、dB表示グラフ、スミスチャートなど簡単にできるし、フィルタ合成やIEEE P370 De-embedding(ZC-2xThru)も使える（MATLABで書かれたものがオリジナル）。(2024.12.05)
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2024年12月6日 NGMNが無線パフォーマンス評価フレームワーク発行、5GAAがC-V2Xのロードマップ発行、Marvellの3nm 1.6Tbps PAM4インターコネクト、Nokiaの2.4Tbps光伝送、Silicon Labsの低消費電力モジュール、Xiaomi 14T Pro分解動画(2024.12.06)
- 高周波回路シミュレータQucsStudioがuSimmicsに名称変更し、バージョンも4.8.3から5.8にアップデートされた。Qucsと区別するためだそうだ。また、Pythonの高周波用ライブラリscikit-rfもv1.5.0にバージョンアップされていた(2024.12.04)
- 日経サイエンス2025年1月号の特集 和算再発見の佐藤賢一さんの記事「算聖 関孝和の実像」に出てきた矢高に対する円弧の２乗の近似式をカシオの高精度計算サイトkeisan.casio.jpの自作式として作った。ものすごい精度であることがよくわかる。(2024.12.03)
- MATLAB Onlineで高周波基板設計用のRF PCB Toolboxを使ってみる。Coupled line バンドパスフィルタやratraceカプラが設計できる。モーメント法（MoM）や有限要素法(FEM)でちゃんと計算してくれているようだ。(2024.12.06)
- MATLAB Onlineで高周波用のRF Toolboxを使ってみる。Touchstoneファイルの読み込み、dB表示グラフ、スミスチャートなど簡単にできるし、フィルタ合成やIEEE P370 De-embedding(ZC-2xThru)も使える（MATLABで書かれたものがオリジナル）。(2024.12.05)
« 数学フィールドワーク＆乱数（ちくま学芸文庫）を買った。数学フィールドワークはゾウリムシ、原子力、二酸化炭素、文学、カオスなど話題が豊富。乱数は学生時代買おうとして買わなかった本の復刻版。乱数の検定ってどうやるんだ？というのが気になっていた。 | トップページ | パンサー尾形さんのNHK笑わない数学、バーチ・スウィンナートン=ダイアー予想（BSD予想）を観てメモ。GeoGebraでいろいろお絵描きしたり、EDSACIIを調べたり（真空管コンピュータ！）、予想の論文見てみたりした。自分でもBSDの計算追いたくなった。 »
コメント