function [indexes,r] = ownOMP( A, y, K)

r=y;
lambda = zeros(size(A,1),1);
indexes = zeros(size(A,1),1);
Anew = zeros(size(A,1),K);
for k=1:K
    w = (A'*r);
    [~, I] = max(abs(w));
    Anew(:,k) = A(:,I);
    indexes(k) = I;
    A(:,I) = [];
    Anewp = pinv(Anew(:,1:k));
    lambda(1:k) = Anewp * y;
    r = y - Anew(:,1:k)*lambda(1:k);
    display(k);
end

for in = 1:size(A,1)
    indexes(in) = lambda(indexes(in));
end

end