lambda         = 633e-9;
z              = 50e-3; % distance between aperture and image plane
yrange         = 1e-3;
N              = 400;

yvec           = linspace(-yrange, yrange,N);

v_vec          = yvec*sqrt(2/lambda/z);

[fs,fc]        = myfresnel(v_vec);

U              = 1/(1+i)*(fc+i*fs+0.5*(1+i));

Irradiance     = abs(U).^2;

figure
plot(yvec*1e3,Irradiance)
title(['Fresnel diffraction observed ' num2str(z*1e3) ' mm from an edge'])
xlabel('Transverse distance from shadow boundary in mm')
ylabel('Irradiance')
grid