On 7/29/07, Daniel Lakeland <dlakelan at street-artists.org> wrote:
> For nonlinear least squares this seems to be a reasonable way to
> go. For linear least squares there are some standard methods using the
> Singular Value Decomposition (SVD) which are numerically stable. I
> don't think we have a numerical SVD in maxima, but perhaps it's time
> to think about getting one, and creating a numerical "generalized
> inverse" function. I know this has been mentioned before.
There is an SVD function (dgesvd) in the lapack package.
See: share/lapack/eigensys.lisp
best
Robert