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

henri

原方程为：
& O: @, z( R$ U7 I9 J$ Q
在用ode45解的时候，第一种写法
5 C4 R* [, G, s/ W6 G& Wfunction dy=t1(t,y)
% H3 i# W" g$ \8 `  o- W/ C5 bdy=zeros(2,1);
- _- a0 g2 t- m- Y2 M/ ~1 g# q* my(2)=dy(1);   
dy(2)=1;

, E6 m" m. e. \- r( ?
: |$ q7 S! l0 N9 b. U当我把第三行等式交换位置后，即第三行改为dy(1)=y(2)

最后两种写法得到的结果完全不一样，貌似第一种写法的解是正确的，第二种写法的解是错误的。请教大神，感谢!

Ferrya

- S8 r1 x) _! P7 m
8 y7 s$ b) F2 B! r
先去搞清楚“赋值”的含义，等号在matlab里面是赋值，不是相等，不能左右互换
* _1 Z* i( _+ h( @a = 1; b = 2; c = 3;
a = b;
disp([a b c])
9 {3 ]9 T5 h( f. M/ v! K
4 j, {/ t# C% i0 R7 l4 [% J# Q8 S
结果是     2     2     3

Ferrya

a = 1; b = 2; c = 3;
b = a;
disp([a b c])
0 l9 h7 a% r% X0 T7 j) I
结果是     1     1     3

Ferrya

正确的应当是
function dy=t1(t,y)
dy=zeros(2,1);
# D0 U5 l0 l$ {, ?dy(1)=y(2);
dy(2)=1;