phases function

Compute Phases