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)
for i=n-1:-1:1
    xm=0;
    for j=i+1:n
        xm=xm+u(i,j).*