Need this coded in Matlab:

Program Gaussian elimination with partial pivoting and back substitution for multiple systems. Your function should take in a coefficient matrix A, and a series of input vectors b as a matrix B. Your function should return a series of vector x as a matrix X, the final upper matrix U and final form of the input vector matrix, which we call C. That is, your code should solve AX=B. Your code should also return the appropriate error message The first line of your function should look like: function (X,U,C)= gaussMultiple Systems (A,B) Program Gaussian elimination with partial pivoting and back substitution for multiple systems. Your function should take in a coefficient matrix A, and a series of input vectors b as a matrix B. Your function should return a series of vector x as a matrix X, the final upper matrix U and final form of the input vector matrix, which we call C. That is, your code should solve AX=B. Your code should also return the appropriate error message The first line of your function should look like: function (X,U,C)= gaussMultiple Systems (A,B)