MATLAB学习笔记之牛顿迭代法

牛顿迭代法的公式如下

给定函数f(x)和近似根值x0，由上式逐步逼近根值；直到误差小于给定值；
代码如下：

function x1=newton(fun,x0,delta)
syms x;%定义符号变量x，用于求导和代入
df=diff(fun(x));%对f(x)进行求导
x1=x0-fun(x0)/subs(df,x0);
while abs(x1-x0)>delta
x0=x1;
x1=x0-fun(x0)/subs(df,x0);
end

这里提一下subs函数，subs函数用于替换表达式中的符号变量，在第四行代码中，

x1=x0-fun(x0)/subs(df,x0);

即用x0的值来替换df表达式中的变量，否则编译器会报错，提示“ 未定义的变量值 ”。