Introduction - If you have any usage issues, please Google them yourself
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.