一个方程的求解，数值解 （R语言编写的）

求解方程的数值解，已经有了第一部计算，现在的问题是 迭代计算 ，因为迭代计算我不太懂，所以程序也不会编
f(x)=0
x(k+1)=x(k)-{J(x(k))}^(-1)*f(x(k)), k=0,1,2….
其中J和f已经计算出来了，这个Newtons迭代函数编写如下？
Newtons<-function (funs, y, ep=1e-5, it_max=100){
index<-0; k<-1
while (k<=it_max){
y1 <- y; obj <- funs(y);
y<- y-solve(obj$J, obj$f);

........
其中要迭代的函数：
funs<-function(x){
。。。
f<- c(。。,。。。) #这是f(.)
J<- matrix(,..,….,…,
nrow=2, byrow=T)
# Jacobi矩阵 j(.)
return(
list(f=f, J=J))
}
要迭代部分函数是正确的，运行问题在于：
错误于drop(.Call("La_dgesv", a, as.matrix(b), tol, PACKAGE = "base")) :
  Lapack例行程序dgesv: 系统正好是奇异的
是不是迭代函数存在问题啊，谢谢了

你这个可能是迭代函数的问题，检查下函数，特别是J和F。

看起来有点像是期权定价B-S的东西。不过没看明白你的Rcode。用到numerical analysis的东西，外加Jacobian变换
R有求解的一些package自己试试看！爱莫能助。有完整的code的话或者简单的说明也许我可以帮你想想算法！
建议试一下nls或者uniroot之类的。或在翻翻R的文档看看！

十分感谢，我刚才给您站内短消息了
2# Isscaliu

谢谢分享。