function ondezqq;
% Dominic Grenier décembre 2002

global d2r;

theta = 30; %degrés
freq = 1e6;
vp = 3e8;
alpha = 0.0012;

AZ = 42;
EL = 40;
theta =theta*d2r;
lambda = vp/freq;
beta = 2*pi/lambda;

xmax = 10^(round(log10(3*lambda)+0.5));
x = linspace(0,xmax,50);
z = x;
tT=0:0.1:1;

[X,Z] = meshgrid(x,z);
ZP = X*sin(theta) + Z*cos(theta);

M=moviein(length(tT));
for i=1:length(tT);
    E = exp(-alpha*ZP).*cos(2*pi*tT(i)-beta*ZP);
    figure(1);
    mesh(Z,X,E)
    caxis([-1 1]);
    axis([0,xmax,0,xmax,-1,1]);
    daspect([xmax xmax 5]);
    view(AZ,EL);
    M(:,i)=getframe;
end
movie(M,3,4)