Defining function file of ODE function

2 次查看（过去 30 天）
Carey n'eville
Carey n'eville 2020-11-23
评论： Carey n'eville 2020-11-23
%We have initial concentration C(0)=5, but I didn't understand how I use this value and also my code does'nt work.
function dCdt=Conc(C,t)
global t;
global C;
k1=0.7;
t=0:1:6;
Vo=1;
Q=Vo;
Cin=200;
dCdt=-(k1+(t/(Vo+Q*t)))*C+((Q*Cin)/(Vo+Q*t));
end
%I get this error, Could you help to me, please:
%Error using /
%Matrix dimensions must agree.
%Error in Conc (line 10)
% dCdt=-(k1+(t/(Vo+Q*t)))*C+mrdivide((Q*Cin),(Vo+Q*t));

请先登录，再回答此问题。

采纳的回答

Stephan
Stephan 2020-11-23
编辑：Stephan 2020-11-23
C0 = 5;
tspan = [0 10];
[t,C] = ode45(@Conc,tspan,C0);
plot(t,C)
function dCdt=Conc(t,C)
k1=0.7;
Vo=1;
Q=Vo;
Cin=200;
dCdt=-(k1+(t./(Vo+Q.*t)))*C+((Q*Cin)/(Vo+Q.*t));
end

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

产品


版本

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by