In Matlab...

Create a function "Lastname_LUCrout.m" to compute the LU factorization of a given matrix. And solve the given system of equations. Your inputs will be the coefficient matrix A and the column of constants b. Your output is the x column of solutions Your solver has three parts a. LU factorization of A b. Forward solve for yin Ly = b c. Backward solve for x in Ux = y Write a sub-function for each part of the solver and call it in the main function Sample Code: function [x] = Lastname_LUCrout (A, b) % Summary of this function goes here % Detailed explanation goes here %Call the factorization sub-function [L U] = factorization(A) % Call the Foward solver for Ly = b [y] = forwardsolve(L, b) %Call the Backward Solver for Ux = y [x] = backwardsolve(U, y) end function [L U] = factorization (A) %Summary of this function goes here %Detailed explanation goes here %Steps to factorize the Matrix goes here end function [y] = forwardsolve (L, b) % Summary of this function goes here %Detailed explanation goes here %Steps to factorize the Matrix goes here end function [x] = backwardsolve (U, y) %Summary of this function goes here %Detailed explanation goes here