filename1 = 'GC4_19deg_250_a.xlsx';
filename2 = 'GC4_19deg_250_b.xlsx';
filename3 = 'GC4_19deg_250_c.xlsx';
filename4 = 'GC4_19deg_500_a.xlsx';
filename5 = 'GC4_19deg_500_b.xlsx';
filename6 = 'GC4_19deg_500_c.xlsx';
filename7 = 'GC4_19deg_1000_a.xlsx';
filename8 = 'GC4_19deg_1000_b.xlsx';
filename9 = 'GC4_19deg_1000_c.xlsx';
filename10 = 'GC4_19deg_2000_a.xlsx';
filename11 = 'GC4_19deg_2000_b.xlsx';
filename12 = 'GC4_19deg_2000_c.xlsx';
%{
filename13 = 'taper_only_a.xlsx';
filename14 = 'taper_only_b.xlsx';
filename15 = 'taper_only_c.xlsx';
%}
sheet1 = 1;
sheet2 = 2;
xlRange = 'A2:B9990';

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);

M21 = 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));
%}