用matlab写牛顿迭代法

以下是使用matlab实现牛顿迭代法的示例代码：

function [root, fx, n] = newton(func, dfunc, x0, tol, maxiter)

% func是原函数，dfunc是一阶导数函数，x0是初始值，tol是容差，maxiter是最大迭代次数

% 初始化
n = 0;
fx = feval(func, x0);
dfx = feval(dfunc, x0);
x = x0;

% 迭代过程
while abs(fx) > tol && n < maxiter
x = x - fx/dfx;
fx = feval(func, x);
dfx = feval(dfunc, x);
n = n + 1;
end

% 输出结果
if abs(fx) <= tol
root = x;
else
root = NaN;
end

end

在使用该函数时，需要提供以下参数：

- func：原函数的句柄（在matlab中，可以使用@符号创建函数句柄，例如：@sin）
- dfunc：一阶导数函数的句柄
- x0：初始值
- tol：容差
- maxiter：最大迭代次数

例如，如果要求解sin(x) = 0的根，可以这样调用函数：

[root, fx, n] = newton(@sin, @cos, 1.5, 1e-6, 100);

其中，@sin和@cos是sin函数和cos函数的句柄，1.5是初始值，1e-6是容差，100是最大迭代次数。函数会返回求得的根root，最终的函数值fx，以及迭代次数n。