x2 = linspace(-2*pi,2*pi,2000);



I = sinc(x2).^2.*sin(x2).^2;

plot(x2,I);

xlabel('x2 in arbitrary units');

ylabel('I in arbitrary units');

grid on;

title('Problem 4.3: Fraunhofer diffraction from a phase-sandwich');