The code below takes the matrix A and creates the upper triangular matrix U.

We have modified the code so that it does banded Gauss elimination p = q when the function is given the value of p as an input argument.

Determine the entires for B through H and L through N.

1 2 3 4 5 6 7 8 for k=1:n-1 end for i=k+1:n m=A (i, k)/A (k, k); for j=k+1:n end end A(i, j) A (i, j)-m*A (k, j);