Write some short code implementing the Jacobi method for a tridiagonal matrix

A and a given vector b. Use "for loops", not matrix-vector or matrix-matrix multipli-

cations, and stop after 10, 000 iterations. Try to be as efficient as you can (i.e., your

code should take into account the fact that A is tridiagonal).

What is the flop count of a single iteration? Explain your answer.