MATLAB中用solve和fsolve求解非线性方程组

clc;
clear;
syms x2
j=1;
su=0;
x1=2;
jie=solve((x1-x1^(-1)*x2^2)/(1-(x1^2+x2^2+x1^(-2)*x2^(-2)-3)/69)==su,x2);

function F = myfun(x)
F = [2*x(1) - x(2) - exp(-x(1));
-x(1) + 2*x(2) - exp(-x(2))];

x0 = [-5; -5];  % Make a starting guess at the solution
options = optimoptions('fsolve','Display','iter'); % Option to display output
[x,fval] = fsolve(@myfun,x0,options) % Call solver