“…Algorithm 1 SQP step0:set β=0.5 ,σ=0.2,ε=1 × 10 −6 ,the initial vector d0 = (1, 1, 1) T ,µ 0 = 0,λ 0 = (0, 0, 0) T ,z 0 = (ε 0 , d 0 , µ 0 , λ 0 ),z 0 = (ε 0 , 0, 0, 0),i = 0 step1:if H(z i ) ≤ 0,stop iteration, else,ψ i = ψ(z i ),ψ i is shown in (34),H(z i ) is shown in (29). step2: Solve the following equations: H(z i ) + H (z i )∆z i = ψz 0 ,then get the solution of the equations:∆z i = (∆ε i , ∆d i , ∆µ i , ∆λ i ) step3:Let m be the smallest non-negative integer m that satisfies the following inequality:…”