 X1 = [.45,0,0,0,0,0,0,0];
Y1 = [.35,.20,.15,.15,.15,.15,.10,.10];
Z1 = [20,30,40,50,60,70,80,90];
Z2 = [40,50,60,70,80,90];
Z3 = [30,40,50,60,70,80,90];
X2 = [.25,0,0,0,0,0,0,0];
Y2 = [.30,.26,.20,.21,.19,.19,.14,.09];
X3 = [.25,0,0,0,0,0,0,0];
Y3 = [.35,.30,.30,.25,.20,.20,.15,.10];
X4 = [.5,.2,.1,0,0,0,0];
Y4 = [0,.25,.18,.25,.25,.25,.20,.20];
X5 = [.7,.6,.5,.3,0,0];
Y5 = [0,0,.10,.20,.16,.16,.15,.15];


SUBPLOT(2,1,1), hold on
plot(Z1,X1,'--b','LineWidth',4,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',10)
plot(Z1,X2,'--y','LineWidth',4,'MarkerEdgeColor','k','MarkerFaceColor','y','MarkerSize',10)
plot(Z1,X3,'--g','LineWidth',4,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',10)
plot(Z3,X4,'--r','LineWidth',4,'MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',10)
plot(Z2,X5,'--k','LineWidth',4,'MarkerEdgeColor','k','MarkerFaceColor','c','MarkerSize',10)
plot(30,.5,'Xk','LineWidth',2,'MarkerSize',10);
plot(40,.7,'Xk','LineWidth',2,'MarkerSize',10);

title('Required Normal Force For Engagement');

ylabel('Force (Newtons)');

legend('Balsa','Rubber','Cork','Particle B.','Hardwood');

SUBPLOT(2,1,2), hold on
plot(Z1,Y1,'--b','LineWidth',3,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',10)
plot(Z1,Y2,'--y','LineWidth',3,'MarkerEdgeColor','k','MarkerFaceColor','y','MarkerSize',10)
plot(Z1,Y3,'--g','LineWidth',3,'MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',10)
plot(Z1,Y4,'--r','LineWidth',3,'MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',10)
plot(Z1,Y5,'--k','LineWidth',3,'MarkerEdgeColor','k','MarkerFaceColor','c','MarkerSize',10)

title('Disengagment (Adhesion) Normal Force');
xlabel('Pin Angle (Degrees)');
ylabel('Force (Newtons)');