Initialize U=[uij] matrix, U(0)

At k-step: calculate the centers vectors C(k)=[cj] with U(k)                                 

Update U(k) , U(k+1)                                                     

If || U(k+1) - U(k)||<     then STOP otherwise return to step 2.

-

Initialize U=[uij] matrix, U(0)

At k-step: calculate the centers vectors C(k)=[cj] with U(k)                                 

Update U(k) , U(k+1)                                                     

If || U(k+1)- U(k)||<     then STOP otherwise return to step 2.