Oyster_phase_matrix_BMP function

Function that calculates phase matrix