% Banded fractional function
function val = FlTb(bands,T) % bands = [0,l1,l2,...,Inf] (um), T = K
  T=T(:);                                      % column vector
  val = diff(FlT(bsxfun(@times,T,bands)),1,2); 
end
