clear; clf;
epsilon = 0.1;
r=linspace(0,1,100);
x = linspace(0,1,100);
y = 0.5*(1-tanh((0.5-x)/(sqrt(2)*epsilon)));
plot(x,y,'k-','linewidth',1.3); hold on
line([0 1],[0 0],'linewidth',1.5,'color','k')
plot(1,0,'k>','markerface','k','markersize',9)
line([0 0],[-0.1 1.2],'linewidth',1.5,'color','k')
plot(0,1.2,'k^','markerface','k','markersize',9)

line([0 1],[1 1],'color','k','linestyle','-')
line([0 0.35],[0.1 0.1],'color','k','linestyle',':','linewidth',1.5)
line([0 0.65],[0.9 0.9],'color','k','linestyle',':','linewidth',1.5)
line([0.35 0.35],[0 0.1],'color','k','linestyle',':','linewidth',1.5)
line([0.65 0.65],[0 0.9],'color','k','linestyle',':','linewidth',1.5)
line([0.35 0.65],[0.5 0.5],'color','k','linestyle','-','linewidth',0.5)
plot(0.64,0.5,'k>','markerface','k','markersize',6)
plot(0.36,0.5,'k<','markerface','k','markersize',6)
text('Interpreter', 'latex','String' ,'$L$' , 'Position' ,[0.45 0.55],'FontSize',23)
text('Interpreter', 'latex','String' ,'$x$' , 'Position' ,[1.04 0.0],'FontSize',23)
text('Interpreter', 'latex','String' ,'$x_\alpha$' , 'Position' ,[0.32 -0.1],'FontSize',23)
text('Interpreter', 'latex','String' ,'$x_\beta$' , 'Position' ,[0.62 -0.1],'FontSize',23)
text('Interpreter', 'latex','String' ,'$0$' , 'Position' ,[-0.07 -0.01],'FontSize',16)
text('Interpreter', 'latex','String' ,'$0.05$' , 'Position' ,[-0.12 0.1],'FontSize',16)
text('Interpreter', 'latex','String' ,'$0.95$' , 'Position' ,[-0.12 0.9],'FontSize',16)
text('Interpreter', 'latex','String' ,'$1$' , 'Position' ,[-0.07 1.0],'FontSize',16)
text('Interpreter', 'latex','String' ,'$\phi^\infty$' , 'Position' ,[-0.12 1.12],'FontSize',23)
axis off
set(gcf,'PaperPositionMode','auto')
print('-deps','Schematic44.eps')


