zc.m

牛顿迭代法

f1=inline('x^3-x^2-1');
f2=inline('3*x^2-2*x');
[x,k]=niudunfa(f1,f2,2)