function []=figure2c()
close all
figure()
hold on
%% Data parameters
ceiling_T_VFT=0.0907;
ceiling_T_A=0.0738;
T_g=mean([ceiling_T_VFT,ceiling_T_A]);

%% Loading data
x_grain=0.03*(T_g);
Ts_fit=0:x_grain:T_g;

% CV
data_name=sprintf('./CV.dat');
fid = fopen(data_name,'r'); numbers=fscanf(fid, '%f'); fclose('all');
Ndata=round(size(numbers,1)/2);
Ts=numbers(2*(1:Ndata)-1);
CVs=numbers(2*(1:Ndata));

%% Fitting
where=logical(Ts<T_g);
Ts_trun=Ts(where);
CVs_trun=CVs(where);
p=polyfit(Ts_trun,CVs_trun,1);
b=p(1); a=p(2);
CV_fits=a+b*Ts_fit;

%% Plotting
scatter(Ts,CVs,100,0.88*[255,165,0]/255,'Marker','<')
plot(Ts_fit,CV_fits,'Color',0.88*[255,165,0]/255,'LineStyle',':')
plot([0,0.15],[1,1],'Color',[0,0,0],'LineStyle',':')
xlim([0,0.15])
ylim([0,2])
xlabel('$T$','Interpreter','latex');
ylabel('$c_{V}$','Interpreter','latex')
print(gcf, '-dpdf', 'figure2c.pdf');
end