Including tridiagonal system solving algorithm (Thomas Algorithm) so we can solve for stencils including other future values
Tested tridiagonal solver, next to try crank-nicolson / backwards euler