Update PE_perfect_xtal.m