<A NAME=cholesky>
<TITLE>cholesky</TITLE></A>
<b><a href=r37_idx.html>INDEX</a></b><p><p>
<B>CHOLESKY</B> _ _ _ _ _ _ _ _ _ _ _ _ <B>operator</B><P>
<P>
<P>
<P>
<P> <H3>
syntax: </H3>
<em>cholesky</em>(<matrix>)
<P>
<P>
<P>
<matrix> :- a positive definite matrix containing numeric entries.
<P>
<P>
<em>cholesky</em>computes the cholesky decomposition of <matrix>.
<P>
<P>
It returns {L,U} where L is a lower matrix, U is an upper matrix,
A = LU, and U = L^T.
<P>
<P>
<P> <H3>
examples: </H3>
<P><PRE><TT>
F := mat((1,1,0),(1,3,1),(0,1,1));
[1 1 0]
[ ]
f := [1 3 1]
[ ]
[0 1 1]
on rounded;
cholesky(F);
{
[1 0 0 ]
[ ]
[1 1.41421356237 0 ]
[ ]
[0 0.707106781187 0.707106781187]
,
[1 1 0 ]
[ ]
[0 1.41421356237 0.707106781187]
[ ]
[0 0 0.707106781187]
}
</TT></PRE><P>Related functions:
<A HREF=r37_0597.html>lu_decom</A>.
<P>
<P>