function z=rab(x,y)& o0 b. K. \- @6 v& d$ b7 d6 C
z(1)=y(2);+ J) ?4 W% N4 l
z(2)=sqrt(1+y(2).^2)./(2*(100-x));8 H4 S3 O/ V" W8 \$ q
这是描述一元二阶微分方程组的函数,2 H4 r" i' t U& I
# V2 _+ w+ J/ C( Q8 l" _5 j y0=[0 0];1 Z# y* n+ r0 A0 U! s* F
[x,y]=ode23('rab1',0,99.99,y0);/ l# T) y% o9 T
plot(x,y(:,1),'o') " _3 P- H. y) B2 f8 ^这是写入matlab的代码,为啥运行后显示错误呢,提示说rab函数必须返回一个列向量?- S! D, \1 D4 X- }! t, Z
( B% ^3 _* l ~& p2 A$ Z