ラプラス方程式をある差分化で数値計算すると高木関数になる例（F((2i+1)/^(2k+1))=0.5*(F(i/2^k)+F((i+1/2^k)) + Ck）

だいぶ前に書いたものの再掲：

昔、”無限・カオス・ゆらぎ”という本を読んで非常に面白かったのを覚えている。

その中でいちばんおもしろかったのは、差分化の方法によっては

元の微分方程式とは似ても似つかないカオス解が出ると言うところで、

その一つはロジスティック方程式の中心差分。

さらにもう一つ、”変なラプラシアン”というのが書いてあった。

（こちらの文献でも同様な話が読めます。

http://repository.kulib.kyoto-u.ac.jp/dspace/bitstream/2433/92310/1/KJ00004774034.pdf

）

ラプラス方程式は、ΔF(x)=0ですが、

一次元だと単にd² F / dx² =0。

単純な差分化だとF(x) = 0.5 * (F(x+Δx) + F(x-Δx)) ということになる。これをメッシュをちょっと工夫して

F((2i+1)/2^k+1)=0.5*(F(i/2^k)+F((i+1）/2^k)) + C_k

というように書きなおしてみる。

つまりあるkの点列が得られたとき、

次のk+1ではある点のデータは前のkの値で、

その点と隣り合ったデータを使って求めると見る（マルチグリッド）

でCk=1/2^(k+1)と選ぶ。Ckはk→∞では消えるはずの項。

（山口さんは、これはポアソン方程式じゃなくて変なラプラシアンだと言いきってますが。）

でkをだんだん増やして行った時のExcelのグラフをGIFアニメにしたのがこちら。

しかもこれはおなじくCkをk→∞で消える項1/4^kと選ぶと、

なだらかになるという。。。

単純な数値計算でも恐ろしいことになる例は多い。。。