
function[w,y,e,J]=aaselfcorrectinglms(x,dn,mu,M,I)
%function[w,y,e,J]aaselfcorrectinglms(x,dn,mu,M,I);
[w(1,:),y(1,:),e(1,:)]=aalms(x,dn,mu,M);
for i=2:I%I=number of iterations, I<8-10 is sufficient;
    [w(i,:),y(i,:),e(i,:)]=aalms(y(i-1,:),dn,mu,M);
end;
J=e.^2;%J is an Ixlenght(x) matrix;
