gexianfa.m

牛顿迭代法

function [x,k]=gexianfa(f,x1,x2,ep)
% 求解非线性方程的割线法，其中
% f(x) --- 表示非线性方程函数
% x1,x2 --- 表示迭代初值
% ep --- 表示精度指标
% k --- 表示循环次数
format long
if nargin <4  ep=1e-5;    end
k=1;
y1=feval(f,x1);
y2=feval(f,x2);
x(k)=x2-(x2-x1)*y2/(y2-y1);
y(k)=feval(f,x(k));
k=k+1;
x(k)=x(k-1)-(x(k-1)-x2)*y(k-1)/(y(k-1)-y2);
while abs(x(k)-x(k-1))>ep
    y(k)=feval(f,x(k));
    x(k+1)=x(k)-(x(k)-x(k-1))*y(k)/(y(k)-y(k-1));
    k=k+1;
end