写一段MATLAB代码解不动点迭代法解非线性方程组function [ f,j ] = funs( x ) f=[x(1)^2+x(2)^2-5;(x(1)+1)*x(2)-(3*x(1)+1)]; j=[2*x(1) 2*x(2);x(2)-3 x(1)+1];

首先，我们需要定义一个不动点迭代法的函数：

function [x,iter] = fixed_point_iteration(f,x0,tol,max_iter)
% f: 函数句柄，表示需要求解的非线性方程组
% x0: 初始解向量
% tol: 迭代收敛的容忍度
% max_iter: 最大迭代次数

iter = 0;
x = x0;

while iter < max_iter
    x_new = f(x); % 计算新的解向量
    if norm(x_new - x) < tol % 判断是否满足收敛条件
        break;
    end
    x = x_new;
    iter = iter + 1;
end

end

然后，我们可以调用该函数来解决非线性方程组：

[x,iter] = fixed_point_iteration(@(x)funs(x),[0;0],1e-6,100);
disp(['解向量为：[',num2str(x(1)),',',num2str(x(2)),']']);
disp(['迭代次数为：',num2str(iter)]);

在这里，我们将 `funs` 函数作为参数传递给了不动点迭代法的函数。最终，我们可以得到方程组的解向量和迭代次数。