Create a function gaussian elimination that performs the gaussian elimination of linear system of the form Ax = b. The function should return the associated upper triangular matrix U and the modified right-hand side f. The function header should look something like:

function [U,f] = gaussian_elimination (A,b)

How would you write the following pseudocode shown below in MATLAB? For example, how would you write the code to find p and to interchange rows?

Gaussian Elimination with Backward Substitution To solve the n n linear system INPUT number of unknowns and equations n; au i Sn and 1 sj sn+1 ed matrix A = [aij], where 1 OUTPUT solution x1, x2, xn or message that the linear system has no unique solu- tion. Step 1 For i = 1, . . . , n-1 do Steps 2-4. (Elimination process.) Step 2 Let p be the smallest integer with 1