看板Statistics
标 题用 R 写二分法的问题
发信站无名小站 (Tue Mar 13 23:39:31 2007)
转信站ptt!Group.NCTU!grouppost!Group.NCTU!wretch
delta0=0
delta1=1
while(abs(
L(delta0,B)$Ld-delta0)>0.00001) 黄色是我写好的fun.
{
phi1=delta0
phi2=delta1
phim=(phi1+phi2)/2
if((L(phi1,B)$Ld-phi1)*(L(phim,B)$Ld-phim)<0)
{
delta0=phi1
delta1=phim
}
else
{
delta0=phim
delta1=phi2
}
if (abs( L(phim,B)$Ld-phim )<0.0001)
break
}
delta0
delta1
---------------- 以上是我的code --------------------------------
我的问题是...
> delta0
[1] 0.9972603
> delta1
[1] 0.9972603
(delta0,delta1是我按stop才停止程式得到的值)
怎麽会 delta0 = delta1 呢??
都是一直 相加除以2 下去~ 怎麽会出现相等呢??
是不是 R 只取到小数点以下第七位~ 要怎麽修改呢?
谢谢
--
夫兵者不祥之器物或恶之故有道者不处君子居则贵左用兵则贵右兵者不祥之器非君子
之器不得已而用之恬淡为上胜而不美而美之者是乐杀人夫乐杀人者则不可得志於天下
矣吉事尚左凶事尚右偏将军居左上将军居右言以丧礼处之杀人之众以哀悲泣之战胜以
丧礼处之道常无名朴虽小天下莫能臣侯王若能守之万物将自宾天地相合以降甘露民莫
之令而自均始制有名名亦既有夫亦将知止知止59-115-229-233.dynamic.hinet.net海