def _calculate_Phi_x(self):
self.Phi_x = scipy.ones(N_CODON, dtype='<STR_LIT:float>')<EOL>for codon in range(N_CODON):<EOL><INDENT>for pos in range(<NUM_LIT:3>):<EOL><INDENT>self.Phi_x[codon] *= self.phi[pos][CODON_NT_INDEX[pos][codon]]<EOL><DEDENT><DEDENT>self.Phi_x /= self.Phi_x.sum()<EOL>
Calculate `Phi_x` (stationary state) from `phi`.