Oyster_phase_matrix_TIF function

Function that calculates phase matrix