clc
clear all
close all

lambda=633e-9;            % longueur d'onde
z=1e-3;                   % distance à parcourir
N=1024;                 % longueur du vecteur
X=500e-6;                   % valeur maximale de la coordonnée latérale

x1=(-N/2+1:N/2)/N*X;      % axe des x de -0.1mm à 0.1mm
dx1=x1(2)-x1(1);          % pas d'échantillonnage dans le plan z=0

% profil gaussien dans le plan z=0
w=50e-6;
y1=exp(-(x1/w).^2);
plot(x1,y1)

sy1=fftshift(y1);
[sy2, dx2]=Fresnelprog1D_TFD_unique(sy1,dx1,lambda,z);
y2=fftshift(sy1);

x2=(-N/2+1:N/2)*dx2;
hold on
plot(x2,abs(y2),'r-')
grid on