clc;close all; 

set(0,'DefaultFigureColor',[1 1 1])
v_simulation=[
0.3512
0.365
0.2273
0.1747];
v_experiment=[
0.23
0.282
0.216
0.147];
T_simulation=[462,492,510,537.83];
T_experiment=[482,521,524,527.8];
composition=[100 75 50 25];
yyaxis left

plot(composition,v_simulation,'o','MarkerSize',20,'MarkerEdgeColor','blue','LineWidth',2.5)
hold on
error_v=[0.041658,0.08120,0.05162,0.03031];
ylabel('Propagation Speed (mm/min)');
xlabel('Cellulose Content (%)');

h=errorbar(composition,v_experiment,error_v,'d','MarkerEdgeColor','blue','LineWidth',2.5);
set(h, 'MarkerSize', 28);

h.CapSize = 12;
%[lgd, icons, plots, txt] = legend('show');
ylim([0.0 0.4])
yyaxis right
% scatter(composition,T_simulation,200,'r','o');
% scatter(composition,T_simulation,200,'r','o');
%scatter(composition,T_simulation,200,'r','o');
ylim([455 600])
plot(composition,T_simulation,'o','Linewidth',5,'MarkerSize',20,'MarkerFaceColor','red','MarkerEdgeColor','red')
hold on
plot(composition,T_experiment,'d','Linewidth',5,'MarkerSize',20,'MarkerFaceColor','red','MarkerEdgeColor','red')

%plot(cell2mat(composition(1:4)),cell2mat(T_simulation(1:4)),'-s')

ylabel('Temperature (°C)')
xlabel('Cellulose Content (%)');
set(gca,'FontSize',25)

xlim([17 108])
set(gca,'XTick',(25:25:100))
legend('Speed from simulation','Speed from experiment','Temperature from simulation','Temperature from experiment');

hFig = figure(1);
set(hFig, 'Position', [0 0 1000 800])

export_fig validation.pdf