k_modes function

K-modes