« インスタントラーメン発明記念館に行ってきた。 | トップページ | 一文橋のフジタカの看板もウインクする。 »

2009年6月16日 (火)

拡散方程式の解をkeisan.casio.jpにUP!

Wolfram|AlphaやMicrosoft Math3.0に浮気していたけれど、久しぶりに高精度計算サイトkeisan.casio.jpに自作式をUPした。これです↓

拡散方程式の解

一次元拡散方程式 ∂φ/∂t= D∂^2φ/∂x^2の解を計算します。
初期条件は高さA,幅εの矩形とします。

計算結果はこんな感じ。

18736000237

どうやって計算しているかというと、初期条件が矩形なら、誤差関数erf(x)を用いて

φ(x,t) = A/2 (erf((x+ε/2)/√(4Dt)) - erf((x-ε/2)/√(4Dt)) )

と厳密に解ける。ここで

erf(x) = (2/√π)∫[0,x] exp(-y^2) dy

です。誤差関数はkeisan.casio.jpにちゃんと入っているので簡単。

« インスタントラーメン発明記念館に行ってきた。 | トップページ | 一文橋のフジタカの看板もウインクする。 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック


この記事へのトラックバック一覧です: 拡散方程式の解をkeisan.casio.jpにUP!:

« インスタントラーメン発明記念館に行ってきた。 | トップページ | 一文橋のフジタカの看板もウインクする。 »

最近の記事

最近のコメント

2019年10月
    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 31    
フォト
無料ブログはココログ