下面是程序：
function dy=vdp4(t,y);
dy=zeros(10,1);
%定义：y(1)=Ta,y(2)=dTa/dx,y(3)=Tc,y(4)=dTc/dx,y(5)=Te,y(6)=dTe/dx,
%y(7)=Ts,j;y(8)=dTs,j/dx;y(9)=Ts,j+1;y(10)=dTs,j+1/dx;
%微分方程组的描述
dy(1)=y(2);
dy(2)={(6e4)*[y(1)-333]+(4e7)*[y(1)-y(7)]-(1.25e8)*[y(5)-y(1)]}*(3e-004)^2/1.6;
dy(3)=y(4);
dy(4)={(4.2e-003)*0.6^2+2.72*y(3)*0.6-200*0.6+2.76*y(3)+(7e4)*[y(3)-333]+(4e7)*[y(3)-y(9)]+(1.25e8)*[y(3)-y(5)]}*(3e-004)^2/1.6;
dy(5)=y(6);
dy(6)={(1.25e8)*[y(5)-y(1)]-(1.25e8)*[y(3)-y(5)]}*(1.3e-004)^2/0.43;
dy(7)=y(8);
dy(8)={(5e4)*[y(7)-333]+(5e6)*[y(7)-298]-(4e7)*[y(1)-y(7)]}*(4e-003)^2/188.1;
dy(9)=y(10);
dy(10)={(6e4)*[y(9)-333]+(4e6)*[y(9)-298]-(4e7)*[y(3)-y(9)]}*(4e-003)^2/188.1;

[t,y] = ode23t(@vdp4,0 ,1,[357 0 358 0 359 0 353 0 353 0]');
plot(t,y(:,1),'-',t,y(:,3),'-.',t,y(:,5),'.',t,y(:,7),':',t,y(:,9),'- -')

下面是出现的错误：

??? Error using ==> funfun\private\odearguments
When the first argument to ode23t is a function handle, the tspan argument must have at least two elements.

Error in ==> ode23t at 220
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...

Error in ==> vdp4main4 at 1
[t,y] = ode23t(@vdp4,0 ,1,[357 0 358 0 359 0 353 0 353 0]');

郁闷的是第二个错误，怎么会出现在ode23函数里？

请高手指教！谢谢！

function dy=vdp4(t,y);
dy=zeros(10,1);
%定义：y(1)=Ta,y(2)=dTa/dx,y(3)=Tc,y(4)=dTc/dx,y(5)=Te,y(6)=dTe/dx,
%y(7)=Ts,j;y(8)=dTs,j/dx;y(9)=Ts,j+1;y(10)=dTs,j+1/dx;
%微分方程组的描述
dy(1)=y(2);
dy(2)=((6e4)*[y(1)-333]+(4e7)*[y(1)-y(7)]-(1.25e8)*[y(5)-y(1)])*(3e-004)^2/1.6;
dy(3)=y(4);
dy(4)=((4.2e-003)*0.6^2+2.72*y(3)*0.6-200*0.6+2.76*y(3)+(7e4)*[y(3)-333]+(4e7)*[y(3)-y(9)]+(1.25e8)*[y(3)-y(5)])*(3e-004)^2/1.6;
dy(5)=y(6);
dy(6)=((1.25e8)*[y(5)-y(1)]-(1.25e8)*[y(3)-y(5)])*(1.3e-004)^2/0.43;
dy(7)=y(8);
dy(8)=((5e4)*[y(7)-333]+(5e6)*[y(7)-298]-(4e7)*[y(1)-y(7)])*(4e-003)^2/188.1;
dy(9)=y(10);
dy(10)=((6e4)*[y(9)-333]+(4e6)*[y(9)-298]-(4e7)*[y(3)-y(9)])*(4e-003)^2/188.1;

[t,y] = ode23t(@vdp4,[0 1],[357 0 358 0 359 0 353 0 353 0]');