function plotit_trap b = 30; n1 = 1.15; input = 0; n2 = 1.0:0.01:1.3; s = size(n2); a1 = 45; a2 = [45 0 -45]; for i = 1:3, for t = 1:s(2), r(t) = n2(t); d(t,i) = deviation(input,a1,a2(i),b,n1,n2(t)); end; end; plot(r,d) line([1.0 1.40],[0 0],'LineStyle','-.','Color','k'); line([1.15 1.15],[-30 10],'LineStyle','-.','Color','k'); text(1.01,17,'First tapered plate index of refraction = 1.15',... 'FontSize',10,'FontWeight','demi'); text(1.01,13,'First block angle = 45 degrees',... 'FontSize',10,'FontWeight','demi'); text(1.01,15,'Interface angle = 30 degrees',... 'FontSize',10,'FontWeight','demi'); text(1.2,-15,'---- Second block angle = + 45 degrees','color','b',... 'FontSize',10); text(1.2,-17.5,'---- Second block angle = 0 degrees','color',[0 .65 0],... 'FontSize',10); text(1.2,-20,'---- Second block angle = - 45 degrees','color','r',... 'FontSize',10); xlabel('Second tapered plate index of refraction','color','b',... 'FontSize',12,'FontWeight','demi'); ylabel('Angular deviation from horizontal','color','b',... 'FontSize',12,'FontWeight','demi'); title('Angular Deviation Provided by "Light Guiding Device"','color','r',... 'FontSize',14,'FontWeight','demi');