Re: Implementing a Jacobi iterative method for Ax=b

[Juan Pablo Carbajal]

On Mon, Oct 29, 2012 at 2:16 PM, Joza <address@hidden> wrote:
> Here are A and b if you want to test it out!
>
> A = [ 3 12 0 -1 0 0; 4 0 31 1 0 0; 2 1 0 0 17 -3; 27 2 0 0 0 1; 0 0 0 -1 1
> 11; 0 0 0 24 -1 0 ]
>
> (6x6)
>
> b = [ 39 117 12 98 14 55]^T
>
> (6x1)
>
>
>
> --
> View this message in context: 
> http://octave.1599824.n4.nabble.com/Implementing-a-Jacobi-iterative-method-for-Ax-b-tp4645833p4645868.html
> Sent from the Octave - General mailing list archive at Nabble.com.
> _______________________________________________
> Help-octave mailing list
> address@hidden
> https://mailman.cae.wisc.edu/listinfo/help-octave

Hi,

Use Agora to post your code, it is easier to track the changes
http://agora.octave.org/

I posted an improved version of your file.
http://agora.octave.org/snippet/JZaP/

You must always check if the algorithm will converge, in order to know
if you should expect an answer or not. Try
demo jacobi_iter

Cheers