我想用R语言对两个向量的所有元素的乘积求和,但出了点问题。这是我的数据定义:
> alpha <- 1/24
> a <- c(-5, -2, 1)
> b <- c(alpha*3, alpha*2, 1-5*alpha)
然后我在尝试:
> result <- sum(a*b)
> result
[1] 8.326673e-17
我希望结果为零,但事实并非如此 为了更好地解释,我想计算这个总和:
(a[1]*b[1]) + (a[2]*b[2]) + (a[3]*b[3])
应该是:
(-5*alpha*3)+(-2*alpha*2)+(1*(1-alpha*5)) = (-5*1/24*3)+(-2*1/24*2)+(1*(1-5*1/24)) = 0