Write a program to find the solution of linear system using Gaussian elimination with

partial pivoting followed by back substitution. The program should return the solution

x$,$ given matrix A and vector b$.$ The first few lines of your program could be

function $[$x$]=$gausselim$($A$,$b$)$

$[$n n$]=$size$($A$)$; r$=(1$:n$)$;

The code above sets n equal to the dimension of the matrix and initializes the pivot

vector r$.$ Use your program to solve the following system

$$

$$

$$

x$13$x$2+7$x$3=2,$

$2$x$1+4$x$23$x$3=1,$

$3$x$1+7$x$2+2$x$3=3$