如何用matlab解微分方程并作图

第一个是线性方程，可以直接求解析解
用dsolve
求数值解：
%%%%%%%%%%%%%%%%M-file
以下要单存一个文件在工作路径下%%%%%%%%%%%%%function
dotx=myfun(t,y)
qin=50；q=50;Cain=1.3;V=1000;K=0.05;
dotx=(qin*Cain-q*y)/V-K*y;
%%%%%%%%%%%%%%%%M-file
以上要单存一个文件在工作路径下%%%%%%%%%%%%%在command
窗口
[t,y]=ode45('myfun'，[0,t_final],y_int)%y_int
初值
第二个相同：
求数值解：
%%%%%%%%%%%%%%%%M-file
以下要单存一个文件在工作路径下%%%%%%%%%%%%%function
dotx=myfun2(t,y)
qin=50；q=50;Cain=1.3;Cbin=3;V=1000;K=0.05
dotx=[(qin*Cain-q*y(1))/V-K*y(1)*y(2);(qin*Cbin-q*y(2))/V-K*y(1)*y(2)];
%%%%%%%%%%%%%%%%M-file
以上要单存一个文件在工作路径下%%%%%%%%%%%%%在command
窗口
[t,y]=ode45('myfun2'，[0,t_final],[y1_int,y2_int])%yi_int是
初值