大家帮忙看下程序

chenxf86542583

function Y=chenxf(x)
%dx/dt=a*x-b*x^2-m*x^2/(n+x^2)-h(t);受免疫项监视的Logistic模型;
%h(t)是周期性治疗外力
%m*x^2/(n+x^2)为免疫项
N=3360; %计算步数
L=0.05; %计算步长
a=1.0;%肿瘤细胞线性增长率
b=0.10;%环境营养限制
m=1.50;%m为免疫因子
n=2000;%ε是免疫系统激活阈值
t0=0;%初始时间
x0=50;%初始肿瘤细胞数量
h0=0;%治疗强度
h=100;
for k=1:N;
    %求解微分方程
k0=a*x-b*x^2-m*x^2/(n+x^2)-100;
k1=a*(x+L*k0/2)-b*(x+L*k0/2)^2-m*(x+L*k0/2)^2/[n+(x+L*k0/2)^2]-100;
k2=a*(x+L*k1/2)-b*(x+L*k1/2)^2-m*(x+L*k1/2)^2/[n+(x+L*k1/2)^2]-100;
k3=a*(x+L*k2/2)-b*(x+L*k2/2)^2-m*(x+L*k2/2)^2/[n+(x+L*k2/2)^2]-100;
x=x+(k0+2*k1+2*k2+k3)*L/6;
T=t0+k*L; %输出t为演化时间
H=100;
Y=x;     %肿瘤细胞数密度
plot3([t0 T],[h0 H],[x0 Y]);
end

运行后  ans=NaN，
程序有什么问题么？

ChaChing

每个人专业不同, LZ给什麽x?
建議樓主說清楚些, 並看下本版規則!