function [x]=LU(A, b)
n=size(A,1);
Q=A
for q=0:n-2
for i=1:n-q
for j=1:n-q
if i==1
if j~=1
Q(i,j)=Q(i,j)/Q(1,1)
end
end
if i~=1
if j~=1
Q(i,j)=Q(i,j)-Q(i,1)*Q(1,j)
end
end
end
end
r=q+1;
for w=1:n-q
t=q+1;
for e=1:n-q
A(r,t)=Q(w,e)
t=t+1;
end
r=r+1;
end
if n-(q+1)>1
Q=ones(n-(q+1),n-(q+1))
r=q+2;
for w=1:n-(1+q)
t=q+2;
for e=1:n-(q+1)
Q(w,e)=A(r,t)
t=t+1;
end
r=r+1;
end
end
end
U=eye(n)
L=eye(n)
for q=1:n
for w=1:n
if q>=w
L(q,w)=A(q,w)
else
U(q,w)=A(q,w)
end
end
end