#=================================================#
# Dies ist der Maple source zum                   #
# Cholesky-Verfahren zur L-D-L^T-Zerlegung        #
#                                                 #
# Erwartet wird, dass die Dimmension n, die       #
# Matrix a, und der Diagonalvektor d bereits      #
# sind. Der unter untere Teil von a wird mit L    #
# ueberschrieben.                                 #
#=================================================#

# cholesky-Zerlegung
for j from 1 to n do 
  d[j]:=a[j,j]-sum(a[j,jjs]^2*d[jjs], jjs=1..j-1);
  for i from j+1 to n do
    a[i,j]:=(a[i,j]-sum(a[i,jjs]*d[jjs]*a[j,jjs], jjs=1..j-1))/d[j];
  od:
od:  
a:=evalm(a);   # zeige das Zwischenergebnis
d:=evalm(d);
