phi.calculate function

Update pheromone levels for each decision node