

%
% draw a biomembrane
%

t = 0:.01:1;

sx = t;
sy = sin(4*pi*t);

mx = cos(2*pi*t);
my = sin(2*pi*t);

s = 1/sqrt(2);
dx = [s 6-s];
dy = dx;

lx = [0 4];
ly = [0 0];

for h=1:3,

    plot(mx,my+(h-1)*6)
    hold on
    plot(lx+1,ly+(h-1)*6)
    plot(mx+6,my+(h-1)*6)
    plot(lx+7,ly+(h-1)*6)
    plot(mx+12,my+(h-1)*6)

end

for h=1:2,

    plot(ly,lx+1+(h-1)*6)
    plot(ly+6,lx+1+(h-1)*6)
    plot(ly+12,lx+1+(h-1)*6)

    plot(dx,dy+(h-1)*6)
    plot(6-dx,dy+(h-1)*6)
    plot(dx+6,dy+(h-1)*6)
    plot(-dx+12,dy+(h-1)*6)

end

for j=1:20,

    [x,y] = ginput(1);
    num = input(' Bar number: ');
    text(x,y,num2str(num),'fontsize',16)

end

text(-.3,-.2,'1','fontsize',16)
text(6-.3,-.2,'2','fontsize',16)
text(12-.3,-.2,'3','fontsize',16)
 
text(-.3,6-.2,'4','fontsize',16)
text(6-.3,6-.2,'5','fontsize',16)
text(12-.3,6-.2,'6','fontsize',16)

text(-.3,12-.2,'7','fontsize',16)
text(6-.3,12-.2,'8','fontsize',16)
text(12-.3,12-.2,'9','fontsize',16)

axis equal
grid
%axis off
hold off

