忍者ブログ
20080511~ 13と7と11の倍数の論理積は13と7と11の積の倍数である。 和ァ・・・
[4485] [4484] [4483] [4482] [4481] [4480] [4479] [4478] [4477] [4476] [4475]
水素原子様の波動関数の一部である、球面調和関数の
緯度方向に関する微分方程式


を解くと、ルジャンドルの陪多項式が得られます。


これを、数値的に解いてみることにします。


まず、上述の微分方程式を差分方程式に変えます。
zの関数P(z)のzによる1階微分は

このように色々定義できますが、2階微分との絡みがあるため、今回は3つ目を採用しようと思います

また、2階微分のほうは

このように差分化されますね。


これを踏まえて微分方程式を差分方程式として解いていってみましょう

元の式はこのようになり

今回はP1とP0を初期値として与えてP2を求めるようにしたいので、この式をP2の式に変形します。





今日はおねむなんでこの辺で。
発散項どうにかならんかなぁ
そもそもの微分方程式の両辺に1-z^2をかけちゃったらどうなるやろか

拍手[0回]

PR

コメント


コメントフォーム
お名前
タイトル
文字色
メールアドレス
URL
コメント
パスワード
  Vodafone絵文字 i-mode絵文字 Ezweb絵文字


忍者ブログ [PR]
カレンダー
03 2024/04 05
S M T W T F S
1 2 3 4 5 6
7 8 9 10 11 12 13
14 15 16 17 18 19 20
21 22 23 24 25 26 27
28 29 30
ブログランキング
ブログランキング参戦中
にほんブログ村 アニメブログ 深夜アニメへ
にほんブログ村 漫画ブログ SF・ファンタジー漫画へ
にほんブログ村 科学ブログ 自然科学へ
よかったらポチッとお願いします^^
最新CM
[12/30 buy steroids credit card]
[09/26 Rositawok]
[03/24 hydraTep]
[03/18 Thomaniveigo]
[03/17 Robertaverm]
最新TB
プロフィール
HN:
量子きのこ
年齢:
43
性別:
男性
誕生日:
1981/04/04
職業:
WinDOS.N臣T
趣味:
妄想・計算・測定・アニメ
自己紹介:
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます
例:H→7番目、XP→15(P)×26+23(X)=413番目。
A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいズコー
バーコード
ブログ内検索
アクセス解析