model_setup
a: sub-diagonal (len n-1), A[i, i-1] b: main diagonal (len n), A[i, i] c: super-diagonal (len n-1), A[i, i+1] d: RHS (len n)