From help-request at octave dot org Fri Oct  7 06:37:12 2005
Subject: Re: Problem of polyfit
From: Tetsuro KURITA <tkurita at mac dot com>
To: Fredrik Lingvall <fl at ifi dot uio dot no>
Cc: Octave_post Help <help at octave dot org>
Date: Fri, 7 Oct 2005 20:36:11 +0900


On 2005/10/07, at 17:32, Fredrik Lingvall wrote:

>
> Can you post the output of  norm( inv(X'*X)*X'*y - (X'*X)\X'*y )?

The result is as follows.

octave> norm(inv(X'*X)*(X'*y) - (X'*X)\(X'*y))
warning: matrix singular to machine precision, rcond = 3.02246e-112
warning: attempting to find minimum norm solution
ans = 50.504

When calculating (X'*X)\(X'*y), warning message raise.


=======================================================
  Tetsuro KURITA
   E-mail: tkurita at mac dot com
=======================================================



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------