[問題] 迴圈取值

[easy1612 (up)]

請問各位大大
小弟得程式碼如下

1.clc;
2.n=2.168;
3.X=[];
4.b=[];
5.for L=0:11
6.M21=[1 L ; 0 1];
7.M22=[1 0; 0 1];
8.M23=[1 11-L;  0 1];
9.M24=[1 0 ; 0 1/n];
10.M25=[1 8/n ;0  1 ];
11.M26=[1 0 ;0 n];
12.M27=[1 5; 0 1];
13.M28=[1 0; -2/150 1];
14.M2t=M21*M22*M21*M23*M24*M25*M26*M27*M28*M27*M24*M25*M26*M23;
15.A=M2t(1,1);
16.B=M2t(1,2);
17.C=M2t(2,1);
18.D=M2t(2,2);
19.det(M2t);
20.q1=(A-D-((D-A)^2+4*B*C)^0.5)/(2*C);
21.q11=1/q1;
22.R=1/real(q11);
23.W=(-(1.064*(10^-18))/(pi*(imag(q11))))^0.5;

24.b=[L,W]

25.X=[X,b'];
26.end
27.X

理想為迴圈內取出L和W

問題出在24行
出來的b值只有L，W都是0
確定W是有值

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.42.207
※ 編輯: easy1612        來自: 140.115.42.207       (01/01 22:13)
推 sunev:把23行的分號去掉，看看情況，再把21的也去掉。
→ easy1612:修正過後結果與修正前相同
→ Gwaewluin:結果當然是相同的，去掉分號只是要你檢查途中的計算結果