Added the function delta_tign_calculation that calculates matrix of delta(t)