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