注意：在通常情况下，matlab的工作区中不会显示函数内的变量。

以如下的SEIR模型为例：

function SEIR(xN,xE,xI,xR,xr,xr2,xB,xB2,xa,xy)

%第一阶段fun1(837300,0,1,0,10,10,0.045,0.045,0.1,0.9)
N = xN;     %总人数                                                             
E = xE;        %潜伏人数                                                          
I = xI;          %感染人数                                                                                                          
R = xR;          %康复人数 
S = N - I-R;       %易感人数  

r = xr;        %一个感染者的接触人数  
r2 =  xr2 ;     %一个潜伏者的接触人数 
B = xB;%x1;%0.045;      %接触感染者，易感到潜伏的传染率                                          
B2 = xB2; %接触潜伏者，易感到潜伏的传染率
a = xa;      %潜伏到感染的传染率                                                   
y = xy;      %感染者的治愈率                                                          


T = 1:150;%天数
for i = 1:length(T)-1 %也就是1：149
   
     S(i+1) = S(i) - r*B*S(i)*I(i)/N-r2*B2*S(i)*E(i)/N ;
     E(i+1) = E(i) + r*B*S(i)*I(i)/N(1)+r2*B2*S(i)*E(i)/N(1)-a*E(i);
     I(i+1) = I(i) + a*E(i) - y*I(i);
     R(i+1) = R(i) + y*I(i);
     Y(i+1)=I(i+1)+I(i)+ R(i+1);
    
end


subplot(1,2,1);
plot(T,Y,'r');grid on;%绘制多条曲线并显示网格

hold on;
subplot(1,2,2);
plot(T,I,'r');grid on;%绘制多条曲线并显示网格
hold on;
end

执行代码发现工作区不显示数组Y,I,T的值： 

解决方法：

以函数输出的形式定义希望输出的变量：

function [A1,B1,C1]=SEIR(xN,xE,xI,xR,xr,xr2,xB,xB2,xa,xy)
% 
%   xr,xr2大于0     xB,xB2,xa,xy都属于[0，1]
%第一阶段fun1(837300,0,1,0,10,10,0.045,0.045,0.1,0.9)
N = xN;     %总人数                                                             
E = xE;        %潜伏人数                                                          
I = xI;          %感染人数                                                                                                          
R = xR;          %康复人数 
S = N - I-R;       %易感人数  

r = xr;        %一个感染者的接触人数  
r2 =  xr2 ;     %一个潜伏者的接触人数 
B = xB;%x1;%0.045;      %接触感染者，易感到潜伏的传染率                                          
B2 = xB2; %接触潜伏者，易感到潜伏的传染率
a = xa;      %潜伏到感染的传染率                                                   
y = xy;      %感染者的治愈率                                                          


T = 1:150;%天数
for i = 1:length(T)-1 %也就是1：149
   
     S(i+1) = S(i) - r*B*S(i)*I(i)/N-r2*B2*S(i)*E(i)/N ;
     E(i+1) = E(i) + r*B*S(i)*I(i)/N(1)+r2*B2*S(i)*E(i)/N(1)-a*E(i);
     I(i+1) = I(i) + a*E(i) - y*I(i);
     R(i+1) = R(i) + y*I(i);
     Y(i+1)=I(i+1)+I(i)+ R(i+1);
    
end

% 
%经典SEIR模型的差分方程（间隔为1天）
A1 = Y;
B1 = I;
C1 = T;
subplot(1,2,1);
plot(T,Y,'r');grid on;%绘制多条曲线并显示网格

hold on;
subplot(1,2,2);
plot(T,I,'r');grid on;%绘制多条曲线并显示网格
hold on;
end

改变的地方在这：

 参考文章