Defining function file of ODE function

%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));
 

Answer 1

[Stephan]

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

[Carey n'eville]

Wow thank you so much!!