Create PE_perfect_xtal.m