Gaussian Elimination help in getting upper triangular matrix

 programming in MATLAB that could work with any size of N*N 2D matrix where N>1. The values for the coefficients and dependent values (the right hand side column vector) should be represented as matrices in MATLAB. You should write a M-file (script or function) that computes and displays the upper triangular matrix. Your code should then use the upper triangular matrix to solve the set of linear simultaneous equations. You are given some test equations below. You should be able to demonstrate that you code can solve these equations.  
a +2b+3c +2d +e = 24