Answer the following question using MATLAB or octave and provide your code:

Modify the function file gauss seidel.m to produce a new function file that implements the Jacobi method and use it to solve:

12x₁ 3x₂ +4x₃ 2x₄ = 12

2x₁ +10x₂ x₃ - 20x₄ = 15

x₁ x₂ + 20x₃ + 4x₄ = 7

x₁ +x₂ 20x₃ 3x₄ = 5

File gauss_seidel.m

function xnew=gauss_seidel(A,b,xold)n=size(A)(1);At=A;xnew=xold;for k=1:nAt(k,k)=0;endfor k=1:nxnew(k)=(b(k)-At(k,:)*xnew)/A(k,k);end endfunction