function [x] = my_tridiag_solve_03_04(a,b,c,f) %risolve il sistema tridiagonale tridiag(b,a,c)x=f %mediante il calcolo della fattorizzazione LU con %l'algoritmo di Thomas [l,u,z]=my_tridiag_lu_03_04(a,b,c); %%1)Risolvere il sistema Ly=f, L=L(l,1) %[y]=my_lower_03_04[L,f] % DA FARE [y]=my_lower_sparse_03_04[l,f] %%2) Risolvere il sistema Ux=y, U=U(u,z) %[x]=my_upper_03_04[U,y] % DA FARE [x]=my_upper_sparse_03_04[u,z,y] end