Scratch(プログラム言語)で佐藤・テイト予想(sin^2θ予想)を確認。 #数楽
NHKの数学ミステリー白熱教室、ラングランズプログラム関連のScratchプロジェクトを作ってきています。
(その1)
(その2)
でさらに進んで、さっき書いた”海辺の美女の問題”で出てきた日本評論社の”数学100の問題”という本に出てきた佐藤・テイト予想をみてみようと思った。
黒木玄さんのTweetと難波完爾さんの資料も参考にしました。。
やったことは
η(τ)^2*η(5τ)^2=(1-x)^2*(1-x^2)^2…(1-x^5)^2*(1-x^10)...
を展開し、素数p=2n+1の場合の係数a_pから
x^2 - a_p*x+p=0の解を求め、その偏角θの分布が
sin^2(θ)になるかどうか見てます。
を展開し、素数p=2n+1の場合の係数a_pから
x^2 - a_p*x+p=0の解を求め、その偏角θの分布が
sin^2(θ)になるかどうか見てます。
計算したのは7000項(素数14000まで)です。難波さんの論文に合わせた。
もちろんScratchでも計算できるんですが、項が多いと時間かかるので別の方法(ExcelのVBA使った。一瞬。)で計算して、Scratchでその係数読みこませ、それをグラフ(猫のヒストグラム)にした。こういうの。
おお、sin^2θによく近似できてる。
« 海辺の美女の問題(結婚相手の最初の4割捨てろ理論?)をカシオの高精度計算サイトに作った(林先生が驚く初耳学?) #初耳学 | トップページ | 「モナドの領域」(筒井康隆)を読んだ。ザ・クレーター(手塚治虫)を思い出した。 »
「パソコン・インターネット」カテゴリの記事
- Google ColabのJulia言語でFPUT問題(Fermi–Pasta–Ulam–Tsingou、非線形結合した振動子が最初に与えたモードに戻る再帰現象)をDifferentialEquations.jlの2階の常微分方程式ソルバーDynamicalODEProblemでシンプレクティック8次のKahanLi8で計算、振動子の動きも動画にしてみる。(2025.05.22)
- PythonでFDTD法で電磁界シミュレーションできるopenEMSを使う(2)例題にあるマイクロストリップパッチアンテナ(MSA)を計算する。Sパラメータや入力インピーダンスだけでなく近傍界から遠方界の変換で指向性も計算できる。電流分布も動画で見る。給電はLumpedポートが使える。(2025.05.19)
- PythonでFDTD法で電磁界シミュレーションできるopenEMSを使う(1)例題にあるマイクロストリップラインのノッチフィルタ(スタブ)を動かして電磁界分布を動画で見てみる。CSXCADでモデルは確認できるし、ParaViewで電磁界分布が見られる。Sパラメータも計算できる。(2025.05.14)
« 海辺の美女の問題(結婚相手の最初の4割捨てろ理論?)をカシオの高精度計算サイトに作った(林先生が驚く初耳学?) #初耳学 | トップページ | 「モナドの領域」(筒井康隆)を読んだ。ザ・クレーター(手塚治虫)を思い出した。 »
コメント