Please solve using Wolfram Mathematica coding:

Suppose A is a 50 x 50 tridiagonal matrix with entries a_ii=0.5 on the diagonal, a_i,i-1= a_i,i+1=0.25 on the super diagonal and sub diagonal and zeros elsewhere. Take b to have i^th entry equal to 1/i. Starting with x₀=(1,1,.,1)

Write a program using the Gauss-Siedel method to the 50 x 50 tridiagonal matrix. Build a corresponding 100 x 100 tridiagonal matrix and apply Gauss-Seidel. How does the size change affect the relative normed error at 50, 500 and 5000 iterations? Use the solution from Gaussian elimination to calculate the relative normed error.