This function computes the qsi matrix for specified indices and time points.
find_qsi(is, tau_no_Inf)
is
Indices of nodes.
tau_no_Inf
Unique time points excluding infinity.
qsi matrix.