请教大神，这个是为什么呢？

paprill

% z$ j0 E2 e0 r  ~& t" m$ i( M
请教大神，这个是为什么呢？？ode45 无法执行赋值，因为左侧的大小为 2×4，右侧的大小为 2×5？非常感谢！

function  f=boat1(t,x)
4 h1 c, r6 Z7 a: Y4 ?if(norm(x)>1e-5)   %限制条件，保证在x足够小的情况下可以结束计算
; C, n, Y: J6 B6 U- H2 ?    f=[-2*x(1)/sqrt(x(1).^2+x(2).^2);1-2*x(2)/sqrt(x(1).^2+x(2).^2)]
. j2 e6 B; X7 velse
    f=[0,0]
. c: f# D  x+ u$ k, W
/ V2 l) p9 X& ]& Oend


y0=[1000 0];      %微分方程的初值
[t,x]=ode45(@boat1,[0,1000],y0)      %求解微分方程的数值解
' @* W2 d8 U+ n0 m1 ysubplot(1,2,1),plot(t,x)；         %画出关于x(1)和x(2)的系统时间响应曲线
plot(x(:,1),x(:,2))          %画出关于x(1)和x(2)的系统相平面曲线
axis([-1 1000,0 1000])          %限定区间
! M4 R; d( i3 Y5 |4 g- Utext(10,15,'B')         %标注B点
text(1000,0,'A')         %标注A点
3 k8 J' l7 C' g  K5 hxlabel('X')
7 d4 P( {' b! f0 Sylable(‘Y’)
/ i2 z, L7 y; c( n

无法执行赋值，因为左侧的大小为 2×4，右侧的大小为 2×5。

/ R% i  f' t: X1 \" ^7 l, L/ g" Q! _出错 ode45 (line 488)

& |( R7 c, V# U        yout(:,idx) = yout_new;

Demyar

你的代码有点问题。