jacobi iterative method

Solve a system of linear equations, Ax = b, where A is a random, diagonally dominant matrix, and b is A times a vector of ones