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

% Calcul de la propagation de Fresnel selon l'algorithme de 
% convolution / 3 TFD

N=length(Ein);
p=[(0:(N/2)-1) (-N/2:-1)];
x=p*deltax;
hz=exp(1i*pi/(lambda*z)*x.^2);
Eout=deltax*...
        exp(1i*2*pi/lambda*z)/sqrt(1i*lambda*z)*...
        ifft(fft(Ein).*fft(hz));