% plot of four functions
N=200;
a=-0.5;
b=0.5;
z=linspace(a,b,N);

g = exp(-4*pi^2*z.^2);
f = (16*pi^2*z.^2+1) .* exp(-8*pi^2*z.^2-0.444);
h = 0.5/pi*cos(2*pi*z);
s = 0.5/pi*sin(2*pi*z);

horz = [-0.5 0.5];
verz = [0 0];

horf = [0 0];
verf = [-0.5 1.1];

plot(z,h,':',z,s,':',horz,verz,horf,verf);
hold on
hg=plot(z,g);
hf=plot(z,f,'.');
hold off;
set(hg,'linewidth',2);
set(hf,'linewidth',3);


set(gca,'box','off');
set(gca,'xlim',horz,'ylim',verf);
set(gca,'xcolor','black');
set(gca,'ycolor','black');
set(gca,'zcolor','black');
set(gca,'aspect',[nan 1]);

ht1=text(-0.5,-0.1,'-0.5','horizontal','center','vertical','baseline');
ht2=text(0.5,-0.1,'0.5','horizontal','center','vertical','baseline');
ht3=text(0.03,-0.1,'WAVELENGTH','fontsize',7,'vertical','baseline');
