Nel caso di matrici tridiagonali il metodo più efficiente è quello di Thomas:
due passate una in avanti e una all'indietro e hai risolto.
il linguaggio più adatto è il Fortran poi trovi i listati (si fa per dire: < 10 righe) già bell'e pronti, oltre che in Fortran, anche in Pascal, Basic, C, forse anche Java