Implemented the J_ij matrix.

According to my calculations, the J_ij matrix has been added. It is
O(dimension^2) and thus, there is feedback of printed percentages added
here.