%%This script will try to propagate an acoustic wave through a two-d
%%phononic crystal

I=400;              J=400;
i=1:I;              j=1:J;
im=2:I-1;           jm=2:J-1;

p=zeros(I,J);           pho=zeros(I,J);
vx=zeros(I,J);          vy=zeros(I,J);