find_deltak function

Compute δk\delta_k