% compute Ort. projection recursively


function A=ls_rec(A,Syyinv,Syy1inv,y,x,N);


 A = A + x*y'/N*Syy1inv-A*y*y'*Syyinv/(N-1+y'*Syyinv*y);
 