filename1 = 'fin200_period400_250_b.xlsx';
%{
filename2 = 'GC1_250_b.xlsx';
filename3 = 'GC1_250_c.xlsx';
filename4 = 'GC1_500_a.xlsx';
filename5 = 'GC1_500_b.xlsx';
filename6 = 'GC1_500_c.xlsx';
filename7 = 'GC1_1000_a.xlsx';
filename8 = 'GC1_1000_b.xlsx';
filename9 = 'GC1_1000_c.xlsx';
filename10 = 'GC1_2000_a.xlsx';
filename11 = 'GC1_2000_b.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);
%}

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