用MATLAB二分法求解非线性方程组，不用fzero函数

二分法是一种简单的数值求解方法，可以用于求解非线性方程组。其基本思想是：将待求解区间逐步缩小，直到满足一定精度要求为止。具体步骤如下：

1. 选定一个初始区间[a, b]，使得f(a)和f(b)异号（即方程组在a和b处的函数值符号不同）。

2. 将区间[a, b]平均分成两部分，得到中点c=(a+b)/2。

3. 计算f(c)的符号，如果f(c)的符号与f(a)相同，则将区间[a, c]作为新的待求解区间；否则将区间[c, b]作为新的待求解区间。

4. 重复步骤2和3，直到得到满足一定精度要求的解。

以下是使用MATLAB实现二分法求解非线性方程组的示例代码：

%% 定义方程组的函数
function y = fun(x)
y(1) = x(1)^2 + x(2)^2 - 1;
y(2) = x(1) - x(2)^2;

%% 初始化参数
a = [0, 0]; % 初始区间左端点
b = [1, 1]; % 初始区间右端点
tol = 1e-6; % 精度要求

%% 进行二分法迭代
while norm(b-a) > tol
c = (a+b)/2; % 中点
if sign(fun(c)) == sign(fun(a))
a = c;
else
b = c;
end
end

%% 输出结果
disp(['x1 = ', num2str(c(1)), ', x2 = ', num2str(c(2))]);