% Boundary x = 2 from semi-infinite strip Slide 32

xx = linspace(0,10,1000);
F = @(w,c) 4*(1-cos(5*w)).*sin(c*w)./(pi*w);
y = @(x) integral(@(w) F(w,x),0,50);
Y = arrayfun(y,xx);
plot(xx,Y,'b-','LineWidth',1.5);
hold on;
plot([0 5],[2 2],'k-','LineWidth',1.0);
plot([5 10],[0 0],'k-','LineWidth',1.0);
grid;

xlim([0,10]);
ylim([-1,3]);
fontlabs = 'Times New Roman';  
xlabel('$x$','FontSize',16,'FontName',fontlabs,'interpreter','latex');  
ylabel('$f(x)$','FontSize',16,'FontName',fontlabs,'interpreter','latex'); 
mytitle = '$ a= 100$';
title(mytitle,'FontSize',16,'FontName',fontlabs,'interpreter','latex'); 
set(gca,'FontSize',16);        

print -depsc test_eg.eps