c <- function(i){
(5-i)/4
}
p <-function(x){
for(i in 6:m){
k <-1:4
1-sum(c((c(1)-x)^4/(c(1)*prod(c(1)-c(k[-1]))),
(c(2)-x)^4/(c(2)*prod(c(2)-c(k[-2]))),
(c(3)-x)^4/(c(3)*prod(c(3)-c(k[-3]))),
(c(4)-x)^4/(c(4)*prod(c(4)-c(k[-4]))),
sum((c(i)-x)^4/(c(i)*prod(c(i)-c(k))))))
}
}
i是不能等于5的 不然 c(5)=0, 1/c(5)= inf
程序有问题 改改ing!!!!!!!!!!!!!!!!!!