Input: A , B , w ( 0 ) = [ z 1 ( 0 ) , z 2 ( 0 ) , α ( 0 ) , β ( 0 ) ] T , k m a x .

for k = 0 , 1 , 2 , , until convergence

Choose τ ( k ) .

Solve M ( k ) u ( k ) = B 1 w ( k ) inexactly,

P 1 M ( k ) u ( k ) P 1 B 1 w ( k ) τ ( k ) .

end for

Output: u ( k m a x ) .