Shun-Lien-Chuang-Physics-of-Photonic-Devices-Wile_Shun lien Chua

clear clear figure clc i=1; for i=1:101; x(i)=(i-1)*0.01; E1(i) = 0.36 + 0.505*x(i) + 0.555*x(i)^2; a1(i) = x(i)*5.6533 + (1-x(i))*6.0584; E2(i) = 0.36 + 2.35 *x(i) + 0.24 *x(i)^2; a2(i) = x(i)*5.66 + (1-x(i))*6.0584; end % plot(x,E1); % TT = text(0.65,1.5,sprintf('Eg(GaAs)=%1.2feV',E1(101))); % set(TT,'fontsize',12) % TT = text(0.05,0.3,sprintf('Eg(InAs)=%1.2feV',E1(1))); % set(TT,'fontsize',12) % xlabel('mole fraction') % set(gca,'fontsize',12) % ylabel('Eg[eV] for In(1-x)Ga(x)As') % set(gca,'fontsize',12) % plot(x,E2); % TT = text(0.65,2.9,sprintf('Eg(AlAs)=%1.2feV',E2(101))); % set(TT,'fontsize',12) % TT = text(0.05,0.2,sprintf('Eg(InAs)=%1.2feV',E2(1))); % set(TT,'fontsize',12) % xlabel('mole fraction') % set(gca,'fontsize',12) % ylabel('Eg[eV] for Al(x)In(1-x)As') % set(gca,'fontsize',12) % plot(a1,E1); % hold on % plot(a1(1),E1(1),'O'); % plot(a1(101),E1(101),'O'); % hold off % TT = text(5.65,1.47,sprintf('a(GaAs)=%1.2fA',a1(101))); % set(TT,'fontsize',12) % TT = text(5.65,1.55,sprintf('Eg(GaAs)=%1.2feV',E1(101))); % set(TT,'fontsize',12) % TT = text(5.9,0.27,sprintf('a(InAs)=%1.2fA',a1(1))); % set(TT,'fontsize',12) % TT = text(5.9,0.35,sprintf('Eg(InAs)=%1.2feV',E1(1))); % set(TT,'fontsize',12) % xlabel('lattice constant[A]') % set(gca,'fontsize',12) % ylabel('Eg[eV] for In(1-x)Ga(x)As') % set(gca,'fontsize',12) % plot(a2,E2); % hold on % plot(a2(1),E2(1),'O'); % plot(a2(101),E2(101),'O'); % hold off % TT = text(5.7,2.7,sprintf('a(AlAs)=%1.2fA',a2(101))); % set(TT,'fontsize',12) % TT = text(5.7,2.85,sprintf('Eg(AlAs)=%1.2feV',E2(101))); % set(TT,'fontsize',12) % TT = text(5.9,0.15,sprintf('a(InAs)=%1.2fA',a2(1))); % set(TT,'fontsize',12) % TT = text(5.9,0.3,sprintf('Eg(InAs)=%1.2feV',E2(1))); % set(TT,'fontsize',12) % xlabel('lattice constant[A]') % set(gca,'fontsize',12) % ylabel('Eg[eV] for Al(x)In(1-x)As') % set(gca,'fontsize',12) plot(a1,E1,'--k',a2,E2,'k'); hold on plot(a1(1),E1(1),'O'); plot(a1(101),E1(101),'O'); plot(a2(1),E2(1),'O'); plot(a2(101),E2(101),'O'); hold off TT = text(5.65,1.55,sprintf('a(GaAs)=%1.2fA',a1(101))); set(TT,'fontsize',12) TT = text(5.65,1.7,sprintf('Eg(GaAs)=%1.2feV',E1(101))); set(TT,'fontsize',12) TT = text(5.9,0.20,sprintf('a(InAs)=%1.2fA',a1(1))); set(TT,'fontsize',12) TT = text(5.9,0.35,sprintf('Eg(InAs)=%1.2feV',E1(1))); set(TT,'fontsize',12) TT = text(5.7,2.7,sprintf('a(AlAs)=%1.2fA',a2(101))); set(TT,'fontsize',12) TT = text(5.7,2.85,sprintf('Eg(AlAs)=%1.2feV',E2(101))); set(TT,'fontsize',12) xlabel('lattice constant[A]') set(gca,'fontsize',12) ylabel('Eg[eV]') set(gca,'fontsize',12)