zhidingdiedai.m

来自「牛顿迭代法」· M 代码 · 共 19 行

function [x,k]=zhidingdiedai(f,x0,ep)
% 求解非线性方程的指定迭代法，其中
%   f(x) --- 需要求根的函数
%   x0 --- 初始点
%   ep --- 精度，当|(x(k)-x(k-1)|<ep时，终止计算。省缺为1e-5
%   k --- 迭代次数
format long
if nargin <3  ep=1e-5;    end
k=1;
y0=feval(f,x0);
x(k)=(2-exp(x0))/10;
y(k)=feval(f,x(k));
k=k+1;
x(k)=(2-exp(x(k-1)))/10;
while abs(x(k)-x(k-1))>ep
    y(k+1)=feval(f,x(k));
    x(k+1)=(2-exp(x(k)))/10;
    k=k+1;
end