function z=rab(x,y)2 ^- Y( n* W' Z, N. ]
z(1)=y(2);& t* d: }6 f' i, F
z(2)=sqrt(1+y(2).^2)./(2*(100-x)); y( A& t! ^3 i这是描述一元二阶微分方程组的函数,1 G8 ?' d) M) ?8 o' Z9 n
F! L8 _6 j" y. g5 l: o7 c1 Y/ Q y0=[0 0]; 5 y& M4 M8 z# S. s* c[x,y]=ode23('rab1',0,99.99,y0);4 c s6 u3 F0 \ |2 ]7 b
plot(x,y(:,1),'o') # H$ Y( A0 c# j$ h这是写入matlab的代码,为啥运行后显示错误呢,提示说rab函数必须返回一个列向量?5 D. ?" H9 W! s; x9 p0 N