Jacobi Method; Simple system of equations

Click For Summary
SUMMARY

The discussion focuses on the Jacobi Method for solving a simple system of equations, specifically targeting the solution of x=2 and y=1. The coefficient matrix T is defined as T = [[0, -1/2], [1, 0]] and the constant vector c as c = [2.5, -1]. The iterative formula x_{k+1} = T x_k + c is utilized, starting from the initial guess of (0,0). It is noted that achieving accuracy to within 2 decimal places requires approximately 19 iterations, indicating a slow convergence rate.

PREREQUISITES
  • Understanding of the Jacobi Method for iterative solutions
  • Familiarity with matrix representation and operations
  • Basic knowledge of MATLAB for programming iterations
  • Concept of convergence in numerical methods
NEXT STEPS
  • Implement the Jacobi Method in Python using NumPy
  • Explore convergence criteria for iterative methods
  • Learn about alternative methods such as Gauss-Seidel for comparison
  • Investigate the impact of different initial guesses on convergence speed
USEFUL FOR

Students and professionals in mathematics, engineering, and computer science who are working with numerical methods for solving systems of equations, particularly those interested in iterative techniques like the Jacobi Method.

RogerDodgr
Messages
20
Reaction score
0
I know this is simple, and I am missing something obvious. I'm suposed to use the "jacobi method"; and with each iteration it should be getting closer and closer to the solution (x=2 and y=1, which it is not). Could someone explain what I'm doing wrong, or how to start?
http://www.sudokupuzzles.net/IMG_0031.jpg
 
Last edited by a moderator:
Physics news on Phys.org
If you do it using matrices, you'll get the coefficient matrix

[tex]\displaystyle T=\begin{pmatrix} 0 & -\frac{1}{2}\\ 1 & 0 \end{pmatrix}[/tex]
and
[tex]c = \begin{pmatrix} 2.5\\-1\end{pmatrix}[/tex]

Then you can evaluate each iteration via

[tex]\displaystyle x_{k+1} = T x_k + c[/tex]

where k+1 is your iteration number, and by looking at your work, you've chosen (0,0) to be your initial guess. I'm not too sure the numbers match up with what you've shown, but give it a try.
 
Using a small program in MatLab, I found that in order to be accurate to within 2 decimal places, it requires about 19 iterations. 3 decimal places took 26 iterations, and 4 decimal places takes about 46. Needless to say, it doesn't converge very quickly.
 

Similar threads

Replies
8
Views
2K
  • · Replies 2 ·
Replies
2
Views
4K
Replies
10
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 5 ·
Replies
5
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
2
Views
1K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K