function [x, l, u] = thomas(a,b,c,z) % Opis: % thomas resi sistem s tridiagonalno matriko z LU razcepom brez % pivotiranja. % % Definicija: % [x,l,u] = thomas(a,b,c,z) % % Vhodni podatki: % a glavna diagonala matrike, % b diagonala nad glavno diagonalo matrike, % c diagonala pod glavno diagonalo matrike, % z stolpec, ki predstavlja desno stran sistema. % % Izhodni podatki: % x resitev sistema A x = z, kjer je matrika A dolocena z diagonalami % a, b, c, % l diagonala pod glavno diagonalo enic matrike L v LU razcepu A % (ostali elementi L so enaki 0), % u glavna diagonala matrike U v LU razcepu matrike A (diagonala nad % glavno diagonalo je enaka b, vsi ostali elementi U so enaki 0). end