作者rooster1125 (rooster)
看板MATLAB
標題[問題]新手勘根問題
時間Tue Apr 24 00:10:25 2012
這學期修了數值分析
第一個作業是二分法 F(x)=(x-1/3)(x-sqrt(3))(x-pi)
教授要我們第一步做roots scanning
以下是我的程式碼]
for a=0:0.1:10
b=a+0.1;
f1=(a-1/3)*(a-sqrt(3))*(a-pi);
f2=(b-1/3)*(b-sqrt(3))*(b-pi);
if f1*f2<0
a
end
end
跑出來的答案是
a =
0.3000
a =
1.7000
a =
3.1000
我的問題是要如何能夠讓答案跑出
c1 =
0.3000
c2 =
1.7000
c3 =
3.1000
這樣的形式
或 c =
0.3000 1.7000 3.1000
向量的形式
以便接下去寫二分法
不好意思 第一次學寫程式
感謝幫忙!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.33.54.146
1F：推 vincent59:for loop外 設1個變數 tt=1; 04/24 01:04
2F：→ vincent59:判斷式裡改為 c(tt)=a; tt=tt+1; 04/24 01:06
3F：→ vincent59:最後 c為1*3的列向量。 04/24 01:07
4F：→ rooster1125:感謝!!! 04/24 01:37
5F：→ rooster1125:但可以請問c(tt)=a; tt=tt+1; 這句話的意思嗎? 04/24 01:45
6F：→ rooster1125:不好意思 我剛翻書了解了! 感謝!! 04/24 01:52