各位好,
因為不曾寫過聯立微分方程式,所以帶不出結果來 Q_Q
還請多指教。
我先是定義一些常數與四則運算式
Kc:= 0.05
Ft:= Fa+Fb+Fc
k:= 0.7
Ct0:= 0.2
ra:=-k*Ct0*(Fa/Ft - Ct0/K* Fb/Ft* Fc/Ft)
kc:= 0.2
接著有三條變數都是V的微分方程式,分別是:
d(Fa)/dV = ra
d(Fb)/dV = -ra - (kc*Ct0*Fb)/Ft
d(Fa)/dV = -ra
個別起始條件是 Fa(V=0) = 10, Fb(V=0) = Fc(V=0) = 0
而我欲想的寫法是:
DSovle[{Fa'[V]==ra, Fb'[V]==-ra-(kc*Ct0*Fb)/Ft, Fc'[V]==-ra,
Fa[0]==10, Fb[0]==0, Fc[0]==0}, {Fa, Fb, Fc}, V]
不過很明顯是跑不出結果 (V=??)
請問我的寫法錯誤在哪呢?