Mercurial > hg > octave-lyh

diff scripts/linear-algebra/expm.m @ 9819:84398271118c

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
fix typo in expm
author Jaroslav Hajek <highegg@gmail.com>
date Sat, 14 Nov 2009 07:02:30 +0100
parents f0c3d3fc4903
children d1978e7364ad
line wrap: on
line diff
--- a/scripts/linear-algebra/expm.m
+++ b/scripts/linear-algebra/expm.m
@@ -139,7 +139,7 @@
   ## inverse balancing.
   d = diag (d);
   r = d * r / d;
-  r = r(p, p);
+  r(p, p) = r;
   ## Inverse trace reduction.
   if (trshift >0)
     r *= exp (trshift);