for i=1:n
   for j=i:n
       s1=0;
       s2=0;
       for k=1:i-1
           s1=s1+L(i,k)*U(k,j);
           s2=s2+L(j,k)*U(k,i);
       end;
       U(i,j)=A(i,j)-s1;
       L(j,i)=(A(j,i)-s2)/U(i,i);
   end
end
disp('Macierz A:')
A
disp('Macierz L:')
L
disp('Macierz U:')
U

Z = zeros(n);
X = diag(Z);
Z = diag(Z);

for i=1:n
    s = 0;
    for j=1:n
        s = L(i,j)*Z(j) + s;
    end
    Z(i) = Y(i) - s;
end
for i=n:-1:1
    s = 0;
    for j=1:n
        s = U(i,j)*X(j) + s;
    end
    X(i) = (Z(i) - s)/U(i,i);
end
disp('Wyniki:')