function [x] = my_tridiag_solve(a,b,c,f)
%risolvere sistema tridiagonale mediante il calcolo 
%della fattorizzazione lu
% tridiag(a,b,c)x=f
[l,u,z]=my_tridiag_lu_31_03(a,b,c);
%%risolvere Ly=f%%%%%%%%%%%%%%%
%%modificare my_triangular_lower 
%%%per accettare in ingresso L codificata sparsa
%%%risolvere Ux=y%%%%%%%%%%%%%%%%
%%%modificare my_triangular_upper 
%%%%per accettare in ingresso U codificata sparsa
%%%%Analizzare il costo computazionale
end