同事们，我的表格有两列 - y和x，其中包含此变量的值：
y x
4 3
5 3.5
6.5 4.5
7.6 5.2
9.9 7.4
我需要计算由表格形式确定的此函数的二阶导数。 我知道二阶导数的近似值是：
（y_（i-1）-2y_i + y_（i + 1））/（（x_i-x_i-1）^ 2）
我想为每个x_i得到二阶导数的向量，之后我想打印这个向量。
F.e。在x=3.5点，二阶导数为(4-2*5+6.5)/((3.5-3)^2)
我的输入数据是df=data.frame (y,x)和k - 列数，其中函数（x - 列）的参数是。但我的功能不起作用。
My_fun=function(data,k) {
d1=data[,k] #to get the column with x
d2=data[-k] #to get the column with y
for(i in 2:(nrow(data)-1)) {
c1=(d2[i-1]-d[i]+d[i+1])/((d1[i]-d[i-1])^2)
}
c(NA,c1,NA) #because for first and last element we can not calculate it
}