MATLAB: Is it possible to do this without for cycle

cyclefor

T = zeros(M,K);
for k = 1:K
    T(Y==k,k) = 1;
end

% data
Y = [1 1 3 2 3 1 1 2]  % M = numel(Y)
K = 3
% engine
T = double((Y(:) - (1:K))==0)

ix=([logical(0);(diff(A(1:N))~=0)]);
L(find(id))=A(id);
locais=L(2:end);

I've got a meeting...gotta' run for now, sorry...

Undoubtedly can vectorize remaining fairly much, too.

One thing is it doesn't appear the for k=1:j B(k)=locais(k)'; end B=B(:);

loop does anything that

B=locais;

doesn't other than perhaps the column/vector orientation. I'd say create B to begin with and fix the orientation there as wanted/needed...

MATLAB is 1-based indexing, you light change the for-loop to start at 1

for i= 1:iter_max

or shift the storage index by 1

d(i+1)=d_n;