Input: m point x 1 , ⋯ , x m and error constant ε > 0
Output: x ( T ) is an 1 − O ( ε ) approximation to model (4)
1: A ← [ ( x 1 ) T , ⋯ , ( x m ) T ] .
2: λ = ε 2 , γ = 1 , and ‖ A x k + b ‖ 2 ‖ b ‖ 2 < 10 − 5 , θ = − 0.25 , s = 90 , t = ( 1 + θ ) 2 s μ λ max ( A T A ) , μ ∈ ( 0 , 4 )
3: x ( T ) , y ( T ) ← A 1 ( A , n , λ , γ , θ , t , s )
4: return x ( T )