[問題]新手勘根問題

[rooster1125 (rooster)]

這學期修了數值分析

第一個作業是二分法 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
推 vincent59:for loop外 設1個變數 tt=1;
→ vincent59:判斷式裡改為    c(tt)=a; tt=tt+1;
→ vincent59:最後 c為1*3的列向量。
→ rooster1125:感謝!!!
→ rooster1125:但可以請問c(tt)=a; tt=tt+1; 這句話的意思嗎?
→ rooster1125:不好意思 我剛翻書了解了! 感謝!!