Input:, , , for, and 0 otherwise.

While:

Compute by Equation(15);

Update as the eigenvectors corresponding to the largest eigenvalues of the matrix;.

Update;.

;

Output resultant and.