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