cornu_spiral_fresnel_integrals.m

% The Cornu spiral
%Principles of Optics, MaxBorn & Emil Wolf,7th edition, Cambridge University Press
%Page No. 481, Fig 8.37
% T.Venkata Rama Krishna
% E-mail: [e-mail]
clc;
clear all;
echo on;
n = 1000;
x = zeros(1,n); 
y = x;
t = linspace(0,4*pi,n);
for i=1:n-1
x(i) = quadl(inline('cos(x.^2)'),t(i),t(i+1),1e-6);
y(i) = quadl(inline('sin(x.^2)'),t(i),t(i+1),1e-6);
end
x = cumsum(x);  % Fresnel cosine integral
y = cumsum(y);  % Fresnel sine integral
plot([-x(end:-1:1) 0 x], [-y(end:-1:1) 0 y])
grid;
axis equal