[問題] 聯立微分方程式求數值解

[theye (Momiji)]

各位好，
因為不曾寫過聯立微分方程式，所以帶不出結果來 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=??)

請問我的寫法錯誤在哪呢？










--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 112.105.125.179
推 Frobenius:Ft是Ft:=Fa+Fb+Fc=10+0+0=10還是Ft:=Fa[V]+Fb[V]+Fc[V]?
推 Frobenius:k跟K一樣嗎?
推 Frobenius:其他式中的Fa,Fb,Fc是Fa[V],Fb[V],Fc[V]還是常數?
推 Frobenius:DSolve[...,{Fa, Fb, Fc},V]應改為Fa[V],Fb[V],Fc[V]
→ theye:謝謝費伯尼斯!
→ theye:藉費伯尼斯大大的建議，我的方程式已是可運算的狀態
→ theye:但這幾天還是試不出:給定最終條件，其他函數的值
→ theye:如給定Vf:=500，找 Fa[Vf]=??
→ theye:看課本都寫得很簡單的樣子  哈哈