x = -1:.01:1;

c1 = 0.5*(1+2*pi^3*x.^2).*exp(-pi^3*x.^2);
c2 = exp(-pi^3*x.^2/2);
c3 = 0.2*sin(2*pi*x);
c4 = 0.2*cos(2*pi*x);

clf
plot(x,c1,'r--','linewidth',2);hold on 
plot(x,c2,'linewidth',2);
plot(x,c3,'-.',x,c4','-.'); 
xlabel('WAVELENGTH')     
axis([-1.0000    1.0000   -0.2000    1.1])
h1 = gca;
set(h1, 'Box','off')
set(h1, 'YColor','black')
