From help-octave-request at bevo dot che dot wisc dot edu Sat Jan  4 03:01:26 2003
Subject: solving Ax=b
From: Sven Khatri <skhatri at cds dot caltech dot edu>
To: help-octave at bevo dot che dot wisc dot edu
Date: Sat, 4 Jan 2003 01:00:26 -0800

Hi All,

I'm looking to solve a problem that is equivalent to solving for x in
the system of linear equations:
  Ax = b
where A \in R^{n \times n}, x \in R^n, b \in R^n. But the catch is
that n is too large to explicitly construct A (even in Matlab's sparse
implementation) but I can construct a function f so that f(x) = Ax. Is
there an easy way (within octave) to solve for x? I can solve the problem
by introducing \tilde{A} = I - A, so that the solution to 
  x = \tilde{A}x + b
is a solution to the above problem and \tilde{A} is a contraction and
solving the problem iteratively but this computation is VERY slow.

thanx...
Sven

PS I hope you all are comfortable with the above latex notation



-------------------------------------------------------------
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
-------------------------------------------------------------