clear; clf;
epsilon = 0.03;
r=linspace(0,1,100);
x = linspace(0,1,100);
y = tanh((0.25-abs(r-0.55))/(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],[-1.2 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 1],[1 1],'color','k','linestyle','--')
text('Interpreter', 'latex','String' ,'$\phi(r,t)$' , 'Position' ,[-0.07 1.4],'FontSize',16)
text('Interpreter', 'latex','String' ,'$r$' , 'Position' ,[1.04 0.0],'FontSize',16)
set(gca,'fontsize',13)
plot(0.3,0,'ks','markersize',12,'markerface','k')
plot(0.8,0,'ks','markersize',12,'markerface','k')
text('Interpreter', 'latex','String' ,'$R_2$' , 'Position' ,[0.3 -0.15],'FontSize',16)
text('Interpreter', 'latex','String' ,'$R_1$' , 'Position' ,[0.81 -0.15],'FontSize',16)
text('Interpreter', 'latex','String' ,'$0$' , 'Position' ,[-0.07 -0.01],'FontSize',16)
text('Interpreter', 'latex','String' ,'$1$' , 'Position' ,[-0.07 1.0],'FontSize',16)
text('Interpreter', 'latex','String' ,'$-1$' , 'Position' ,[-0.15 -1.0],'FontSize',16)
text('Interpreter', 'latex','String' ,'$1$' , 'Position' ,[0.975 -0.16],'FontSize',16)
line([0.69 0.92],[-0.25 -0.25],'color','k','linewidth',1.3)
line([0.69 0.92],[0.25 0.25],'color','k','linewidth',1.3)
line([0.69 0.69],[-0.25 0.25],'color','k','linewidth',1.3)
line([0.92 0.92],[-0.25 0.25],'color','k','linewidth',1.3)
axis off
axis([0 2 -1.5 1.5])
axes('position',[0.6 0.25 0.28 0.55])
r1=linspace(0,1,1000);
x1 = linspace(0,1,1000);
y1 = tanh((0.25-abs(r1-0.55))/(sqrt(2)*epsilon));
plot(x1,y1,'k-','linewidth',1.3); hold on
plot(x1(794),y1(794),'ko','linewidth',1.3)
plot(x1(807),y1(807),'ko','linewidth',1.3)
line([0 1],[0 0],'color','k','linewidth',1.5)
plot(0.8,0,'ks','markersize',12,'markerface','k')
text('Interpreter', 'latex','String' ,'$R_1$' , 'Position' ,[0.798 -0.05],'FontSize',16)
line([x1(794) x1(794)],[0 y1(794)],'color','k','linestyle','--')
line([x1(807) x1(807)],[0 y1(807)],'color','k','linestyle','--')
text('Interpreter', 'latex','String' ,'$r_{m_1}$' , 'Position' ,[0.792 -0.04],'FontSize',16)
text('Interpreter', 'latex','String' ,'$r_{m_1 +1}$' , 'Position' ,[0.803 0.04],'FontSize',16)
text('Interpreter', 'latex','String' ,'$\phi (r_{m_1},t)$' , 'Position' ,[0.793 0.18],'FontSize',16)
text('Interpreter', 'latex','String' ,'$\phi (r_{m_1 +1},t)$' , 'Position' ,[0.7935 -0.16],'FontSize',16)
axis([0.79 0.81 -0.25 0.25])
set(gca,'fontsize',13)
axis off
line([0.79 0.81],[-0.25 -0.25],'color','k','linewidth',1.1)
line([0.79 0.81],[0.25 0.25],'color','k','linewidth',1.1)
line([0.79 0.79],[-0.25 0.25],'color','k','linewidth',1.1)
line([0.81 0.81],[-0.25 0.25],'color','k','linewidth',1.1)
annotation(figure(1),'line',[0.488 0.6],[0.45 0.25],'linewidth',1.1);
annotation(figure(1),'line',[0.488 0.601],[0.587 0.8],'linewidth',1.1);
print('-deps','Schematic43.eps')


