Using Octave or Matlab. Define a 100x100 square matrix A and the column vector b by:

Aij = Iij + 1/(-j+1), bi=1+2/i, 1<=i,j<=100

where Iij is the 100x100 identity matrix(i.e. 1 on the main diagonal and 0 everywhere else). Solve Ax=b for x using both Gauss-Seidel and A\b construct.Do not give the whole vector x in you output, but only x2,x50 and x99