% Schematic21.m %
clear all;clf;
for i=1:40
    points(i,1)=1+0.5*cos(i/(39)*2*pi);
    points(i,2)=1+0.5*sin(i/(39)*2*pi);
end
plot(points(:,1),points(:,2),'k-','linewidth',1.5)
values=points';
patch(points(:,1),points(:,2),[0.8627 0.8627 0.8627])
hold on
text(1.35,0.15, 'fluid 1','FontSize',20)
text(0.75,1.00,'fluid 2','FontSize',20)
axis ([0 2 0 2])
N=51;
for i=1:N
    aline(i,1)=1.0;
    aline(i,2)=2/(N-1)*(i-1);
    bline(i,1)=(2/(N-1)*(i-1)*0.70+0.30);
    bline(i,2)=2/(N-1)*(i-1);
end
annotation('arrow',[0.514285714285714 0.692857142857143],[0.840476190476191 0.838095238095238],'HeadStyle','plain');
annotation('arrow',[0.514285714285714 0.651785714285714],[0.7833333333334 0.780952380952381],'HeadStyle','plain');
annotation('arrow',[0.514285714285714 0.608928571428571],[0.721428571428572 0.719047619047619],'HeadStyle','plain');
annotation('arrow',[0.521428571428571 0.421428571428571],[0.314285714285714 0.311904761904762],'HeadStyle','plain');
annotation('arrow',[0.519642857142857 0.380357142857143],[0.261904761904762 0.261904761904762],'HeadStyle','plain');
annotation('arrow',[0.517857142857143 0.333928571428571],[0.201380952380953 0.202380952380953],'HeadStyle','plain');
t = 0.15;
hold on
plot(values(1,:),values(2,:),'k','linew',1.5);
[m,n]=size(values);
line([2.0 2.0],[ 2 0 ],'color','k','linewidth',1);
line([0 2],[ 2 2 ],'color','k','linewidth',1);
text('Interpreter','latex','String','$-U$','position',[0.7 -0.1],'FontSize',18)
text('Interpreter','latex','String','$U$','position',[1.2 2.1],'FontSize',18)
annotation('arrow',[0.248214285714286 0.851785714285714],[0.161904761904762 0.164285714285714],'HeadStyle','plain');
annotation('arrow',[0.248214285714286 0.248214285714286],[0.164285714285714 0.954761904761905],'HeadStyle','plain');
text('Interpreter','latex','String','$x$','position',[2.0+t -0.1],'FontSize',18)
text('Interpreter','latex','String','$y$','position',[-0.15 2.0+t],'FontSize',18)
text('Interpreter','latex','String','$H$','position',[1.9 -0.1],'FontSize',18)
text('Interpreter','latex','String','$H$','position',[-0.18 1.95],'FontSize',18)
text('Interpreter','latex','String','$0$','position',[ -0.175 -0.1 ],'FontSize',18)
axis image;axis off;box on;axis([-t 2+t -t 2.0+t]);
set(gcf,'PaperPositionMode','auto')