> with(LinearAlgebra):
> #Gauss-Elimination mit Ueberschreiben von A, ohne Pivot
> gauss1:=proc(A) 
> local j,k,m,n; 
> # option trace
> n:=RowDimension(A); 
> for m from 1 to n do 
>     for j from m+1 by 1 to n do 
>        A[j,m]:=A[j,m]/A[m,m]; 
>        for k from m+1 by 1 to n do A[j,k]:=A[j,k]- A[j,m]*A[m,k] end do 
>     end do  
> end do; 
> evalm(A) 
> end proc;

> A:=Matrix([[6.,7.,0,1,2,5],[13,4,2,4,5,3],[6,2,3,6,8,9],[1,3,0,0,2,4],[1,1,3,5,7,9],[0,0,4,5,0,2]]);
> LU:=gauss1(A);

>