Create plotPsa.m