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

henri

原方程为：
! R* R% G0 X1 S% D4 k
在用ode45解的时候，第一种写法
function dy=t1(t,y)
dy=zeros(2,1);
y(2)=dy(1);   
2 y% j' J* d5 x! q# {& {- j/ u) n# ^dy(2)=1;
0 N& w7 D" j3 U* ~' N3 h
. F( R1 I: \, ^, r0 e7 i
2 j+ g; j- J+ C. J2 }: m4 |当我把第三行等式交换位置后，即第三行改为dy(1)=y(2)
( a0 n( c( h) G8 W! ^$ B0 Q
最后两种写法得到的结果完全不一样，貌似第一种写法的解是正确的，第二种写法的解是错误的。请教大神，感谢!

Ferrya

2 J$ L5 ]1 W) y: L
, }5 }/ y- X. Y# W+ K先去搞清楚“赋值”的含义，等号在matlab里面是赋值，不是相等，不能左右互换
* w- |- g( l) r; ja = 1; b = 2; c = 3;
7 Z' w" p5 x2 M' M0 }/ a) Q* pa = b;
disp([a b c])
- y# y/ E" i) {0 Y4 i7 Q
) {& K9 B0 Z0 H. [# Q! p! w
结果是     2     2     3

Ferrya

a = 1; b = 2; c = 3;
b = a;
7 N: _" W1 }8 l: ddisp([a b c])
  l/ S* q8 B, X. e/ P
结果是     1     1     3

Ferrya

正确的应当是
6 _1 P( I5 W3 P. {! q+ `# N2 N( H$ ffunction dy=t1(t,y)
; a) O9 }+ W( ]: \" Kdy=zeros(2,1);
dy(1)=y(2);
0 M; a1 e4 A0 t1 tdy(2)=1;