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

henri

原方程为：

在用ode45解的时候，第一种写法
function dy=t1(t,y)
" N; f7 N6 a. [+ R+ S  y" Fdy=zeros(2,1);
y(2)=dy(1);   
dy(2)=1;


当我把第三行等式交换位置后，即第三行改为dy(1)=y(2)
0 L6 n0 w. v% o6 {# U& z( z
最后两种写法得到的结果完全不一样，貌似第一种写法的解是正确的，第二种写法的解是错误的。请教大神，感谢!

Ferrya

# X1 G8 s+ c' R

1 s' @) |/ R- B) t8 I先去搞清楚“赋值”的含义，等号在matlab里面是赋值，不是相等，不能左右互换
  R/ {+ t( M" B) Q2 ba = 1; b = 2; c = 3;
/ l4 `) B; {# P" h  K5 q! Ua = b;
disp([a b c])
7 c! v9 b! L+ Z$ Y& z
" g+ D8 ?, x+ D4 E. o' C* A
结果是     2     2     3

Ferrya

a = 1; b = 2; c = 3;
b = a;
disp([a b c])
* _3 _% @+ j0 i% Z( P, ^! K
结果是     1     1     3

Ferrya

正确的应当是
function dy=t1(t,y)
dy=zeros(2,1);
dy(1)=y(2);
$ ~% U$ I8 v  Z; m# Sdy(2)=1;