
 
 c = ones(5,1);
 a = -4*ones(5,1);
 b = ones(5,1);
 y = ones(5,1);
 my_tdma(c,a,b,y)
 
 A = diag(b(2:end), -1) + diag(a) + diag(c(1:end-1),1);
 A\y