%% ondessonores.m
% Periode spatiale et pas de la grille
Lx= 1.; k= 2*pi/Lx; dX= 0.1; dY= 0.1;

% Domaine physique
pr= [-dX, 3*Lx + dX, -dY, 0.6*Lx + dY];

% Grille
X= (0 : dX : 3*Lx); Y= ( 0 : dY : 0.6*Lx); [Xp,Yp]= meshgrid(X,Y);

% Periode temporelle, frequence et pas de temps
T= 2.; omega= 2*pi/T; dt= T/20;

% Amplitude du deplacement
U= 0.06;

% Animation
for t= 0:dt:8*T
    xp= Xp - U*sin(k*Xp - omega*t) ; yp= Yp;
    plot(xp,yp,'ko','Markersize',4,'MarkerFaceColor','k');    
    axis(pr); xlabel('x'); ylabel('y');
    pause(0.05);
end
