Input: m point $x^1\mathrm{,}\cdots \mathrm{,}{x}^m$ and error constant $\epsilon >0$

Output: $x^{\left(\text{T}\right)}$ is an $1-O\left(\epsilon \right)$ approximation to model (4)

1: $A\leftarrow \left[{\left(x^1\right)}^{\text{T}},\cdots ,{\left(x^m\right)}^{\text{T}}\right]$ .

2: $\lambda =\frac{\epsilon }{2},\gamma =1$ , and $\frac{{\Vert A{x}^k+b\Vert }_2}{{\Vert b\Vert }_2}<{10}^{-5},\theta =-0.25,s=90,$ $t=\frac{{\left(1+\theta \right)}^2}{s\mu }{\lambda }_{\max }\left(A^{\text{T}}A\right),\mu \in \left(0,4\right)$

3: $x^{\left(\text{T}\right)},y^{\left(\text{T}\right)}\leftarrow A1\left(A,n,\lambda ,\gamma ,\theta ,t,s\right)$

4: return $x^{\left(\text{T}\right)}$