function [Eout] = Fresnelprog2D_spectre_angulaire (Ein,deltax,deltay,lambda,z)

% Calcul de la propagation de Fresnel selon l'algorithme de la propagation
% du spectre angulaire

[N,M]=size(Ein);
p=[(0:(M/2)-1) (-M/2:-1)];
q=[(0:(N/2)-1) (-N/2:-1)]';
[pp,qq]=meshgrid(p,q);

fx=pp/(M*deltax);    % définition des axes des fréquences
fy=qq/(N*deltax);    %
H=exp(2*1i*pi*z/lambda*(1-0.5*lambda^2*(fx.^2+fy.^2)));
Eout=ifft2(fft2(Ein).*H);