貌似第一种写法的解是正确的，第二种写法的解是错误的。

henri

原方程为：

在用ode45解的时候，第一种写法
function dy=t1(t,y)
dy=zeros(2,1);
% {& w* `0 c2 h0 a+ ]# ly(2)=dy(1);   
: a0 N3 D' T2 S: T  Pdy(2)=1;
2 X/ G1 R. W, J" P
, W0 W* Y* A( N" p! |% s. x
# m$ x! y8 }% f" U  [; _/ i% S; }1 x当我把第三行等式交换位置后，即第三行改为dy(1)=y(2)
  H: k0 S8 N: }4 |
最后两种写法得到的结果完全不一样，貌似第一种写法的解是正确的，第二种写法的解是错误的。请教大神，感谢!

Ferrya

. ?$ a" `. h& e6 K
  j: A' G2 d4 j
; d1 R1 T8 z7 C" h3 |先去搞清楚“赋值”的含义，等号在matlab里面是赋值，不是相等，不能左右互换
1 v" e2 ^  @) k8 S& K5 ta = 1; b = 2; c = 3;
a = b;
& z; h6 m% p* G1 [disp([a b c])
) ~/ `+ d) x1 Z7 C/ V

0 C, c/ h$ ^8 M8 Z' I& p1 q结果是     2     2     3

Ferrya

a = 1; b = 2; c = 3;
# ]9 [1 d' }0 }, M! u7 |' xb = a;
! L+ P3 X+ g- @+ ddisp([a b c])
  @  t$ R9 }5 A! H8 c+ i
结果是     1     1     3

Ferrya

正确的应当是
) Z1 [/ u. f  p6 @* g4 D0 V9 {) nfunction dy=t1(t,y)
2 W; z3 U6 r/ gdy=zeros(2,1);
dy(1)=y(2);
dy(2)=1;