% Use this program to reproduce Fig 3.7 from text
clear all
close all
n = 0;
 for x = 0:.05:4
     n = n+1;
     sx(n) = quadl('fresnels',.0,x);
     cx(n) = quadl('fresnelc',.0,x);
end
 plot(cx)
 x=0:.05:4; 
 plot(x,cx,'k',x,sx,'k--')
 grid
 xlabel ('x')
 ylabel ('Fresnel integrals: C(x); S(x)')
 legend('C(x)','S(x)')