$*1.$ Write a MATLAB function which computes the solution of $Ax=g$ using a

normalized Cholesky decomposition, when $A$ is a tridiagonal, Toeplitz, positive

definite matrix of size $nn,$ with $a$ on the diagonal and $b$ on the superdiagonal and

subdiagonal. The first line of your code should be: function $x=$ trid$($a$,$b$,$g$).$ Show

how you tested your code for accuracy.

Coding instructions: If an array grows inside a "for" loop, pre$-$allocate all necessary

memory. Use only two "for" loops and do not use any logical tests. Do not create or

store any two$-$dimensional arrays and do not use Matlab's chol function. Store only

what is necessary and reuse variable names by overwriting.