clear
clc

load('LSGO-DEA_areal_strain.mat')

ydata =     {          p3_S,   UoB_S(:, 1),   UoB_S(:, 2),   UoB_S(:, 3)};
xdata =     {          p3_V,         UoB_V,         UoB_V,         UoB_V};
lineSpecs = {          'o-',          's-',          's-',          's-'};
% colors =    { [0.2 0.5 1.0], [0.0 0.8 0.0], [0.0 0.6 0.0], [0.0 0.4 0.0]};
colors =    { [0.2 0.5 1.0], [0.9 0.2 0.0], [0.7 0.1 0.0], [0.5 0.0 0.0]};

fwidth = 500;
fheight = 500;

xmin = 4800;
xmax = 10200;

%% voltage-strain plot

f = figure('pos', [100 100 fwidth fheight]);
ax = axes;

hold on
  

for i = 1:4
    x = xdata{i};
    y = ydata{i};
    
    h = plot(x, y, lineSpecs{i});
    h.Color = colors{i};
    h.MarkerEdgeColor = colors{i};
    h.MarkerFaceColor = colors{i};
    h.LineWidth = 1.5;
    h.MarkerSize = 10;
end

axis([xmin xmax 0 9.5])

box on
% grid on

xlabel('\itVoltage\rm [V]')
ylabel('\itAreal strain\rm [%]')

lgndlbl = {'3-segment', 'UoB, 1st cycle', 'UoB, 2nd cycle', 'UoB, 3rd cycle'};
legend(lgndlbl, 'Location', 'northwest');
legend boxoff

ax.FontSize = 16;
ax.LineWidth = 1.5;
%ax.FontWeight = 'Bold';



