请问我下面的程序那里出错了，用的是牛顿环法。
s1=nls(y~b0+x1^b1+x2^b2,data=c1,start=list(b0=1,b1=2,b2=1))
summary(s1)
r自带的结果：
Estimate Std. Error t value Pr(>|t|)
b02.25309 0.30297 7.437 1.78e-09 ***
b10.50375 0.07178 7.018 7.67e-09 ***
b20.10824 0.19833 0.546 0.588
我自己写的程序：
b0=1;b1=2;b2=1;
w=y-b0-x1^b1-x2^b2
exp1=expression(b0+x1^b1+x2^b2)
c2=data.frame(z1=eval(D(exp1,"b0")),z2=eval(D(exp1,"b1")),z3=eval(D(exp1,"b2")))
e0=1;e1=1;e2=1
while(e0^2+e1^2+e2^2>=0.05){
reg2=lm(w~z1+z2+z3-1,data=c2)
e0=coef(reg2)[1];e1=coef(reg2)[2];e2=coef(reg2)[3]
b0=b0+e0;b1=b1+e1;b2=b2+e2;
c2=data.frame(z1=eval(D(exp1,"b0")),z2=eval(D(exp1,"b1")),z3=eval(D(exp1,"b2")))
w=residuals(reg2)
}
b0
b1
b2
运行结果：
b0
0.439291
> b1
1.519478
> b2
0.3713611
>
请高手赐教！
万分感谢！