Write an R function to solve the following systems of equations using the Jacobi iterative method discussed in class. The code should check to ensure the system will converge and if not, it should report a warning message.

A1 = [3 1 1 0; 1 5 -1 2; 1 0 3 1; 0 1 1 4];

A2 = [2.5 1 1 0; 1 4.1 -1 2; 1 0 2.1 1; 0 1 1 2.1];

A3 = [2 1 1 0; 1 3.5 -1 2; 1 0 2.1 1; 0 1 1 2.1];

b1 = b2 = b3 = [1 4 -2 1]';