« 遅延カーマック・マッケンドリック方程式をルンゲクッタ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版)をようやく見てきた。 »

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

コメント

コメントを書く

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

トラックバック


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

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

最近の記事

最近のコメント

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