function [L,U] = lubasic(A) 

% LU razcep matrike A brez pivotiranja,
% v tej osnovni varianti imamo tri vgnezdene zanke,
% delamo s posameznimi elementi

% V starih verzijah Matlaba je bila ta koda pri izklopljeni optimizaciji z
% feature accel off pocasnejsa od blocnih variant lubasic2 in lubasic3
% ki namesto BLAS1 uporabljata BLAS2 oz. BLAS3

% Bor Plestenjak

n = size(A,2);
for j = 1 : n-1 
   for i = j+1 : n
      A(i,j) = A(i,j)/A(j,j);
      for k = j+1 : n
         A(i,k)=A(i,k)-A(i,j)*A(j,k);
      end
   end
end   
U = triu(A);
L = eye(n)+tril(A,-1);

