解非线性方程组问题！！！请各位指教啊
高手们，我刚接触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=;
我用ODE45求解
>> tstop=50;
>> y0=;
>> =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, ...
出现这些错误是什么原因啊 时间范围应该i是区间，你只给出了终止时间怎么行呢
页:
[1]