Try this and let me know how poorly / well it works:
(%i50) load(linearalgebra)$
(%i53) moore_penrose_pseudoinverse(matrix([a,b],[2*a,2*b]));
(%o53)
matrix([a/(5*b^2+5*a^2),(2*a)/(5*b^2+5*a^2)],[b/(5*b^2+5*a^2),(2*b)/(5*b^2+5*a^2)])
(1) Unlike most functions in linearalgebra, moore_penrose_pseudoinverse
doesn't autoload;
assuming that moore_penrose_pseudoinverse is at least useable, that's a
bug.
(2) There is no user documentation for moore_penrose_pseudoinverse; and
that's a bug too.
Barton
maxima-bounces at math.utexas.edu wrote on 08/24/2009 11:51:17 AM:
>
> Does anyone know if there is some implementation for calculating a
> pseudo inverse in Maxima? I've found this link - "http://
> www.math.utexas.edu/pipermail/maxima/2007/008247.html" - but it
> seems that the implementation is not exactly for Maxima and it only
> handle numerical matrices...