Maxima で綴る数学の旅

紙と鉛筆の代わりに、数式処理システムMaxima / Macsyma を使って、数学を楽しみましょう

数学

-数学- q級数を高速に展開するプログラム

桜の花びら 色々と恥ずかしいコードなのですが、とりあえずここに公開してみます。このコードは全て著者が書いたものです。ライセンスはとりあえずGNU GPL v2としました。商用で使いたいからBSDライセンスを希望する、などの方はご相談ください(いないと思…

-数学- q級数の計算をもっと高速化

お花見 q級数の展開計算をしたい、と思ったのです。以前に書いた記事 ではMaximaのプログラムで書いてみたわけですが、1000次まで正確に計算するのに8分もかかるという遅さでした。これでは色々と試すには遅すぎですよね。 そこで今回全部Lisp言語で書き直し…

-数学- 「7は合同数」をMaximaで確認してみよう

にほんブログ村 tsujimotterさんの以下の記事に触発されてこの記事を書きました。いつも興味深い記事を投稿されているtsujimotterさんには改めて多謝です。 すべての辺の長さが有理数であるような直角三角形の面積になる数、それが合同数です。どんな自然数…

-数学- 虚2次無理数をモジュラー変換で基本領域に移す

キティちゃんのチョコレート。 にほんブログ村 2次無理数αを根として持つ整数係数既約2次多項式方程式の判別式 D(α)(以下、略して2次無理数の判別式と呼びます)が保型性を持つ、というお話を紹介し、その証明をしてみました。 保型性といえば基本領域で…

-数学- 2次無理数の判別式の保型性の証明

ご機嫌です。 ある2次無理数 \( x \) に対して元の2次式の判別式の値を\( D\left(x\right) \)とします。今回は実数\( a,b,c,d \)が\( a\,d - b\,c=1 \)を満たすとき、 $$D\left(x\right)=D\left(\frac{b+a\,x}{c\,x+d}\right)$$ が成り立つことを示します。 …

-数学- 2次無理数の判別式の保型性

犬鍋 お久しぶりでございます。 あまりに久しぶりの記事です。数学も久しぶりで、リハビリを兼ねてこちらの本を読んでいます。 数論入門 (現代数学への入門) 作者: 山本芳彦 出版社/メーカー: 岩波書店 発売日: 2003/11/11 メディア: 単行本 購入: 1人 クリ…

-数学- NHK 数学ミステリー白熱教室「ラングランズ・プログラムへの招待」

NHKのEテレで 数学ミステリー白熱教室「ラングランズ・プログラムへの招待」 という番組が放送されるそうです。 4回シリーズで第1回は2015年11月13日の午後11時から。つまり今夜ですよ。皆さん、必ずご覧ください。 ラングランズ・プログラムといえば、フ…

-数学- ディリクレ級数の計算 係数が約数のa乗和とb乗和の積の場合

お昼寝 引き続き、黒川教授のこちらの本からです。 ラマヌジャンζの衝撃 (双書―大数学者の数学) 作者: 黒川信重 出版社/メーカー: 現代数学社 発売日: 2015/08 メディア: 単行本 この商品を含むブログを見る 整数nの約数のa乗和と同じnの約数のb乗和をかけた…

-数学- ディリクレ級数の計算 係数が約数のk乗和の場合

にほんブログ村 黒川教授がこの8月に出版された本: ラマヌジャンζの衝撃 (双書―大数学者の数学) 作者: 黒川信重 出版社/メーカー: 現代数学社 発売日: 2015/08 メディア: 単行本 この商品を含むブログを見る から、 \(\sum_{n=1}^{\infty }{\frac{\sigma_{…

-数学, Maxima入門- 数論的関数のディリクレ級数、texput関数の使い方

ディスクキャッチ! にほんブログ村 ラマヌジャンζの衝撃 (双書―大数学者の数学) 作者: 黒川信重 出版社/メーカー: 現代数学社 発売日: 2015/08 メディア: 単行本 この商品を含むブログを見る 黒川教授の上記の本を読んでいて、以下の式が出てきました。 $$ …

-数学- ラマヌジャン ζの衝撃 (ゼータの衝撃)

にほんブログ村 【楽天ブックスならいつでも送料無料】ラマヌジャンζの衝撃 [ 黒川信重 ] この本、購入しました。以前こちらの記事 -数学- アイゼンシュタイン級数関連の記事の参考文献 - Maxima で綴る数学の旅で雑誌「現代数学」に黒川教授がゼータ関数に…

-数学- ゼータ関数の非自明零点と行列の固有値(その4)今回使ったMaximaのパッケージについて

波飛沫に霞む にほんブログ村 <a href="http://maxima.hatenablog.jp/entry/2015/08/16/012243" data-mce-href="http://maxima.hatenablog.jp/entry/2015/08/16/012243">-数学- ゼータ関数の非自明零点と行列の固有値(その1)エルミート行列の固有値を求める - Maxima で綴る数学の旅</a>および(その2)ではグラフ描画のためのdrawパッケージ以外にもdistrib, descriptive, lapackというパッケージ…

- ゼータ関数の非自明零点と行列の固有値(その3)物語

お友達 にほんブログ村 エルミートランダム行列の固有値の間隔分布とゼータ関数の非自明零点の間隔分布、あまりに出自の異なる2つの概念が、グラフに書くと完全に一致している、、、。ということはちょうどゼータの非自明零点を固有値とするエルミート行列…

-数学- ゼータ関数の非自明零点と行列の固有値(その2)ゼータの非自明零点の間隔分布を求める

ボール取ってきました にほんブログ村 前の記事ではランダムエルミート行列の固有値を求め、その間隔分布をグラフにプロットし、理論分布との一致を確認したのでした。 <a href="http://maxima.hatenablog.jp/entry/2015/08/16/012243" data-mce-href="http://maxima.hatenablog.jp/entry/2015/08/16/012243">-数学- ゼータ関数の非自明零点と行列の固有値(その1)エルミート行列の固有値を求め</a>…

-数学- ゼータ関数の非自明零点と行列の固有値(その1)エルミート行列の固有値を求める

にほんブログ村 量子力学で登場する波動関数の記述に現れる演算子はエルミート行列なのですが、この挙動を調べるためにランダムエルミート行列というものが物理の方では研究されています。行列の性質を知るためには特性方程式やその根である固有値が重要にな…

-数学- n番目の素数を求める公式

プールで泳ぐパピヨン にほんブログ村 ご無沙汰しておりました。 常日頃フォローしているブログに以下のようなツィートが紹介されていました。 「素数をnで表す式は発見されてない」みたいなの見る度に別にイライラはしないがnで表す式はいくらでも見つかっ…

-数学- 不思議な対称性「2つの曲線は等しい」

にほんブログ村 $$ f\left(t\right):=\exp \left(i\,t\right)-\frac{\exp \left(6\,i\,t\right)}{2}+\frac{i\,\exp \left(\left(-14\right)\,i\,t\right)}{3} $$ この指数関数の3つの項からなる関数は複素平面上で閉曲線を描きます。一方、コメントに書かれ…

-数学- 「不思議な対称性」へのコメントについて

にほんブログ村 記事 <a href="http://maxima.hatenablog.jp/entry/2015/06/21/223648" data-mce-href="http://maxima.hatenablog.jp/entry/2015/06/21/223648">-数学- 不思議な対称性 - Maxima で綴る数学の旅</a>に膨大なサイズの式によるコメントをいただきました。コメントの式と記事で紹介した複素関数の式が同じ曲線を表していることを証明せよ、というコメントでした。 最初は失礼ながらスパム…

-数学- 不思議な対称性

にほんブログ村 Creating Symmetry: The Artful Mathematics of Wallpaper Patterns 作者: Frank A. Farris 出版社/メーカー: Princeton University Press 発売日: 2015/06/02 メディア: Kindle版 この商品を含むブログを見る とあるブログで上記の本の紹介…

-数学- 手触り感のある数学、さわれるゼータ関数

にほんブログ村 数学って抽象度が高くなると何をやっているのか分からなくなることがよくあります。そんな時は具体例を計算したり、図やグラフを書いたり、、、。そんな手触り感のある数学は楽しいですね。 このオブジェは手触り感のある数学の愛好家(専門…

-Maxima入門/数学- グラフ理論パッケージの応用 スペクトルグラフ理論

にほんブログ村 グラフ理論パッケージと行列関連の機能を使うと、色々とグラフの性質を計算できるそうです。特に行列の固有値から色々な性質がわかる理論があり、「スペクトルグラフ理論」と呼ばれているようです。 なんだか響きもカッコイイですよね。「ス…

-数学- q級数の計算の高速化

// にほんブログ村 詳しい(そして実に面白いのですが)ことは、tsujimotterさんのブログ記事: &lt;a href="http://tsujimotter.hatenablog.com/entry/6xx-xy-yy" data-mce-href="http://tsujimotter.hatenablog.com/entry/6xx-xy-yy"&gt;6xx + xy + yy の…

-数学- アイゼンシュタイン級数関連の記事の参考文献

// にほんブログ村 アイゼンシュタイン級数関連の記事を書いてきました。 もともとの動機は「現代数学」、 現代数学 2014年 10月号 [雑誌] 出版社/メーカー: 現代数学社 発売日: 2014/09/12 メディア: 雑誌 この商品を含むブログを見る という雑誌に東工大の…

-数学- モジュラー形式とラマヌジャンの不思議な等式

// メディテレーニアンハーバー にほんブログ村 アイゼンシュタイン級数の保型性を使うと、ラマヌジャンの不思議な等式(%o1)を示すことができます。 (%i1) sum(n^5/(exp(2*%pi*n)-1),n,1,inf)=1/504; $$ \tag{%o1} \sum_{n=1}^{\infty }{\frac{n^5}{e^{2\,\p…

-数学- リプシッツ公式を導く

// にほんブログ村 リプシッツ公式(%o1)はオイラーのコタンジェント公式から導くことができます。 (%i1) lip:sum((z-n)^(-k),n,-inf,inf)=(-2*%pi*%i)^k/(k-1)!*sum(n^(k-1)*q^n,n,1,inf); $$ \tag{%o1} \sum_{n=-\infty }^{\infty }{\frac{1}{\left(z-n\rig…

-数学- 約数のk乗和の母関数

// ラウンジ にほんブログ村 ランベルト級数の式変形と関係して約数のk乗和の母関数が登場しました。その辺の計算をしてみます。 ではまずMaximaでの表示をきれいにするための準備です。powerdispl:trueとすると級数を次数の昇冪順に表示します。またtexput(…

-数学- アイゼンシュタイン級数のフーリエ展開をリプシッツ公式から導く

// 疾走 にほんブログ村 重さkのアイゼンシュタイン級数のフーリエ展開の係数は約数のk乗の和が登場するのがとても面白いです。どうしてそんなことが起こるのでしょうか。 まず、約数のk乗和関数をimaximaやMaxima on Androidでカッコよく表示するためのおま…

-数学- アイゼンシュタイン級数のフーリエ展開

// 東京駅ライトアップ にほんブログ村 アイゼンシュタイン級数はいろいろな形で数論と結びついています。すでに見たように、m,nが(0,0)を除くすべての整数を走る、代わりにm,nが互いに素なすべての整数を走るように変形することができ、その時にはゼータ関…

-数学- アイゼンシュタイン級数がモジュラー形式であること

// 模擬レース にほんブログ村 いつものようにアイゼンシュタイン級数を(%i1)で定義します。 (%i1) Ez1:Ez[k](z)=lsum(1/(m*z+n)^k, mn,Z^2);$$ \tag{%o1} \mathrm{Ez}_{k}(z)=\sum_{\mathrm{mn}\in{Z^2}}{\frac{1}{\left(m\,z+n\right)^{k}}} $$(%i2) modu:…

-数学- アイゼンシュタイン級数 (その3)

// ウォームアップ にほんブログ村 アイゼンシュタイン級数にはz->z+1やz->-1/zなどの変数変換をしてもあまり変化しません。(これらの変換である種の対称性があります)。 アイゼンシュタイン級数の定義は以下の通りでした。ただし\( \left(m, \,n\right)=\…

This is a test.
天気