function x=eqznormali(A,b)

X=A'*A;
[L,U] = lu(X);
b1=A'*b;
D=diag(diag(U));
sqrtD=sqrt(D);
L=L*sqrtD;
y=L\b1;
x=L'\y;