% Drug absorption model - basic
k0 = 0.01;
k1 = 0.035;
[Tout, Cout] = ode45(@drugAbsorption, [0, 200], [200, 0, 0],[],k0, k1);
function Cout = drugAbsorption(t, ICs, k0, k1)
A = ICs(1);
B = ICs(2);
L = ICs(3);
dAdt = -k0*A;
dBdt = k0*A - k1*B;
dLdt = k1*B;
Cout = [dAdt; dBdt; dLdt];
end
如題，ode45 時定義function時，用A B作變數時為什么需要先設定A和B初始狀態，就是code里面的A = ICs(1);B = ICs(2)這兩部分，我看官方ode45里都沒有這一步，官方的例子在下面。
function dydt = vdp1(t,y)dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];[t,y] = ode45(@vdp1,[0 20],[2; 0]);第一次發帖，謝謝各位解答啦！
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/262655.html
標籤：其他開發語言