function x=matrixLU(A,b,n)
[L,U] = lu(A);
y(1)=b(1);
for i=2:n
    ym=0;
    for j=1:i-1
        ym=ym+y(j).*L(i,j);
    end
    y(i)=b(i)-ym;
end
x(n)=y(n)/U(n,n);
for i=n-1:-1:1
    xm=0;
    for j=i+1:n
        xm=xm+U(i,j).*x(j);
    end
    x(i)=(y(i)-xm)./U(i,i);
end