Can someone please translate this into C programming?

Function [x] = Jacobi (A, b) %JACOBI the shorthand method that is way easier. n = length (b); U = triu (A, 1) L = tril (A, -1); D = diag(diag (A)); Di = inv (D); eps = 1e-3; err = 1; i = 1; x = zeros (n, 1); while err > eps xn = Di * (b - (L+U)*x); err = max (abs(xn-x)); disp (['Iteration' num|2str (i) 'error' num2 str (err)]); i = i + 1; x = xn; end end