filename1 = 'ext_GC1_17deg_MMI2_a.xlsx';
filename2 = 'ext_GC1_17deg_MMI2_b.xlsx';
filename3 = 'ext_GC1_17deg_MMI2_c.xlsx';
filename4 = 'ext_GC4_19deg_MMI2_a.xlsx';
filename5 = 'ext_GC4_19deg_MMI2_b.xlsx';
filename6 = 'ext_GC4_19deg_MMI2_c.xlsx';
filename7 = 'extext_GC1_17deg_MMI2_a.xlsx';
filename8 = 'extext_GC1_17deg_MMI2_b.xlsx';
filename9 = 'extext_GC1_17deg_MMI2_c.xlsx';
% filename10 = 'extext_GC4_19deg_MMI2_a.xlsx';
filename11 = 'extext_GC4_19deg_MMI2_b.xlsx';
filename12 = 'extext_GC4_19deg_MMI2_c.xlsx';
filename13 = 'taper_only_a.xlsx';
filename14 = 'taper_only_b.xlsx';
filename15 = 'taper_only_c.xlsx';
sheet1 = 1;
sheet2 = 2;
xlRange = 'A2:B10003';

M1 = xlsread(filename1, sheet2, xlRange);
M2 = xlsread(filename2, sheet2, xlRange);
M3 = xlsread(filename3, sheet2, xlRange);
M4 = xlsread(filename4, sheet2, xlRange);
M5 = xlsread(filename5, sheet2, xlRange);
M6 = xlsread(filename6, sheet2, xlRange);
M7 = xlsread(filename7, sheet2, xlRange);
M8 = xlsread(filename8, sheet2, xlRange);
M9 = xlsread(filename9, sheet2, xlRange);
% M10 = xlsread(filename10, sheet2, xlRange);
M11 = xlsread(filename11, sheet2, xlRange);
M12 = xlsread(filename12, sheet2, xlRange);
%{
M13 = xlsread(filename13, sheet2, xlRange);
M14 = xlsread(filename14, sheet2, xlRange);
M15 = xlsread(filename15, sheet2, xlRange);
%}

% M13 = cat(1, M1(1,2), M2(1,2), M3(1,2), M4(1,2), M5(1,2), M6(1,2), M7(1,2), M8(1,2), M9(1,2), M10(1,2), M11(1,2), M12(1,2));

%{
plot(X1,Y1,'g');
hold on

plot(X2,Y2,'r');
hold on

plot(X3,Y3,'k');
hold on

plot(X4,Y4,'m');
ylabel ('Total loss (dB)')
hold on
    
yyaxis right
    %ylim ([5 15])
    plot(c,b, 'b');
    ylabel ('dB/mm');
    hold on
    
title('SiWG023-C3: Horizontal-slot, strip-waveguide loss , GC design: 1, devices: c');
    xlabel('Wavelength (nm)');
    ylabel('Loss (dB/mm)');
    legend('WG length = 250um','WG length = 500um', 'WG length = 1000um', 'WG length = 2000um', 'Loss per mm (dB/mm)');
    hold on

    
disp(b);
disp((X1(5)));
disp(X2(5));
disp(X3(5));
%}