Jacobi

Add application "Jacobi Method" for performance evaluation.
b = A[:, -1].reshape((DIM,))