MATLAB: Solving the Ordinary Differential Equation

ordinary differential equation

I am not sure how to solve these systems of differential equation. However, the final graph representation of the result is two exponential curves for and in respect to time.
Also, with =, the variable ks and BP are all constant.

Best Answer

EDITED
use dsolve()
or
Alternate method using ode45:
Screen Shot 2018-11-15 at 11.17.17 AM.png
tspan=[0 1];
y0=[0;0];
 
[t,x]=ode45(@myod,tspan,y0)
plot(t,x)
lgd=legend('Cp(t)','Cr(t)')
lgd.FontSize=20
 
function dxdt=myod(t,x) 
tau=2;
ks=3;
BP=6;
k1=5;
k2=7;
x(1)=exp(-t)/tau; %x(1)->Cp
dxdt=zeros(2,1); 
dxdt(1)=k1*x(1)-(k2/(1+BP))*x(2); %x(2)->Cr
dxdt(2)=k1*x(1)-k2*x(2);
end
Related Question