linalg[leastsqrs] in Maple V R4

I know that I am using an admittedly ancient version of the software,
but it suits my modest personal purposes for the most part.

I am wondering if anyone has has any experience with the leastsqrs
procedure, in either older or newer incarnations, only to discover that
it seems painfully inaccurate on repeat application.

I have tried to use linalg[leastsqrs] to find the least squares
solution of an fairly hefty sized overdetermined system, about 80
equatons in 16 unknowns. Repeat application on the exact same
matrix/vector input combination yields different results with the loss
of many if not most of the significant digits. I am not talking just a
few digits lost to accumulated rounding error, I am talking a couple of
dozen digits when Digits is set pretty high. The problem gets worse the
larger the linear system is, and using the 'optimize' argument seems to
slow things down and doesn' add any value.

This is driving me wonky. Is this a bug or a feature? Is it present in
newer versions of Maple? Would I be better off using Svd() and
computing the least squares solution myself from that output? I really
don't mind "black box" routines if they give reproducible results, and
linalg[leastsqrs] in my hands does not.

Many thanks in advance,

Les

.

Relevant Pages

  • Re: Is the question answered?
    ... the reason it has not is that Mr. bimms is unwilling to ... All he can do is repeat his ... To claim that the digits of pi are random is at best completely misleading. ... If suddenly the pool balls started jumping randomly around the table, ...
    (talk.origins)
  • Re: Assembly Language - Mathematics WITHOUT maths coprocessor
    ... My math teacher made the mistake of saying that you could always tell these numbers apart from numbers that never repeat, because the ones that would repeat would repeat after just a few digits. ... I went home and wrote a little program on my TRS-80 to try to find a closer fraction, and I eventually came up with 355/113, which is accurate to six decimal places. ...
    (alt.lang.asm)
  • Re: Divisibility by 7 & 13
    ... This is but one method to check divisibility by 7. ... Again, in this method, we are to repeat the ... process till all the digits are exhusted. ... mathematicians to make sure it is correct; if it is either unknown or ...
    (sci.math)
  • Re: Number Theory
    ... What are the possible remainders when dividing a whole number by 7? ... decimal answer doesn't repeat after 8 digits, ... to be, why isn't the period longer than 6; you've taken it to be, why ...
    (sci.math)