%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%   非线性方程组的数值求解  %%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clc;       % 清除干净屏幕
clear all; % 清除内存

alpha = 1; % 初值1估计值
beta = 1;  % 初值2估计值
eps = 1e-5;   % 计算允许误差

for jj=1:100
    f1=sin(alpha)+3*sin(beta)-2.9873;  % 方程1
    f2=cos(alpha)+2*cos(beta)-1.9142;  % 方程2
    
    if(abs(f1)<eps && abs(f2)<eps)   % 如果解满足，&&表示并且
        alpha   % 在屏幕显示解1
        beta    % 在屏幕显示解2
        break;  % 终止计算
    else   % 否则转入下面的修正计算
        J=[ cos(alpha)  3*cos(beta);
            -sin(alpha) -2*sin(beta)];  %  雅可比矩阵
        dx=-inv(J)*[f1;f2];   % 计算修正值
        alpha=alpha+dx(1);    % 解1进行修正
        beta=beta+dx(2);      % 解2进行修正
    end  % 结束if
end % 结束for