T =[26.533 26.574 26.614 26.729 26.838 26.970 27.14 27.92 28.90 29.92 31.80 34.53 42.57];
xz = csvread('TempCurve_data_0403_2019.csv');
xy=xz(:,9); al=xy; alpha =al; alphac =200.0; y=alpha/alphac;
F = csvread('FreqMaster_Single_f_data_0412_2019.csv'); zf =F(:,9);
FF= csvread('FreqMaster_Squaree_ff_data_0412_2019.csv'); zff =FF(:,9);
f =zf.*10^6; ff=zff.*10^12;
fc = 1.118752*10^11*((T- 26.50)./299.660).^1.90625;
fs =f./fc;
x =log(fs);
plot(x, y,'-kh',...
'LineWidth',0.5,...
'MarkerEdgeColor','k ',...
'MarkerFaceColor','r ',...
'MarkerSize',7 )