function [centros,lrate]=clust(InpPat,Tarpat,n_centros)

[l,c]=size(InpPat);

% Valores iniciais para o algoritmo
[centros,vk,lr]=initoakmr(InpPat,n_centros,eps);
  
for x=1:l,
   [centros,vk,lr]=optkmeansr(InpPat(x,:),centros,vk,lr);
   lrate(x)=lr;
end