代码
root<-function(x){
a=0
b=1
while( a != b ){
a = b
b=0.5*(b + x/b)
}
b
}
问题
上面这段代码用于求x的root值。在这里有一个十分不专业的问题, 那就是a != b的写法。
对于numeric的变量来说,使用 != 或者是 = = 是非常不专业的。针对整型数,这会返回正确的结果,但是由于机器精度,对于浮点数来说,他就“可能”会返回一个错误结果,于是在上述代码中,该问题则会导致程序陷入While的无限循环。
解决方式
abs(a-b)<exp(-30)
对于向量形式的两个变量,则可使用 max(abs(a-b))<exp(-30) 来返回两个向量中的所有元素是否完全一致
PS:对于exp(n)中的n,可以按需选用。例如说,如果两个浮点数的差距只有10的负十几次方的话,n可以改为-40,-50,-60…,省事一点的话,可以直接选到-100之类的,省去麻烦。