W_ps=tsalesman(S_in,N_ce,W_ps,Nu_m,Iter,S_wi)

% ---------------------------------------------------------------------
% (see tp3ej1.m)
% ---------------------------------------------------------------------

function W_ps=tsalesman(S_in,N_ce,W_ps,Nu_m,Iter,S_wi)
    Out=zeros([N_ce 1]);
    Datos=2*rand(500,S_in)-1;    
    paux1=[1:N_ce]-fix(N_ce/2); 
    for it=1:Iter
        for i=1:size(Datos,1)
            for j=1:N_ce
                Out(j)=norm(Datos(i,:)-W_ps(j,:));
            end
            [aux1,jm]=min(Out);            
            aux1=abs(circshift(paux1',jm-fix(N_ce/2))');
            aux2=(aux1<=S_wi);
            aux1=exp(-(aux1/(S_wi+1e-4)).^2);
            aux1=(aux1.*aux2)';
            for k=1:S_in
                W_ps(:,k)=W_ps(:,k)+Nu_m.*aux1.*(Datos(i,k)-W_ps(:,k));
            end
            
            
            %--------------------
            statusbar('refresh');
            %--------------------
            
            
        end
    end