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.