请问大家，e初始值迭代范围为［0.5，2］，迭代步长为0.01，这里怎么改stata代码啊，有大佬做过吗，或者有大佬用matlab估计过CES供给面标准方程吗，有偿求代码

nlsur (lnY_norm = log({xi}) + 1/(({e}-1)/{e})*log(((1-{a})*(nL_norm*exp(t_bar*{gammaL}/{lammaL}*(t_norm^{lammaL}-1)))^(({e}-1)/{e})) + ({a}*(nK_norm*exp(t_bar*{gammaK}/{lammaK}*(t_norm^{lammaK}-1)))^(({e}-1)/{e})))    ) ///  (lLs = log(1-{a})  +(({e}-1)/{e})*log({xi}) - (({e}-1)/{e})*log(nY_norm/nL_norm) + (({e}-1)/{e})*(t_bar*{gammaL}/{lammaL}*(t_norm^{lammaL}-1))) ///
(lKs = log({a})  +(({e}-1)/{e})*log({xi}) - (({e}-1)/{e})*log(nY_norm/nK_norm) + (({e}-1)/{e})*(t_bar*{gammaK}/{lammaK}*(t_norm^{lammaK}-1))) ///
, initial(e 0.7 xi 1 lammaL 1 lammaK 1  a 0.5)