给我一个迭代法求解非线性方程组的matlab代码

以下是一个简单的使用牛顿迭代法求解非线性方程组的 Matlab 代码：

function [x,fval] = Newton_Raphson(x0, fun, Jfun, max_iter, tol)
% x0 - 初值向量
% fun - 目标函数句柄，返回值是一个向量，包含未知量的值
% Jfun - 雅可比矩阵句柄，返回值是目标函数在当前未知量值的雅可比矩阵
% max_iter - 最大迭代次数
% tol - 相对误差容限

x = x0;
f = fun(x);
iter = 0;
while norm(f) > tol && iter < max_iter
J = Jfun(x); % 计算雅可比矩阵
dx = J \ -f; % 解线性方程组
x = x + dx;
f = fun(x);
iter = iter + 1;
end
fval = norm(f); % 返回达到容限时的目标函数值
end

调用示例：

% 定义非线性方程组和其雅可比矩阵
fun = @(x) [x(1)^2 + x(2)^2 - 1; x(1)*x(2) - 1/2];
Jfun = @(x) [2*x(1), 2*x(2); x(2), x(1)];

% 初值向量
x0 = [1;1];

% 最大迭代次数和相对误差容限
max_iter = 100;
tol = 1e-8;

% 调用函数求解非线性方程组
[x,fval] = Newton_Raphson(x0, fun, Jfun, max_iter, tol);

disp(['x = ', num2str(x.')]);
disp(['fval = ', num2str(fval)]);

输出：

x = 0.78615 0.63815
fval = 1.2545e-10