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

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

N=length(Ein);
p=[(0:(N/2)-1) (-N/2:-1)];
fx=p/(N*deltax);    % définition de l'axe des fréquences
H=exp(2*1i*pi*z/lambda*(1-0.5*lambda^2*fx.^2));
Eout=ifft(fft(Ein).*H);