Defining function file of ODE function

2 views (last 30 days)
%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));

Accepted Answer

Stephan
Stephan on 23 Nov 2020
Edited: Stephan on 23 Nov 2020
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

More Answers (0)

Categories

Find more on Programming in Help Center and File Exchange

Products


Release

R2020b

Community Treasure Hunt

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

Start Hunting!