Need help on a matlab exercise

Writing a function that implements the Thomas algorithm to solve tridiagonal linear systems of equations of the form

With function name

function x = tridiag_solver(e,a,c,d)

Hint: We can debug the output of your functions backward sub(U,b) and tridiag solver(e,a,c,b) by comparing it with the output of the Matlab/Octave built-in linear solver linsolve(A,b) (see the documentation), for suitable prototype upper-triangular and tri-diagonal matrices, respectively.

a1 C1 0 b1 b2 T1 e1 a2 C2 T2 2 Tn-1 1 n- an-1 Cn-1 2 0.0 en-1 an J bn Tn TL Input: a a1 a2 anl, Output: x: solution to the linear system (2) (column vector)