#include #include # define N 5 jacobi(A,b,x,tol,n,itr) double A[N][N], b[N], x[N], tol[2]; int n, itr[2]; { double LU[N][N], D[N]; double xnew[N], w[N], s1, s2, ep; int i, j, k; for (i=0; i