for k = 1 : n-1
		for i = k+1 : k+l-k%l
			Iik = A[i,k] / A[k,k]

			A[i,k] = 0.0 #Iik # mnozniki - rozklad LU zapisany w jednej macierzy
			for j = k+1 : min(k+l, n)
				A[i,j] = A[i,j] - Iik * A[k,j]
			end
			b[i] = b[i] - Iik * b[k]
		end
	end
	return (A, b)