This scr ipt will attempt to find the band gaps in a phononic crystal by carrying out a Plane Wave Expansion of the various quantities using square array method. There are 2 programs, one for short and the other one is for acoustic long