用matlab，选取一个初始值x0，用Newton法求解方程x+e^x-2=0，使近似解的误差不超过0.5*10^-8，的完整代码

syms x
f = x*exp(x)-2;
df = diff(f);
x0 = 1; % 初始值
tol = 0.5e-8; % 误差限制
iter = 0;
while abs(subs(f,x,x0)) > tol
x0 = x0 - subs(f,x,x0)/subs(df,x,x0);
iter = iter + 1;
end
disp(['近似解为：',num2str(x0)]);
disp(['迭代次数为：',num2str(iter)]);