Mercurial > hg > octave-nkf

diff scripts/linear-algebra/expm.m @ 8757:79576d40acb6

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
fix unscaling in expm
author Jaroslav Hajek <highegg@gmail.com>
date Mon, 16 Feb 2009 12:45:15 +0100
parents e07e93c04080
children eb63fbe60fab
line wrap: on
line diff
--- a/scripts/linear-algebra/expm.m
+++ b/scripts/linear-algebra/expm.m
@@ -149,8 +149,8 @@
   endfor
 
   ## inverse balancing.
-  ds = diag (s);
-  r = ds * r / ds;
+  d = diag (d);
+  r = d * r / d;
   r = r(p, p);
   ## Inverse trace reduction.
   if (trshift >0)