please help me to simplify this code bu using matlab. This code is about using Jocobi Method to solbe the 16*16 matrix. And the original question is at third image. thank you so much for helping!!!!!

q function x=Jacobi(A,b,x) n=rank (A): Rmax=Inf: tol=le-5; iter=0; Owhile Rmax>tol x old=x for i=1:n sigma=0; for j=l:n sigma=sigma+A(i, j)*x(j); end x(i)=x_old(i)+(b(i)-sigma)/A(i,i); R(i)=b (i)-sigma; end Rmax=max (abs (R)); iter=iter+l; end c = ['No.of Iteration=', num2 str(iter) ); disp (c); end format long e =(-2 i ooooooooooooo 0; 1 -2 i ooooooooooooo; o1 -2 i oooooooooooo; 1 -2 i ooooooooooo; 1 -2 i oooooooo 1 -2 i ooo 1 -2 i ooooo 1 -2 i oooo 1 -2 i ooo o o o o o o o o 1 -2 i o o o o 0; o o o o o o o o o 1 -2 i o o o 0; oo o o o o o o o o 1 -2 1 1 -2 i o0; o o o o o o o o o o o o 1 -2 1 0; 1 -2 1); b=(-1;-2;-2;-2;-2;-2;-2;-2;-2;-2;-2;-2;-2;-2;-2;-1); x=(0;0;0; 0; 0;0;0;0;0; 0; 0; 0; 0; 0; 0; 0); x=Jacobi(A,b,x); 1. (100 Points] Consider a system of linear equations 0 0 0 0 1 -2 1 0 1 -2 1 0 1 -2 0 0 1 ... ... 0 ... 0 0 1 -2 7 9 a : 9 9 7 0 ... 1 -2 1 0 0 ... 1 -2 1 1 -2 | 0 0 0 ... 0 0 0 1 Solve the system by Gauss elimination Solve the system by TDMA. Solve the system by Jacobi iteration. Solve the system by Gauss-Seidel method. Solve the system by the SOR with w=1.7