tridiagonal code snippets in matlab