« 遅延カーマック・マッケンドリック方程式をルンゲクッタ8次でExcel VBAで。 | トップページ | 「アバター」(3D版)をようやく見てきた。 »

2010年1月22日 (金)

遅延微分方程式(酵素反応論)をルンゲクッタ8次でExcel VBAで。

次は酵素反応論からの例。方程式は

y1'(x) = I - z*y1(x)

y2'(x) = z*y1(x) - y2(x)

y3'(x) = y2(x) - y3(x)

y4'(x) = y3(x) - 0.5*y4(x)

ここで

z=1/(1+0.0005*y4(x-4)^3)

で遅延している。

計算結果はこちら。

Kouso002 

この文献の結果を再現してますね。だいぶ慣れてきた。次いってみよう。

« 遅延カーマック・マッケンドリック方程式をルンゲクッタ8次でExcel VBAで。 | トップページ | 「アバター」(3D版)をようやく見てきた。 »

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

コメント

コメントを書く

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

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/47106104

この記事へのトラックバック一覧です: 遅延微分方程式(酵素反応論)をルンゲクッタ8次でExcel VBAで。:

« 遅延カーマック・マッケンドリック方程式をルンゲクッタ8次でExcel VBAで。 | トップページ | 「アバター」(3D版)をようやく見てきた。 »

最近のコメント

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