clear all;
n=5;
A=rand(n,n) % Zufallsmatrix 
[L, R]=LRnaiv(A); 
L
R
rest=A-L*R