Hej,
od paru godzin próbuję napisać program w C++ rozwiązujący układ równań z wykorzystaniem dekompozycji LU oraz metod podstawiania w przód i wstecz. O ile pierwsza część nie sprawiła problemu i zakodowałem dekompozycję LU, to utknąłem na drugiej części. W Internecie znalazłem parę implementacji, w tym również C++, ale nie do końca je rozumiem a adaptacja ich do mojego kodu na głupiego nie działa poprawnie. Przyklejam mój kod: https://pastebin.com/squpVkzN
Funkcja forward_subsitution
to moja próba implementacji podstawiania w przód, oczywiście nieudana
edit: chętnie przyjmę również review mojego kodu, na codzień nie koduje w C++