解非线性方程组问题！！！请各位指教啊

基于MATLAB

高手们，我刚接触MATLAB有些可能会问一些很简单的问题
这是一个弹簧摆的极坐标数学模型
function dy=zhendong2(t,x)
r0=1;
m=0.1;
k=100;
g=9.8;
f1=x(3);
f2=x(4);
f3=-(k/m)*(x(1)-r0)+g*cos(x(2))+x(1)*x(4)^2;
f4=(-g*sin(x(2))-2*x(3)*x(4))/x(1);
dy=[f1;f2;f3;f4];
我用ODE45求解
>> tstop=50;
>> y0=[1;5];
>> [t,y]=ode45('zhendong2',tstop,y0,[]);
??? Index exceeds matrix dimensions.

Error in ==> F:\ProgramFiles\matlab\work\zhendong2.m
On line 6  ==> f1=x(3);

Error in ==> F:\ProgramFiles\matlab\toolbox\matlab\funfun\private\odearguments.m
On line 104  ==> f0 = feval(ode,t0,y0,args{:});

Error in ==> F:\ProgramFiles\matlab\toolbox\matlab\funfun\ode45.m
On line 155  ==> [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, args, ...
出现这些错误是什么原因啊

xionic

时间范围应该i是区间，你只给出了终止时间怎么行呢