Mercurial > octave
diff scripts/linear-algebra/logm.m @ 4334:764229f9a5c8
[project @ 2003-02-19 06:24:02 by jwe]
author | jwe |
---|---|
date | Wed, 19 Feb 2003 06:25:09 +0000 |
parents | |
children | 4c8a2e4e0717 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/scripts/linear-algebra/logm.m Wed Feb 19 06:25:09 2003 +0000 @@ -0,0 +1,36 @@ +## Copyright (C) 2003 John W. Eaton +## +## This file is part of Octave. +## +## Octave is free software; you can redistribute it and/or modify +## it under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## Octave is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. +## +## You should have received a copy of the GNU General Public License +## along with this program; see the file COPYING. If not, write to the +## Free Software Foundation, 59 Temple Place - Suite 330, Boston, MA +## 02111-1307, USA. + +## -*- texinfo -*- +## @deftypefn {Function File} {} logm (@var{a}) +## Compute the matrix logarithm of the square matrix @var{a}. Note that +## this is currently implemented in terms of an eigenvalue expansion and +## needs to be improved to be more robust. +## @end deftypefn + +function B = logm (A) + + if (nargin != 1) + usage ("B = logm (A)"); + endif + + [V, D] = eig (A); + B = V * diag (log (diag (D))) * inv (V); + +endfunction