MATLAB Answers

Defining function file of ODE function

조회 수: 1(최근 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
  댓글 수: 1
Carey n'eville
Carey n'eville 2020년 11월 23일
Wow thank you so much!!

댓글을 달려면 로그인하십시오.

추가 답변(0개)

태그

제품


릴리스

R2020b

Community Treasure Hunt

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

Start Hunting!

Translated by