non_linear_roots.m

数值分析, 同济大学教材<<现代数值数学和计算>>数值分析课程所有算法的matlab代码,所有例程均经过测试,私家珍藏

function non_linear_roots(fun, x0)
tol = 1e-6;
syms x;
it = 0;
x = x0 + tol + tol;
while abs(x - x0) > tol
  it = it + 1;
  if (it > 10000) || (x == Inf) || (x == -Inf) || (x == NaN)
    fun
    disp('Iteration Failed.');
    return;
  end
  x0 = x;
  x = subs(fun, x);
end
fun
x
it, disp('Iteration Success.');