Muss NICHT zu vollbesetzten Matrizen führen! Siehe Tau-Tschebyscheff-Methode: Obere Dreiecksmatrix, die ganz einfach in Tridiagonal-Gestalt gebracht werden kann und über eine LU-Zerlegung SEHR effizient gelöst werden kann. Solange die Zeitschrittweite konstant bleibt bleiben diese Matrizen konstant, wodurch die LU-Zerlegung enorme Vorteile bietet.

Sehr ungenauer Beitrag - werde ich demnächst überarbeiten...