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

henri

原方程为：
2 q' K0 q% {  h+ f( ~
在用ode45解的时候，第一种写法
function dy=t1(t,y)
dy=zeros(2,1);
1 \- b- p& p2 @y(2)=dy(1);   
( d- M/ }  T2 p% Y; A4 l7 wdy(2)=1;


7 y; s+ h! o$ n6 A) d: ?1 D3 k当我把第三行等式交换位置后，即第三行改为dy(1)=y(2)

( |. M; \9 f3 W9 X( E, `( ^$ k9 M最后两种写法得到的结果完全不一样，貌似第一种写法的解是正确的，第二种写法的解是错误的。请教大神，感谢!

Ferrya

+ V- v3 \0 B" ^( |$ h/ U

先去搞清楚“赋值”的含义，等号在matlab里面是赋值，不是相等，不能左右互换
* W  a6 L% B  B2 n1 Ra = 1; b = 2; c = 3;
0 t3 G$ g( p/ `6 k4 o8 J6 {a = b;
disp([a b c])
$ @8 z/ F' x: ^, b$ [
: _9 I; r; k/ C3 ?9 D7 H* B
结果是     2     2     3

Ferrya

a = 1; b = 2; c = 3;
b = a;
disp([a b c])
* a" Q  c+ G# v  B* c0 e# G& M; z6 E9 g
结果是     1     1     3
' r2 V  q' Z1 |6 H  T  L

Ferrya

正确的应当是
function dy=t1(t,y)
dy=zeros(2,1);
dy(1)=y(2);
4 g/ P4 f: D( T- K& X' ]+ r/ g5 Z. ldy(2)=1;