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

paprill

. O- ~2 C. _5 S% ^; z( }
; }) r7 N5 f5 }' e  Z9 K! x请教大神，这个是为什么呢？？ode45 无法执行赋值，因为左侧的大小为 2×4，右侧的大小为 2×5？非常感谢！
$ _8 p! @: f( {. u6 g! @1 F1 H  G. d
% u3 }+ B1 i. z. \; I4 Lfunction  f=boat1(t,x)
if(norm(x)>1e-5)   %限制条件，保证在x足够小的情况下可以结束计算
    f=[-2*x(1)/sqrt(x(1).^2+x(2).^2);1-2*x(2)/sqrt(x(1).^2+x(2).^2)]
else
* s7 h" a! L/ U( \+ B: \    f=[0,0]
0 r" r5 n) |9 V2 H1 I
end


& l3 L. B9 }" {% ]- r2 ?, }y0=[1000 0];      %微分方程的初值
[t,x]=ode45(@boat1,[0,1000],y0)      %求解微分方程的数值解
subplot(1,2,1),plot(t,x)；         %画出关于x(1)和x(2)的系统时间响应曲线
0 P9 k0 S  u$ t$ X; Iplot(x(:,1),x(:,2))          %画出关于x(1)和x(2)的系统相平面曲线
axis([-1 1000,0 1000])          %限定区间
text(10,15,'B')         %标注B点
text(1000,0,'A')         %标注A点
! B0 Z& E* u, G" \3 o5 gxlabel('X')
+ r' C" ~5 S5 i( F; o0 Tylable(‘Y’)

4 {" F3 D" e6 ~/ m! h, P) l- n
无法执行赋值，因为左侧的大小为 2×4，右侧的大小为 2×5。
2 G, F* t$ i! H0 Z0 i* D2 g
出错 ode45 (line 488)
! n& F: h, f7 D4 F
        yout(:,idx) = yout_new;

Demyar

你的代码有点问题。