compute_shapley function

Compute shapley values