作者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