close all;clear all;clc;
f=@(x)[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))];
x0=[-5;-5];
options=optimset('Display','iter');
[x,fval]=fsolve(f,x0,options);
f=@(x)[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))];
x0=[-5;-5];
options=optimset('Display','iter');
[x,fval]=fsolve(f,x0,options);