Re: Implementing a Jacobi iterative method for Ax=b

[c.]

On 28 Oct 2012, at 17:59, Juan Pablo Carbajal wrote:

> do not build the matrix, just keep the vector
> x_new = (b - N*x_old)./M;

Octave is clever enough to do that itself, diagonal matrices are a special 
class:

>> typeinfo (diag (diag (M)))
ans = diagonal matrix


and have their own methods.
c.