%Header Record
Format:TXT
Communication SW:0
Data Type:PG
Capacity:807
File Name:LR ZERLG
Group Name:
Password:
Option1:NL
Option2:
Option3:
Option4:
%Data Record
' QUADR. MATRIX D
' D...A
' LOESUNG ALS
' E...L
' F...R
' MIT L*R=A
\ClrText
" "
" "
" "
\Locate 1,1,"1\->LR NACH GAUSS"
\Locate 1,2,"2\->LR NACH CROUT"
\Locate 1,3,"3\->LR NACH CHOLESKY"
?\->A
A=2\=>\Trn \Mat D\->\Mat D
\Dim \Mat D\->\List 1
\List 1[1]\->M
\List 1[2]\->N
\Identity M\->\Mat F
\Identity M\->\Mat E
\If A\<>3
\Then \For 1\->K \To N
\For K\->J \To N
0\->S
\For 1\->L \To K-1
S+\Mat E[K,L]*\Mat F[L,J]\->S
\Next 
\Mat D[K,J]-S\->\Mat F[K,J]
\Next 
\If K<M
\Then \For K+1\->I \To M
\If \Mat F[K,K]=0
\Then 0\->\Mat E[I,K]
\Else 0\->S
\For 1\->L \To K-1
S+\Mat E[I,L]*\Mat F[L,K]\->S
\Next 
(\Mat D[I,K]-S)/\Mat F[K,K]\->\Mat E[I,K]
\IfEnd
\Next 
\IfEnd
\Next 
\If A=2
\Then \Trn \Mat E\->\Mat Z
\Trn \Mat F\->\Mat E
\Mat Z\->\Mat F
\Trn \Mat D\->\Mat D
\IfEnd
\Else 'A=3
\For 1\->J \To N
\For J\->I \To M
0\->S
\For 1\->L \To J-1
S+\Mat E[I,L]*\Mat E[J,L]\->S
\Next 
\If I=J
\Then \If \Mat D[I,J]<S
\Then "NICHT MOEGLICH"
\Stop
\Else \sqrt(\Mat D[I,J]-S)\->\Mat E[I,J]
\IfEnd
\Else \If \Mat E[J,J]=0
\Then 0\->\Mat E[I,J]
\Else (\Mat D[I,J]-S)/\Mat E[J,J]\->\Mat E[I,J]
\IfEnd
\IfEnd
\Next 
\Next 
\Trn \Mat E\->\Mat F
\IfEnd
\Mat E*\Mat F\->\Mat Z
\For 1\->I \To M
\For 1\->J \To N
\If \Mat Z[I,J]\<>\Mat D[I,J]
\Then \Identity 1\->\Mat Z
"FEHLER"
\Stop
\IfEnd
\Next 
\Next 
\Identity 1\->\Mat Z
"OK"

%End